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ABSTRACT 

It  is  well  known  that  natural  and  laboratory 
processes  can  alter  isotopic  abundances.  This  thesis 
is  an  initial  study  to  find  the  extent  to  which  tellurium 
isotope  ratios  may  be  altered. 

In  the  introduction,  a  brief  review  of  historical 
isotope  studies  is  given  and  some  of  the  natural  mech¬ 
anisms  which  alter  isotope  abundances  summarized.  Isotope 
fractionation  studies  of  sulphur  and  selenium  are  re¬ 
viewed  more  extensively  because  tellurium  is  similar  to 
sulphur  and  selenium. 

The  vibrational  frequencies  are  computed  for 

isotopic  compounds  of  tellurium.  From  these,  vibrational 

partition  function  ratios  are  calculated  and  these  ratios 

are  in  turn  used  to  calculate  equilibrium  constants  for 

iqn  122 

various  Te  -  Te  exchange  reactions  at  different 

temperatures.  Providing  that  the  exchange  can  be  effected, 

l^o  122 

variations  in  the  Te  J  /Te  ratio  in  equilibrium  iso¬ 
tope  reactions  of  up  to  2.6/  at  20°C  are  predicted.  Cal¬ 
culations  are  also  made  of  kinetic  isotope  effects  pre¬ 
dicting  values  as  high  as  1.0175  for  the  ratio  of  isotopic 


rate  constants. 


iv. 


Experimental  verification  of  some  of  the  theo¬ 
retical  predictions  was  obtained.  In  the  chemical 

reduction  of  potassium  tellurite  at  room  temperature, 

1?2  =  4- 

the  Te  0^  ion  was  found  to  be  reduced  1.0071-  0.0003  times 

130  = 

faster  than  the  reduction  rate  of  Te  0^ .  In  a 

122  = 

microbiological  reduction  of  tellurite  ion,  the  Te  0^ 

ion  was  found  to  be  reduced  1.0060  t  0.0005  times  faster 

130  = 

than  the  reduction  rate  of  Te  0^ . 

On  the  bases  of  the  theoretical  predictions,  the 
effects  noted  in  both  reduction  studies  are  related  to 
tellurium  oxygen  bond  breakage. 

Of  four  natural  samples,  three  were  identical  in 
122  130 

their  Te  /Te  ^  ratio;  whereas,  the  fourth  was  enriched 
by  0.4  percent  in  the  lighter  isotope. 


ACKNOWLEDGEMENTS 


The  author  wishes  to  express  his  sincere  thanks 
to  his  supervisor.  Dr.  H.  R.  Krouse,  for  suggesting  this 
project  and  for  his  inspiration,  guidance,  and  assistance 
in  the  course  of  this  investigation.  The  assistance  of 
various  members  of  the  Department  of  Computing  Science; 
the  technical  services  of  P.  Alexander,  glassblower;  and 
assistance  of  R.  McCready  of  the  Department  of  Micro¬ 
biology  in  obtaining  bacterial  cultures  and  medium  are 
also  greatly  appreciated. 

Finally,  I  wish  to  acknowledge  the  assistance 
of  my  wife,  Marilyn,  and  Mrs.  P.  Babet  in  the  preparation 
of  the  manuscript. 


TABLE  OF  CONTENTS 


vi . 


INTRODUCTION  Page  1 

A.  Equilibrium  Isotope  Effects  3 

B.  Kinetic  Isotope  Effects  4 

C.  Mechanisms  of  Isotope  Fractionation 

in  Nature  6 

1.  Physical  processes  6 

a.  Diffusion  6 

b.  Isotopic  vapour  pressures  7 

c.  Miscellaneous  9 

2.  Chemical  processes  9 

a.  Equilibrium  isotopic  exchange  9 

b.  Kinetic  isotope  effects  13 

3.  Biological  processes  13 

D.  The  Group  VIA  Elements  15 

1.  General  chemistry  15 

2.  Isotope  fractionation  studies  18 

a .  Sulphur  l8 

b.  Selenium  26 

3.  Tellurium  28 

a.  History  and  occurrence  28 

•  b.  Uses  30 

c.  Physical  properties  30 

d.  Chemical  properties  3^- 

e.  Absolute  abundances  of  the  isotopes 

of  tellurium  35 

E.  Advantages  of  looking  at  sulphur,,  selenium,, 

and  tellurium  isotopes  simultaneously  in 
natural  deposits  37 

THEORY  4l 

A.  Theory  of  Equilibrium  Isotope  Effects  42 


vii . 


B.  Evaluation  of  Partition  Function  Ratios  45 

C.  Kinetic  Isotope  Effects  49 

D.  Isotope  Fractionation  During  First  Order 

Competing  Reactions  56 

THEORETICAL  CALCULATIONS  59 

A.  Calculation  of  Equilibrium  Exchange 

Constants  59 

B.  Calculation  of  Ratios  of  Isotopic  Rate 

Constants  64 

EXPERIMENTAL 

A.  Extraction  of  Tellurium  from  Natural 

Samples  89 

B.  Chemical  Reduction  of  KoTe0o  ^ 

2  3  7-l 

C.  Reduction  of  KpTeOp  by  Scopulariopsis 

brevicaulis^  ^  73 

1.  Reduction  of  KpTeO^  on  liquid  medium 

2.  Extraction  of  elemental  tellurium 

reduced  from  KpTeOo  in  the  biological 
investigation  ^ 

D.  Preparation  of  TeF^  33 

E.  Mass  Spectrometric  Analysis  80 

RESULTS  AND  DISCUSSION  85 

A.  Chemical  Reduction  of  Potassium  Tellurite  85 

B.  Microbiological  Reduction  Study  88 

C.  Natural  Isotopic  Abundance  Studies  9] 

BIBLIOGRAPHY  93 

APPENDIX  I 

CALCULATIONS  OF  ISOTOPIC  VIBRATIONAL  FREQUENCY 

SHIFTS  FOR  TELLURIUM  -  CONTAINING  MOLECULES  A1 

A.  Diatomic  Molecules  A1 

B.  Polyatomic  Molecules  A5 

TeFg  A7 


viii . 


TeBr^  Al4 

SCTe  Al8 

H2Te  A22 

TeO^  A30 

3 

APPENDIX  II 

IBM  7040  COMPUTER  PROGRAMS  FOR  CALCULATION  OF 
ISOTOPIC  VIBRATIONAL  FREQUENCIES  AND 
PARTITION  FUNCTION  RATIOS  A36 

A.  Calculation  of  Isotopic  Vibrational 

Frequencies  and  Partition  Function 

Ratios  for  Diatomic  Species  A37 

B.  Computation  of  Isotopic  Vibrational 

Frequencies  for  XY^  Molecules  A40 

C.  Computation  of  Isotopic  Vibrational 

Frequencies  for  XY^  Molecules  A43 

D.  Computation  of  Vibrational  Frequencies 

for  XYZ  Molecules  A46 

E.  Computation  of  Force  Constants  for  Nbn- 

linear  XY^  Molecules  A49 

F.  Computation  of  Isotopic  Vibrational 

Frequencies  and'  Partition  Function 

Ratios  for  Non-linear  XY^  Molecules  A55 

G.  Computation  of  Force  Constants,,  Apex 

Angle,  Isotopic  Vibrational  Frequencies, 
and  Partition  Function  Ratios  for 
Pyramidal  XY^  Molecules  A58 

H.  Computation  of  Partition  Function  Ratios 

for  General  Molecules 


A63 


ix. 


LIST  OF  TABLES 


I.  Occurrence  of  Tellurium 


Page 

31 


II.  Absolute  Abundances  of  the  Isotopes  of 

Tellurium  38 


III. 

IV. 


Q£/Qj[  for  Te-containing  Molecules 

Equilibrium  Constants  for  Te2^2^— ■ 
Exchange  Reactions 


6l 


62 


V.  Minimum  Values  of  ^139/^122  65 

VI.  Summary  of  the  Calculated  Ratios  of  Rate 

Constants  for  the  Reduction  of  K^TeO^  67 

VII.  Ion  Species  Resulting  from  the  Ionization 

of  TeFg  84 

VIII.  Te13°/Te122  Fractionation  Study  86 

IX.  Calculation  of  Vibrational  Frequency  Shifts 

for  Te"^2  and  Te222-containing  Molecules  A 6 


I 


LIST  OF  FIGURES 


Page 


Figure 

1. 

0  Growth  Rings  in  the  Jurassic  Belemnite 

12 

Figure 

2. 

34 

The  Distribution  of  SJ  in  Nature 

20 

Figure 

3. 

The  Biological  Sulphur  Cycle 

22 

Figure 

4. 

Variations  in  the  Se^/Se^  Ratio  Found 

For  Natural  Samples 

27 

Figure 

5. 

Fluorine  Line 

78 

Figure 

6. 

Typical  Tellurium  Spectra 

82 

Figure 

7. 

Te122/^1^0  Fractionation  Study 

87 

Figure 

8. 

Variation  of  Isotopic  Composition  of 

Product  with  Percent  Reaction 

89 

Figure 

9. 

Normal  Vibrations  of  an  Octahedral 

XYg  Molecule 

A  8 

Figure 

10. 

Normal  Vibrations  of  a  Tetrahedral 

XY^  Molecule 

A 15 

Figure 

11. 

Relationship  Between  Force  Constants 
for  H^Te 

A2  6 

Figure 

12. 

Normal  Vibrations  of  a  Pyramidal 

XY^  Molecule 

A32 

INTRODUCTION 


The  concept  of  Isotopes,,  new  and  revolutionary  in 
1910*  is  today  basic  to  our  understanding  of  the  struc¬ 
ture  of  matter. 

The  first  direct  evidence  of  the  existence  of  iso¬ 
topes  came  from  studies  of  radioactive  elements  and  their 
decay  products.  Soon  the  chemists  investigating  the 
decay  products  found  pairs  of  radioactive  elements  which 
were  inseparable  and  appeared  to  be  identical  in  their 
chemical  properties.  Soddy  (100)  in  1910  suggested  that 
such  substances  differing  in  their  physical  properties 
but  having  the  same  chemical  properties  be  called  !! iso¬ 
topes11.  He  suggested  that  the  existence  of  isotopes 
was  a  general  phenomenon  and  that  we  could  also  expect 
the  stable  elements  to  be  mixtures  of  these  inseparable 
species.,  or  isotopes. 

J.  J.  Thomson  in  1912  (128)  using  his  positive  ray 
apparatus*  showed  that  neon  wras  indeed  a  mixture  of  at 
least  two  kinds  of  atoms  having  atomic  weights  of  20  and 
22  since  the  positive  ion  beam  formed  from  neon  was  re¬ 
solved  into  two  beams.  This  was  the  first  evidence  for 
the  existence  of  isotopes  among  the  stable  elements. 
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Aston  (2)  and  Dempster  (25)*  subsequently  examined  the 
isotopes  of  many  other  elements.  During  the  two  decades 
to  follow,  practically  all  the  elements  had  been  inves¬ 
tigated.  Today,  the  known  stable  elements  consist  of 
mixtures  of  from  one  to  ten  isotopes  and  these  total 
to  about  300  different  stable  isotopes  in  addition  to 
the  radioactive  species.  The  radioactive  disintegration 
studies  suggested  that  isotopic  abundances  of  an  element 
in  a  deposit  could  be  altered  by  radioactive  decay. 
Richards  (90)  found  that  lead  associated  with  thorium 

ores  had  a  higher  atomic  weight  than  that  associated 

232 

with  uranium  ores.  This  was  reasonable  since  Th  decay 

208  238 

terminates  in  Pd  ' ;  whereas  U  decays  finally  to 
Pb206.  Ey  knowing  the  amount  of  radioactively  produced 
lead  present  relative  to  the  amount  of  undecayed  parent. 


the  age  of  lead  deposits  can  be  determined. 

Similarly,  Aldrich  and  Nier  (l)  observed  that  the 

4o  4o 

amount  of  Ar  associated  with  K  was  up  to  three  times 

40 

the  concentration  of  Ar  in  atmospheric  argon.  The 

40  4o 

amount  of  Ar1  produced  by  the  radioactive  decay  of  K 


has  been  used  to  date  the  time  of  formation  of  potassium 
minerals.  Many  other  radioactive  isotopes  have  been 
used  for  dating  purposes. 
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However,,  despite  all  the  information  about  isotopes 
that  had  been  discovered  prior  to  1930,  it  was  believed 
that  isotopes  were  identical  in  their  chemical  proper¬ 
ties  and  that  the  stable  isotopes  of  the  elements  were 
of  constant  abundance  in  nature.  But  soon  after  the  dis¬ 
covery  of  deuterium,  the  heavy  isotope  of  hydrogen,  by 
Urey,  Brickwedde,  and  Murphy  (112)  in  1932,  it  was 
shown  by  theory  and  experiment  that  the  isotopes  of  the 
light  elements  did  differ  in  their  chemical  properties 
and  that  the  separation  of  these  isotopes  by  chemical 
means  was  possible. 


A.  Equilibrium  Isotope  Effects 

In  1933  Urey  and  Rittenherg  (114)  calculated  the 
equilibrium  constants  for  the  reaction 


Hg(g)  +  r>2(g) 


^  2KD( g) 


at  different  temperatures.  They  found  these  to  differ 
from  unity  which  indicated  that  the  isotopes  of  hydrogen 
were  different  chemically.  The  reaction  was  carried  out- 
in  the  laboratory  at  various  temperatures  by  Rittenberg, 
Bleakney,  and  Urey  (91)  and  the  equilibrium  constants 
determined.  These  were  found  to  be  close  to  those  cal¬ 
culated  theoretically.  This  wTas  the  beginning  of  many 
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isotope  fractionation  investigations  which  compared 
experimental  results  with  theoretical  calculations. 

Urey  and  Greiff  (ll6)  calculated  the  equilibrium 
constants  for  some  of  the  light  elements  such  as  lith¬ 
ium*  boron*  carbon*  nitrogen*  chlorine*  and  bromine. 
Their  results  indicated  equilibrium  fractionations 
varying  from  1 $  to  10$  and  experimental  confirmations 
of  some  of  these  were  soon  obtained.  For  example*  Weber 
et  al  (127)  found  CO^  in  equilibrium  with  water  to  be 


,18 


enriched  in  0'LU  by  4.6$  at  0°C  as  compared  with  4.4$ 
predicted  by  theory.  Other  equilibrium  isotope  effects 
and  data  are  summarized  by  Urey  (118) .  Many  isotope 
exchange  reactions  involving  the  isotopes  of  the  light 
elements  have  been  studied  in  recent  years  and  the 
agreement  betwTeen  theory  and  experiment  has  been  excel¬ 
lent.  To  date*  (with  the  exception  of  this  thesis) 
selenium  is  the  heaviest  element  for  which  exchange  re¬ 
actions  have  been  investigated  (6o). 


E.  Kinetic  Isotope  Effects 

■  —  '  ■  ■■  ■  rn  mm  ■  I  r  ■  H»IM  M  . .  ■  — "  ■  ■  •  — *,r^ 

In  1932  when  Urey  and  Washburn  (113)  were  doing 
their  experiments  with  deuterium  and  protium*  they  dis¬ 
covered  that  partially  electrolyzed  water  was  enriched 
in  deuterium*  protium  having  been  evolved  faster  than 
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deuterium  at  the  cathode.  Two  years  later,  Bach,  Bon- 
hoeffer,  and  Fajans  (5)  showed  that  reacted  over 
three  times  faster  than  with  bromine.  The  same  year 
Farkas  and  Farkas  (37)  reported  similar  results  in  the 
photochemical  reactions  of  and  with  chlorine. 

Urey  and  Teal  (117)  and  Eidinoff  (31)  have  reviewed 
the  extensive  work  done  on  the  relative  reaction  rates 
of  and  -  containing  compounds. 

Kinetic  isotope  effects  were  similarly  found  for 
other  lighter  elements.  In  1949  Lindsay,  McElheran, 
and  Thode  (65)  and  Eigeleisen  and  Friedman  (12)  re¬ 
ported  carbon  isotope  effects  in  the  decomposition  of 
oxalic  acid  and  the  decarboxylation  of  malonic  acid. 
These  results  were  reported  on  further  by  Lindsay, 
Bourns,  and  Thode  (63*64)  in  1951  and  1952  and  were  in 
excellent  agreement  with  theory. 

To  date,  much  literature  is  available  on  the  math¬ 
ematical  calculation  of  isotopic  rate  constants.  In 
1935  Eyring  (36)  set  forth  the  theory  of  absolute  rates 
from  which  the  ratio  of  rate  constants  for  reactions  of 
isotopic  molecules  can  be  calculated,  and  this  has  been 
revised  by  Bigeleisen  (11)  to  an  expression  which  is 
more  convenient  to  use.  Several  authors  (23*  92*  106) 
have  written  reviews  on  kinetic  isotope  effects. 
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C.  Mechanisms  of  Isotope  Fractionation  in  Nature 

Natural  isotope  fractionation  mechanisms  can  be  sub¬ 
divided  into  physical,,  chemical*  and  biological  processes. 
1.  Physical  processes  of  natural  isotope  fraction¬ 
ation 

a.  Diffusion :  The  mass  dependence  of  diffusion 

in  the  solid*  liquid*  or  gaseous  state  can  lead  to  iso¬ 
topic  fractionation  and  this  is  one  of  the  principles 
used  in  commercial  isotope  separation.  Studies  of  nat¬ 
ural  diffusion  phenomena  for  many  isotopes  have  been 
reported.  An  interesting  example  is  the  studies  of  ni¬ 
trogen  isotopes  by  Hoering  (53) •  He  discovered  that  ni- 

14 

trogen  contaminants  in  natural  gas  were  enriched  in  N 
over  the  nitrogen  found  in  the  associated  crude  oil. 

He  suggested  that  migration  of  the  oil  and  gas  through 
porous  deposits  was  responsible.  Migration  by  molecular 
flow  in  which  the  mean  free  path  of  the  gas  is  greater 
than  the  pore  dimensions  results  in  flow  rates  of  the 
isotopic  molecules  which  are  inversely  proportional  to 
the  square  root  of  the  masses.  This  postulate  was  later 
verified  by  a  laboratory  diffusion  experiment  and  he 
was  able  to  produce  to  some  degree  the  fractionation 
found  in  nature. 
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With  lighter  isotopes  and  low  temperature  condi¬ 
tions,  the  isotopic  alteration  can  he  quite  high.  This 
is  exemplified  h y  a  recent  lithium  diffusion  experiment 
in  titanium  dioxide  by  Johnson  and  Krouse  (55) . 

Lithium  was  allowed  to  diffuse  to  a  depth  of  1.7  mm  at 

room  temperature  and  the  innermost  lithium  was  found 

6 

to  be  7  percent  enriched  in  Li  in  comparison  to  the 
undiffused  lithium  showing  that  the  Li^  diffuses  quite 

7 

noticeably  faster  than  the  Li  . 


b.  Isotopic  vapour  pressures:  Friedman  (39) 

suggested  that  the  different  vapour  pressures  of  the 
isotopic  compounds  might  influence  the  natural  abundance 
of  each  isotope.  He  showed  that  there  was  a  gradual 
depletion  of  the  deuterium  content  in  rainfall  samples 


progressing  from  the  Pacific  Coast  of  North  America  in¬ 
land  over  the  Rocky  Mountains.  He  found  further  that 
ocean  waters  near  the  equator  were  found  to  be  enriched 
in  deuterium  in  comparison  to  ocean  waters  polewards. 
Both  these  phenomena  are  c  au  s e  d  by  the  preferential 


evaporation  of  Ho0  over  HDO.  In  the  rainfall  study,  the 
EDO  condenses  faster  enriching  the  coastal  regions  in 
HDO  while  the  remaining  vapour  enriched  in  H^O  moves 


further  inland . 
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The  oxygen  isotopes  of  the  earth's  water  bodies 

1  ft 

behave  similarly.  Epstein  (34)  reported  that  the  C> 
content  of  the  ice  is  very  low  in  the  Arctic  regions 
relative  to  ocean  waters  and  that  the  actual  value  in 
the  ice  cap  varies  in  a  cyclic  fashion  with  depth.  This 
cyclic  phenomenon  which  is  due  to  seasonal  changes  may 
be  explained  by  thinking  of  the  atmosphere  between  the 
equator  and  the  poles  as  a  distillation  column.  Water 
is  evaporated  from  the  oceans  at  the  equator,,  rises  up¬ 
ward  at  the  equator,  and  spreads  out  toward  the  poles. 

In  the  distillation  of  the  water,  the  simple  process 
separation  factor  a  is  given  by  the  ratio  of  vapour  pres¬ 
sures  of  the  two  isotopic  water  species  as 

PHp0l6  1.008  at  25°C 

a  =  — =■— =-tt  = 

PH2010  1.011  at  0°C 

"1  ft  1 8 

Since  H.^O1  has  the  higher  vapour  pressure,  the  0  is 

favoured  in  the  liquid  and  will  be  depleted  in  the  vapour 

As  rain  falls  in  the  different  latitudes,  the  vapour 

spreading  out  towards  the  poles  becomes  more  and  more 

depleted  in  the  heavy  isotope.  We  therefore  have  a  multi 

stage  process.  The  extent  of  the  isotope  fractionation 

which  in  this  case  is  the  depletion  of  0^  in  snow  fall- 
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ing  on  the  Greenland  Ice  caps  will  depend  on  the  number  of 

stages  in  the  process  and  on  the  difference  in  temperature 

between  the  equator  and  the  poles.  Thus.,  the  depletion 
1 8 

of  C>  in  the  snow  and  ice  in  the  Arctic  will,,  in  general., 

be  higher  in  winter  than  summer.  Seasonal  fluctuations 
l8 

in  the  (j  content  with  depths  have  shown  up  in  the  layers 
of  ice  in  the  Greenland  glacier.  From  such  studies,  es¬ 
timates  of  the  age  of  glaciers  can  be  determined  by  simply 
1 8 

counting  Cl  concentration  rings.  Such  fractionations 
caused  by  differences  in  vapour  pressure  are  termed 
"Rayleigh"  fractionations. 

c.  Miscellaneous ;  It  is  expected  that  dif¬ 
ferences  in  the  isotopic  solubilities,  melting  and  sub¬ 
limation  points  might  be  other  mechanisms  of  natural  iso¬ 
tope  fractionation.  Very  little  work  to  date  has  been 
reported  on  these  physical  processes, 

2.  Chemical  processes 

a.  Equilibrium  isotopic  exchange :  Friedman 

suggested  further  that  equilibrium  isotope  exchange  re¬ 
actions  would  produce  fractionation  of  isotopes  in  nat¬ 
ural  processes.  For  example,  a  study  of  fumaroles  in 
Yellowstone  Park  showed  the  deuterium  content  to  be  de¬ 
pleted  by  as  much  as  45$.  He  suggested  that  this  might 
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be  the  result  of  the  equilibrium  exchange  reaction 

K 


H2  +  KDO 


-  HD  +  H20 


The  equilibrium  constant  for  exchange  reactions  is  de¬ 
pendent  upon  temperature,  and  if  this  hypothesis  is  the 
only  contributing  factor,  theoretical  considerations, 
based  on  the  observed  deuterium  depletions,  place  the 
temperatures  of  these  fumaroles  above  400°C. 

In  the  exchange  of  the  oxygen  isotopes  between 
carbonate  and  water,  the  heavy  isotope  O'1  is  favored 
in  the  carbonate  over  that  of  that  water.  This  reaction 


can  be  written  in  the  form 


16* 


1/3  co3  +  h2 


0 


18 


K 


.18= 


-  1/3  C03  +  H20 


16 


The  equilibrium  constant  has  been  determined  by  Clayton 
(20)  to  be  1.0301  at  25°C.  meaning  that  at  25°C.  lime¬ 
stone  and  calcium  carbonate  material  will  be  enriched 
in  O'1  over  that  of  sea  water  by  about  3/.  The  exact 
extent  of  this  enrichment  will  depend  on  the  temperature 
of  the  water  at  the  time  of  deposition.  The  reaction 
has  been  studied  in  detail  by  Urey  (119)  and  he  suggested 

that  it  might  be  used  in  the  determination  of  temperatures 

l8 

and  that  the  accurate  measurement  of  the  O'1  content  of 
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suitable  marine  carbonate  sediments  would  indicate  the 
temperature  of  the  seas  at  the  time  of  deposition. 

According  to  Clayton  the  temperature  T  is  given 
by  the  equation: 

In  K  *  2725  T"2 

As  a  direct  result  of  this  temperature  dependence,  the 

l8 

concept  of  the  0  or  paleo-temperature  scale  was  for¬ 
mulated. 

Studies  were  carried  out  by  Urey  et  al  (120), 

Epstein  et  al  (32),  and  McCrea  (73)  using  well  preserved 
samples  of  belemnites  from  the  U,  S.,  Sweden,  Denmark, 
England,  France,  Holland,  and  Germany  and  these  showed 
that  the  temperatures  which  were  determined  seemed  rea¬ 
sonable  . 

A  fascinating  study  was  made  by  Epstein  (33)  of  the 

j  8 

0~  content  in  the  carbonate  shells  generated  by  marine 
animals.  The  subject  of  investigation  was  a  Jurassic 
belemnite,  an  extinct  animal  which  flourished  about  120 
million  years  ago.  Successive  growth  rings  of  the  shells 
showed  a  cyclic  variation  in  the  heavy  oxygen  (O'1  )  con¬ 
tent  which  were  associated  with  seasonal  variations  in 
temperatures.  Figure  1  shows  these  seasonal  vibrations 

“i  O 

in  temperature  as  indicated  by  the  O'1  content  of  the 
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growth  rings  of  the  belemnite.  Such  an  oxygen  Isotope 
study,  therefore,  not  only  gives  us  information  about 
climatic  conditions  in  various  parts  of  the  world  mil¬ 
lions  of  years  ago,  but  tells  us  something  about  the 
lives  and  habits  of  animals  that  now  are  extinct. 

b .  Kinetic  isotope  effects:  Due  to  the  dif¬ 

ficulty  in  observing  purely  chemical  kinetic  isotope 
effects  in  nature,  very  little  has  been  reported  in  the 
literature.  Numerous  reports  have  been  given  for  kinetic 
isotope  fractionation  in  biochemical  processes  and  will 
be  discussed  later  in  connection  with  sulphur  isotope 
fractionation. 

3.  Biological  processes 

Biological  and  microbiological  isotope  frac¬ 
tionations  are  perhaps  due  to  both  physical  and  chemical 
processes.  Because  of  the  complexity  of  the  mechanisms, 
there  is  insufficient  data  to  verify  whether  isotopes 
are  fractionated  in  such  processes  as  enzyme  complexing 
and  diffusion  through  biological  membranes.  Significant 
isotope  fractionations  have  been  realized  with  carbon, 
nitrogen,  oxygen,  and  sulphur  in  biochemical  processes 
and  the  results  in  many  instances  are  consistent  with 
similar  inorganic  chemical  processes. 
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Nier  and  Gulbransen  (83)  and  two  years  later ,  Nler 
and  Murphy  (84)  found  that  In  general,  limestones  were 
enriched  and  plants  depleted  In  the  heavier  carbon  Iso¬ 
tope.  They  reported  variations  up  to  5$  in  the  C^/C1^ 
abundances.  Many  workers  have  since  re-affirmed  these 
results  and  can  be  found  in  the  literature  (6,22,38,62, 
71,80,86,97,107,110,130,133,134,135).  One  of  the  many 

interesting  results  of  these  investigations  was  that 

12 

both  land  and  marine  plants  tend  to  concentrate  C  , 
although  the  latter  do  this  to  a  lesser  degree.  This 
phenomenon  may  be  attributed  to  differences  in  the  cel¬ 
lular  structures  of  the  two  types  of  plants,  assuming 
that  the  concentration  of  the  C  is  due  to  kinetic  and/ 
or  equilibrium  isotope  effects  in  the  photosynthesis  of 
atmospheric  CO^. 

Much  investigation  of  the  isotope  fractionation 
due  to  bacteria  has  been  carried  out.  The  most  fruitful 
of  these  investigations  have  been  carried  out  with 
sulphur,  and  this  will  be  discussed  in  more  detail  later. 
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D.  The  Group  ^TA  Elements 

1.  General  chemistry  (47)*  (96) 

The  group  VI  A  elements  include  oxygen,  sulphur, 
selenium,  tellurium,  and  polonium  and  are  character:!  zed 
by  6  valence  electrons.  They  differ  in  nuclear  charge, 
in  size,  and  in  the  number  of  inner  electron  "shells”. 
Little  is  known  about  polonium  since  work  on  i.t  has 
been  retarded  because  of  its  high  radioactivity.  It 
occurs  naturally  in  uranium  minerals  such  as  pitchblende 

—  Q 

but  only  to  the  extent  of  5  x  10  per  cent  of  the 
mineral.  Thus  its  occurrence  is  so  small  that  its  . 
significance  in  stable  Isotope  fractionation  processes 
is  negligible. 


0  S  So 


The  Group  VIA  Elements 
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The  firmness  with  which  the  outer  electrons  are 
held  increases  with  increasing  nuclear  charge,  but  is 
diminished  as  the  distance  from  the  nucleus  and  the 
number  of  inner  electron  shells  increases.  Since  the 
two  latter  effects  predominate,  it  may  be  predicted 
that  the  electrons  are  most  tightly  held  by  the  oxygen 
atom  and  most  easily  lost  (or  shared)  by  the  tellurium 
atom.  This  assumption  accounts  for  and  is  in  accord 
with  the  systematic  variation  in  physical  and  chemical 
properties  which  these  elements  show. 

The  elements  of  the  group  readily  form  compounds 
or  ions  by  sharing  or  adding  2  electrons,  provided  by 
atoms  of  other  elements,  to  complete  the  octet  of  the 
inert-gas  structure  thereby  acquiring  an  oxidation  state 
of  -2.  This  tendency  is  most  prominent  in  the  case  of 
the  smallest  atom,  that  of  oxygen,  and  least  in  the  case 
of  the  largest  atom,  that  of  tellurium.  In  addition, 
the  atoms  of  sulphur,  selenium,  and  tellurium  may  lose 
(or  share)  their  valence  electrons  so  as  to  acquire  ox¬ 
idation  states  up  to  +6.  Oxygen  does  not  show  this  tend¬ 
ency,  because  its  electrons  are  so  tightly  bound  that  no 
other  element,  with  the  possible  exception  of  fluorine, 
is  able  to  remove  them.  The  change  in  oxidation  state 
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of  an  element  during  a  chemical  reaction  Is  one  of  the 
factors  which  determines  the  extent  to  which  Isotope 
abundances  may  be  altered  In  the  laboratory  or  by  nat¬ 
ural  processes. 

All  of  the  elements  of  the  group  show  allotropy. 

Just  as  oxygen  can  exist  as  diatomic  0^  and  triatomic 
0^,  the  other  elements  can  be  obtained  In  more  than  one 
form.,  the  forms  differing  either  In  the  number  of  atoms 
per  molecule.,  or  in  the  arrangement  of  molecules  in  the 
solid.  There  is.,  going  from  sulphur  to  tellurium,  an 
increasing  disposition  toward  formation  of  long  chains 
of  atoms  held  together  by  covalent  bonds.  In  general, 
the  larger  the  atom,  the  less  the  tendency  to  form  mul¬ 
tiple  bonds  and  the  greater  the  tendency  of  each  atom 
to  be  bound  to  more  than  one  atom.  Oxygen  stands  alone 
from  the  group  in  being  a  diatomic  gas  at  room  temperature. 

As  would  be  expected  because  of  the  increasing 
number  of  electronic  shells,  the  ionic  radii  increases 
from  1 . 40A  for  oxygen,  1.84a  for  sulphur,  I.98A  for  sele¬ 
nium,  to  2.4A  for  tellurium  (47).  Also  of  interest,  the 
oxidation  potential  for  each  of  these  elements  is  as 
follows:  oxygen,  -1.23  volts;  sulphur,  -0.l4  volts; 

selenium,  +0.40  volts;  and  tellurium,  +0.7'2  volts  (96). 
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Thus,  there  is  a  much  greater  tendency  for  oxygen  to 
form  H^O  than  for  tellurium  to  form  H^Te.  Unlike  H^O 
and  H^S,  H^Se  and  H^Te  are  better  reducing  agents  than 
hydrogen . 

From  the  foregoing  discussion,  it  is  readily  seen 
why  the  properties  of  sulphur,  selenium,  and  tellurium 
vary  from  the  decidedly  non-metallic  nature  of  sulphur 
to  the  mildly  metallic  character  of  tellurium.  For  ex¬ 
ample,  the  free  element  sulphur  is  a  non-conductor  of 
electricity,  selenium  is  a  poor  conductor  (but  its  con¬ 
ductivity  increases  proportionately  to  its  illumination), 
and  tellurium  is  somewhat  better  than  selenium. 

2 .  Isotope  fractionation  studies 

In  considering  isotope  fractionation  in  the 
chemical  group  VIA,  oxygen  has  received  extensive  study 
and  some  examples  have  been  previously  cited.  It  would 
seem,  however,  that  tellurium  isotopes  should  resemble 
sulphur  and  selenium  in  their  behavior  because  of  the 
similarities  of  their  valence  states.  Therefore,  some 
of  the  significant  findings  of  sulphur  and  selenium  iso¬ 
tope  fractionation  studies  are  reviewed  here. 

a.  Sulphur:  Extensive  studies  of  the  varia¬ 


tions  in  the  isotopic  abundances  of  sulphur  have  been 
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carried  out.  It  was  clear  that  differences  In  the  chemical 
properties  of  the  sulphur  Isotopes  did  exist  and  that  isotope 
fractionation  could  be  expected  in  natural  processes. 

For  example.,  calculations  of  the  isotope  exchange 
constant  for  the  reaction 

S320^  +  HgS34  * . ^  s340=  +  HgS32 

gives  K  at  25° C  to  be  1.071*  thus  favoring  the  enrich- 
ment  of  S  in  the  sulphate  by  about  seven  percent  (ill). 

It  was  therefore  suspected  that  a  trend  toward  this  fa¬ 
voured  distribution  existed  in  nature  with  sulphates  en¬ 
riched  and  sulphides  depleted  in  the  heavy  isotope  of 
34 

sulphur,,  . 

The  initial  data  published  by  Thode,  Macnamara,  and 

Collins  (104)  showed  that  the  S^/S^  ratio  in  terrestial 

samples  varied  by  as  much  as  5$.  As  predicted.,  sulphates 

34 

were  generally  found  to  be  enriched  in  S  while  sulphides 
were  depleted  in  the  heavier  isotope.  Figure  2  shows  the 
general  distribution  of  sulphur  isotopes  in  nature.  The 
sulphur  isotope  distribution  data  are  expressed  in  terms 
of  6S^  °/oo  defined  as  follows: 
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FIG.  2  THE  DISTRIBUTION  OF  S34  IN  NATURE 
Taken  from  Thode  (103) 
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( o34  /o32x  / q34 /o32n 

6  S34  0 /oo  =  '  ^sample  '  ^standard 

( S34/S32) 
v  /  >  standard 

Positive  and  negative  6 -values  mean  that  the  samples  are 
enriched  or  depleted  in  respectively  with  respect  to 
the  standard. 

It  was  reported  by  Macnamara  and  Thode  (68)  that  the 
sulphur  in  meteorites  had  a  remarkably  constant  S32/S34 
ratio  which  was  close  to  the  average  S32/S34  ratio  for 
t-errestial  samples.  These  facts  led  them  to  believe  that 
the  meteoritic  S  /SJ  ratio  represented  the  primordial 
abundance  of  terrestial  sulphur  before  any  fractionation 
occurred  or  the  ratio  of  the  sulphur  isotopes  in  the  earth 
at  the  time  that  it  was  formed.  It  is  clear  that  since 
the  earth  was  formed,,  the  sulphur  in  the  crust  has  been 
fractionated  in  biological  and  geological  processes  due  to 
differences  in  the  chemical  properties  of  the  isotopes. 

The  exchange  cited  does  not  take  place  under  normal 
conditions  (124)  and  all  efforts  to  establish  it  under 
laboratory  conditions  for  a  reasonable  length  of  time 
failed.  Therefore,  Tudge  and  Thode  (ill)  suggested  that 
the  following  biological  sulphur  cycle  might  well  provide 
the  mechanism  for  this  fractionation  in  nature  (Fig.  3) . 


NATIVE  SULPHUR 


OXIDATION  AND  AER03IC  BACTERIA 


FIG.  3  THE  BIOLOGICAL  SULPHUR  CYCLE  IN  NATURE 
Taken  from  Thode  (108) 
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According  to  Thode  (108)  the  bacterial  reduction  process 
is  found  to  be  a  unidirectional,  not  an  equilibrium  pro¬ 
cess,  and  isotope  fractionation  occurs  because  the  light- 
isotope  containing  compounds  react  faster  than  the  spe¬ 
cies  containing  the  heavier  isotope.  Thode  goes  on  to 
say  that  in  particular,  the  bacterial  reduction  of  sul¬ 
phate  is  largely  responsible  for  the  isotopic  variations 
which  occur  in  the  oceans,  evaporites,  and  materials  in 
the  sediments. 

Many  studies  have  been  made  of  isotope  fractionation 
in  the  sulphur  cycle.  Anaerobic  reducing  and  photosyn¬ 
thetic  hydrogen  sulphide  oxidizing  bacteria  abound  in  the 
Cyrenacian  lakes  of  Africa.  These  produce  a  sedimentary 
layer  of  elemental  sulphur  (15)  and  analysis  of  samples 
from  these  lakes  gave  direct  evidence  for  the  fraction¬ 
ation  of  sulphur  in  the  sulphur  cycle,  Macnamara  and 

oh 

Thode  (69)  found  the  content  of  the  elemental  sulphur 

to  be  3.2$  less  than  that  of  the  soluble  sulphate  from 
the  same  lake. 

Further  evidence  is  provided  in  the  salt  domes  of 
Louisiana  and  Texas  (107)  where  elemental  sulphur  was  de- 
pleted  by  about  and  sulphide  by  about  5^  in  rel¬ 

ative  to  the  associated  gypsum.  Two  possibilities 


. 
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exist  for  the  formation  of  this  sulphur  by  reduction  of 
the  sulphate:  -  either  the  high  temperature  reduction  of 
the  sulphate  by  organic  matter  (oil  fields  associated 
with  the  salt  domes)*  or  by  the  bacterial  reduction  of 
the  sulphate  with  petroleum  providing  a  source  of  carbon 
and  free  energy  for  life.  The  first  of  these  mechanisms,, 
however,,  seems  incapable  of  causing  the  large  fraction¬ 
ation  observed;  whereas*  the  latter  is  quite  capable  of 
doing  so*  as  shown  by  the  Cyrenacian  Lake  study.  Miller 
(78)  added  further  evidence  to  the  "bacterial  reduction" 
postulate  by  isolating  anaerobic  sulphate-reducing  bac¬ 
teria  in  these  deposits.  Thus  the  production  of  carbon 
dioxide  by  these  sulphate-reducing  bacteria  alters  the 
calcium  sulphate  to  a  porous  calcitic  limestone  with  el¬ 
emental  sulphur  embedded  in  it  and  this  accounts  for  the 
type  of  deposit  that  is  present  in  these  salt  domes  today. 

The  hypothesis  that  sulphate-reducing  bacteria  were 
responsible  for  the  sulphur  isotope  fractionation  has  been 
substantiated  by  laboratory  experiments.  Thode*  Kleer- 
ekoper*  and  McElheran  (105)  showed  that  Desulphovibrio 
desulfuricans  while  reducing  sulphate*  produced  HgS  which 
was  ifo  richer  in  S^2  than  the  nutrient.  In  subsequent 
studies*  Harrison  and  Thode  (44)  found  that  this  fraction- 
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ation  varied  with  reaction  conditions  and  that  the  max¬ 
imum  enrichment  was  2.5^  at  25°C.  Kaplan  and  Rittenberg 
(57)  in  1964  reported  enrichments  as  high  as  4.6$.  They 
showed  also  that  other  micro-organisms  such  as  _S.  cer- 
evisiae  were  capable  of  carrying  out  fractionation  of 
oxidized  sulphur  compounds  in  reduction  experiments  with 
equal  or  greater  efficiency. 

McElheran  in  1951  (74)  was  unable  to  find  any 
fractionation  in  the  oxidation  of  sulphur  to  sulphate 

by  Thiobacillus  but  in  1964  Kaplan  and  Rittenberg  (57) 

32 

reported  an  enrichment  of  1.8$  of  S  in  oxidizing  H^S 

to  SO^  by  chemosynthetic  oxidation  and  an  enrichment  of 

1.0$  in  photosynthetic  oxidation.  Nakai  and  Jensen  (96) 

published  results  in  1964  for  the  bacterial  reduction 

and  oxidation  of  sulphur  using  marine  mud  cultures.  Values 

on  oh 

for  the  ratio  of  the  isotopic  (SD  and  )  rate  constants 
of  about  1.020  were  found  for  the  reduction  experiments, 
and  values  ranged  from  1.0003  to  1.0017  for  the  oxidation 
experiment.  Fractionation  factors  between  sulphate  and 
sulphide  in  bottom  sediments  vary  between  1.03  and  1.04 
for  sea  samples  and  1.005  to  1.009  for  fresh  water  samples, 
indicating  that  the  sulphate  produced  becomes  more  en¬ 
riched  in  S32  than  does  sulphur  or  sulphide. 
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b.  Selenium:  An  investigation  into  the  iso¬ 

topic  fractionation  of  selenium  was  conducted  by  Krouse 
and  Thode  (60).  They  reported  Se^/Se^  ratios  as  found 
in  figure  4.  The  Se^/Se^  ratio  varied  1.5$  in  the 
samples  examined.  The  distribution  bears  some  similar¬ 
ity  to  that  of  the  sulphur  isotopes.  Samples  of  plant 
materials  and  soil  showed  the  largest  variation  in  iso¬ 
tope  ratios;  whereas,  selenides  from  massive  sulphide 
ores  showed  little  or  no  variation  with  respect  to  mete¬ 
or  itic  selenium. 

They  further  reported  that  Se'  0^  was  reduced  1.5$ 

82  = 

faster  than  Se  0^  to  elemental  selenium.  This  result 

is  in  agreement  with  a  simple  theoretical  calculation. 

This  compares  to  a  2.3$  kinetic  isotope  effect  in  the 

S2  =  34  = 

chemical  reduction  of  S  0^  and  S  '0^  reported  by  Harrison 
and  Thode  (44) . 

Further  investigation  into  the  isotope  fractionation 
of  selenium  has  been  carried  out  by  Rees  (89) .  Among  other 
experiments,  he  carried  out  investigations  into  the  kin¬ 
etic  isotope  effects  in  the  reduction  of  SeIV  to  Se°,  and 

Ty  yi 

into  the  isotope  exchange  reaction  between  Se  and  Se  . 

He  reported  a  kinetic  isotope  effect  varying  from  1.6$  at 
4o°C  to  2.0$  at  24°C  for  a  20$  reduction  of  sodium  selenite. 
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and  phase  fractionations  of  2.7$  and  4.2$  respectively 

for  59  -  and  146  -  day  equilibrations  of  SeVI  and  SeIV. 

In  each  case,  the  higher  valence  states  were  enriched 

82 

in  the  heavier  Se 

Because  of  the  numerous  reports  of  isotopic  fraction¬ 
ation  of  sulphur  and  selenium,  a  similar  fractionation 
of  tellurium  isotopes  by  natural  processes  is  strongly 
predicted . 

3 .  Tellurium 

a.  History  and  occurrence  -  The  discovery  of 
tellurium  was  made  late  in  the  eighteenth  century  (78) . 

A  peculiar  mineral  known  as  white  gold-ore  or  grey  gold- 
ore  occurring  in  the  sandstone  near  Zalathna,  Transyl¬ 
vania,  was  considered  by  chemists  and  mineralogists  as 
a  kind  of  alloy  of  antimony  and  bismuth.  Various  forms 
were  given  a  host  of  names  but  not  until  1782  that  F.J. 

n 

Muller  von  Reichenstein  showed  that  there  was  neither  any 
bismuth  nor  antimony  in  the  new  mineral,  and  suspected 
that  a  new  element  had  been  discoverd.  Then  sometime 
after  1789,  M.  H.  Klaproth  extracted  the  unknown  metal 
and  observed  some  of  its  properties  from  which  concluded 
that  the  mineral  contained  a  peculiar  and  distinct  metal 
essentially  different  from  every  other  metallic  substance 
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hitherto  known.  It  was  he  who  named  it  tellurium  after 
the  earth  (Latin  -  tellus) . 

Tellurium  occurs  in  only  a  few  places  and  in  small 
quantities;  it  is  less  abundantly  distributed  than  any 
other  element  in  the  sulphur  family.  According  to  the 
data  of  spectrographic  analysis,,  it  is  absent  in  the 
spectra  of  the  sun  and  the  stars,  however,  it  could  be 
explained  by  the  low  sensitivity  of  the  spectrographic 
method  with  respect  to  these  elements  (98) .  According 
to  F.  W.  Clarke  and  H.  S.  Washington  (19),  the  igneous 

rocks  on  the  earth,'1  s  crust  contain  8  x  10  percent 

-8  -Q 

sulphur,  n  x  10  percent  selenium,  and  n  x  10  percent 

tellurium;  or  tellurium  is  about  as  abundant  as  gold. 

It  occurs  combined  with  gold,  silver,  bismuth,  and  many 

other  metals.  The  chief  tellurium  minerals  are  the  tel- 

lurides,  but  some  tellurites  and  tellurates  exist  in 

natural  form.  These  are  listed  in  Table  I.  Tellurium 

has  not  been  detected  in  soils  (123).  It  is  found  in 

many  volcanic  areas  and  has  been  found  in  "minute"  amounts 

in  some  meteorites  (29) .  It  is  obtained  principally  as 

a  by-product  of  the  refining  of  other  minerals  with  which 


it  combines. 
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b.  Uses  -  To  cite  some  uses,  it  is  used  in 
coloring  glass,  preparation  of  organic  dye-stuffs,  man¬ 
ufacture  of  high  resistance  alloys,  staining  silver,  as 
a  delicate  test  of  sterilization  in  bacteriology,  ident¬ 
ification  of  presence  of  diphtheria  and  other  germs,  in 
motor  fuels  as  an  anti-knock  constituent,  insecticides, 
germicides,  fungicides  and  wood  preservatives,  and  in  the 
manufacture  of  extremely  hard  alloys  having  very  great 
tensile  strength.  Its  main  use  in  recent  years  is  in  the 
semiconductor  field  where  it  has  found  use  as  an  intrin¬ 
sic  photoconductive  detector  (77)  and  also  in  rectifiers 
(96). 


c.  Physical  properties  -  Tellurium  exists  pri¬ 
marily  as  "metallic"  tellurium  which  is  the  most  stable 
allotrope.  This  form  is  grey  black  in  color  with  a  metal¬ 
lic  luster  and  of  hexagonal  crystal  structure.  In  1913 
Cohen  and  Kroner  (21)  claimed  that  tellurium  existed  in 
two  dynamic  allotropes,  and  according  to  Mellor  (76) 
tellurium  is  obtained  in  the  amorphous  form  as  a  black 
powder  when  tellurium  is  precipitated  from  an  acid  solution. 
The  amorphous  form  suffers  no  change  upon  heating  until  it 
is  melted  (m.p.  452°C)  and  then  it  is  transformed  into  the 
very  brittle  and  easily  powdered  white-grey  crystalline  form. 
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Table  I  OCCURRENCE  OF  TELLURIUM  (7 6),  (98) 


Native  Element 
Tellurides 

Altalte 

Arseniotellurlte 

Arsenotellurlte 

Calaverlte 

Coloradoite 

Cslklovalte 

Frohberglte 

Goldf ieldlte 

Goldschmidt It e 

GrGnlingite 

Hedleyite 

Hesslte 

Hydrogen  tellurlde 

Jos8lte 

Krennerlte 

Melonlte 

Montbrayite 

Muthmannite 


Te 


PbTe 

Te^As^O^ 

Te^As^S^ 

AuTe^ 

HgTe 

Bl2TeS2 

FeTe2 

Cu^Sbg ( Sj Te) ^ 
Au2AgTeg 
Bl^S^Te 
Bl^Te^ 

Ag2Te 
H2Te 
Bl^TeS 
( Au,  Ag)  Te 
NigTe^ 

AUgTe^ 

( Ag ,  Au )  Te 
Au  2  Sb  2  P  b^Te  g  0  ^  (~ 


Nagyaglte 
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Oruetite 

Bi^Te^Bi^S^ 

Petzite 

(AgJ,Au)2Te 

Pilasonite 

Bi^Te2 

Rickardite 

Cu^Te^ 

it 

Stutzite 

Ag^Te 

Sylvanite 

AuAgTe^ 

Talapite 

Bi(SJTe)3Ag3 

Telluro  bismuthite 

Bi2Te^ 

Tetradymite 

Bi2Te2 

Vondiestite 

( Ag, Au) ^BiTe^ 

Wehrlite 

Bi^Te^Ag 

Weissite 

Cu2Te 

Mixtures 

Coolgardite 

coloradoite,  pet¬ 
zite,,  calaverite, 
sylvanite 

Henryite 

lead  telluride  and 
pyrites 

Kalgoorlite 

coloradoite,  pet¬ 
zite,  and  tellurium 

White  tellurium 
(Gelberg,  mullerite) 

antimonial  telluride 
of  gold,  silver, 
and  lead 

Tellurites 


Durdenite 

Tellurite 

Fe2(Te03)3.4H20 

TeOg 

Tellurite 
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Dunhamite 
Emmons it e 
Mackayite 


Pb TeO^ 

Fe2 (TeOQ) ^ .2Ho0 
Fe2(TeO  )  nH20 


Tellurates 

Ferrotellurite 

Magnolite 

Montanite 


FeTeO^ 

HggTeO^ 

Bi(0H)2Te04 
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Textbooks  as  recently  as  1953  claimed  a  hexagonal  lat¬ 
tice  and  the  existence  of  two  polymorphic  modifications 
a  and  P  with  a  transition  point  at  354°C.  However,  Bose 
(l4)  in  1941  showed  by  x-ray  diffraction  that  precipi¬ 
tated  tellurium  is  not  amorphous  but  is  hexagonal  in  a 
finely  divided  state,  ie .  it  is  the  same  as  metallic 
tellurium,  and  Dodwell  (27)  proved  in  1945  that  tellurium 
does  not  change  its  structure  when  heated  up  to  360°C 
and  consequently  does  not  have  any  polymorphic  varieties. 
Crystallization  of  tellurium  fused  in  sodium  hydrosul¬ 
phate  and  cooled  in  carbon  dioxide  results  in  tine  formation 
of  a  crystalline  star  resembling  that  of  antimony. 

d.  Chemical  properties  -  In  agreement  with  the 
general  observation  that  the  basic  properties  of  the 
natural  families  of  elements  increases  with  atomic  weight, 
tellurium  exhibits  a  greater  basicity  than  selenium  and 
sulphur . 

H.  Davy  (24)  was  the  first  to  discover  hydrogen 
telluride  ( H^  Te ) ,  which  he  formed  by  the  electrolysis  of 
water  plus  acid  using  a  tellurium  cathode.  H^Te  is  also 
formed  by  the  action  of  acids  on  tellurides,  but  the  elec¬ 
trolytic  process  gives  the  best  yields.  It  is  a  colorless 
gas  with  a  foul  smell  resembling  that  of  hydrogen  sulphide. 
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Heating  of  tellurium  in  the  presence  of  oxygen 
produces  tellurium  dioxide,  TeO^ .  It  forms  small,  color¬ 
less,  octahedral  crystals.  The  hydrated  form  of  telluri¬ 
um  dioxide  is  tellurous  acid,  H^TeO^.  Tellurium  monoxide 
is  formed  when  tellurium  sulphite  is  heated  between  l80° 

and  230°C  in  vacuo  (TeSO^  - >  TeO  +  SO^)  .  It  is  a 

brown,  porous  mass,  stable  in  dry  air  but  decomposes  when 
heated  into  tellurium  and  tellurium  dioxide.  Tellurium 
trioxide  (TeO^)  is  formed  when  telluric  acid  (H^TeO^)  is 
heated  above  360°C  but  below  red  heat.  If  the  tempera¬ 
ture  be  too  high,  the  trioxide  decomposes  into  the  diox¬ 
ide.  TeO^  ^-s  orange -yellow  and  crystalline  in  form. 

Tellurium  combines  with  the  halogens  to  produce  tel¬ 
lurium  halides.  Of  particular  interest  in  this  work  is 
the  relatively  unreactive  tellurium  hexafluoride  (TeF^) 
gas  and  its  possible  use  in  mass  spectrometric  determina¬ 
tion  of  the  tellurium  isotopes.  Both  hexa-halides  and 
tetra-halides  are  possible.  In  the  case  of  fluorine, 

TeFg  is  formed  at  temperatures  reaching  200°C  in  excess 
fluorine,  whereas,  TeF^  is  formed  at  temperatures  below 
l80°C  in  excess  tellurium  (5 6) . 

e .  Absolute  abundances  of  the  isotopes  of 

tellurium  -  The  first  report  of  tellurium 

126 

isotopes  came  in  1925  when  Aston  (3)  discovered  Te  , 
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128  180 

Te  ,  and  Te  using  his  mass  spectrograph.  As  to 
their  relative  abundance,  he  only  reported  that  the 
latter  two  were  of  equal  intensity  and  double  that  of 
the  first.  Then  in  1931  he  reported  Te12^  (4).  Bain- 
bridge  (7)  in  1932  reported  isotopes  122,  123,  124,  and 
127  and  Bartlett  (8)  in  1934  predicted  from  Am/Az  plots 
that  there  should  be  isotopes  Te21^  and  Te22<2.  Dempster 
(26)  in  1936  confirmed  Bainbridge's  and  Bartlett's  re¬ 
sults  using  his  mass  spectrometer  but  that  there  was  no 
Te"*"2^  nor  Te22^.  It  was  not  until  1942  when  R.  Migeotte 
(77)  using  a  Bainbridge  mass  spectrograph  reported  the 
eight  isotopes  with  percentage  abundances.  Four  years 
later,  Williams  and  Yuster  (132)  confirmed  these  eight 
isotopes  giving  relative  abundances.  White  and  Cameron 
two  years  later  reported  similar  results  (131).  The 
latest  values  of  the  isotopes  and  their  percentage  abun¬ 
dances  were  given  by  Hollander  et  al  in  1953  (54).  The 
results  of  these  investigations  are  summarized  in  Table  II. 

As  regarding  the  atomic  mass  of  each  of  the  isotopes, 
Duckworth  et  al  (28)  in  1951  first  reported  the  masses 
Te126,  Te128,  and  Te130  as  125.9427,  127.9471,  and  129.9467 
respectively  and  one  year  later,  Halsted  (43)  gave  the 
atomic  masses  for  all  eight  isotopes .  Up  to  i960,  these 
were  the  accepted  values  and  no  new  values  were  found  in 
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the  literature.  These  values  are  as  given  below: 


Mass  Number 

Atomic  Mass  (amu) 

120 

119.94288 

122 

121.94193 

123 

122.9^368 

124 

123.94278 

125 

124.94460 

126 

125.94420 

128 

127.94649 

130 

129.94853 

E.  Advantages  of  looking  at  sulphur,  selenium,  and  tellurium 
isotopes  simultaneously  in  natural  deposits 

One  of  the  reasons  for  considering  tellurium  for 
an  isotope  fractionation  study  is  to  compare  its  behaviour 
with  results  already  obtained  with  sulphur  and  selenium.  As 
it  has  already  been  shown,  a  considerable  number  of  studies 

O  o  Q 

of  the  SJ  /SJ  variations  have  yielded  much  information  about 
natural  processes.  The  preliminary  study  which  has  been  con¬ 
ducted  with  Se^/Se^  also  looks  promising  in  that  it  shows 
some  similarities  to  the  sulphur  isotope  distribution. 

Whereas  sulphur  exists  in  nature  from  the  -2  to 
the  +6  valence  states,  it  is  very  rare  that  selenium  exists 
in  valence  states  higher  than  +4,  and  tellurium  rarely  exists 


Absolute  Abundances  of  the  Isotopes  of  Tellurium 


Table  II  cont'd 
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in  a  valence  state  as  high  as  +4.  This  results  because 
of  the  different  oxidation  potentials  for  the  elements. 
Further,  similar  chemical  reactions  of  these  elements  have 
different  temperature  dependences,  and  because  of  their 
different  sizes,  the  three  elements  behave  differently 
in  physical  processes  such  as  diffusion. 

If  all  three  elements  can  be  simultaneously  ex¬ 
amined  isotopically,  it  may  be  possible  to  separate  the 
various  physical,  chemical,  and  biological  processes  which 
have  occurred  in  a  deposit. 


THEORY 


An  isotope  effect*  or  isotope  fractionation*  occurs 
when  the  relative  abundances  of  the  isotopes  of  an  element 

in  some  system  are  altered  by  physical  or  chemical  means* 
and  is  a  result  of  slight  differences  in  the  properties  of 
these  isotopes.  Chemically  an  effect  may  occur  when  two 
or  more  molecular  species  containing  the  element  in  question 
co-exist  in  chemical  equilibrium  (equilibrium  isotope  effect) 
or  when  one  or  more  molecular  species  change  their  identity 
in  the  course  of  a  chemical  reaction  (kinetic  isotope  effect) 

The  principles  of  isotopic  fractionation  are  now 
well-known  and  it  is  possible  to  calculate*  from  a  knowledge 
of  the  fundamental  vibrational  frequencies  of  the  compounds* 
the  extent  to  which  the  isotopes  might  be  expected  to  frac¬ 
tionate  in  a  naturally  occurring  process. 

Looking  at  the  members  of  the  group  VIA  elements* 
there  is  a  mass  difference  of  two  or  a  percentage  mass 
difference  of  13$  in  the  case  of  the  two  oxygen  Isotopes 
0  D  and  0  ;  a  mass  difference  of  two  or  a  percentage  mass 

difference  of  7$  in  the  case  of  S  and  S  ;  a  mass  differ¬ 
ence  of  six  or  percentage  mass  difference  of  8$  in  the  case 

of  Se^°  and  Se^'c;  and  a  mass  difference  of  eight  or  percentag 

130  122 

mass  difference  of  7$  in  the  case  of  Te  ^  and  Te 
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Significant  isotope  effects  have  been  found  for  the  two 
oxygen,  sulphur,  and  selenium  isotopes  respectively,  and 
therefore,  some  fractionation  of  the  tellurium  isotopes 
is  expected. 


A.  Theory  of  Equilibrium  Isotope  Effects 

Although  chemical  equilibrium  may  be  considered  as 
that  condition  of  a  chemical  system  where  the  rate  at 
which  a  forward  reaction  is  taking  place  is  exactly  balanced 
by  the  rate  of  the  reverse  reaction,  the  theory  of  equili¬ 
brium  isotope  effects  is  best  considered  in  terms  of  the 
statistical  mechanics  of  a  system  of  particles  in  dynamic 
equilibrium. 

Urey  and  Rittenberg  (ll4),  (115)*  were  the  first 
to  show  by  theoretical  calculations  that  marked  differences 
in  the  equilibrium  constants  of  isotopic  exchange  reactions 
should  exist.  Urey  and  Grieff  (ll6)  applied  statistical 
mechanics  to  the  calculation  of  equilibrium  exchange  con¬ 
stants  in  terms  of  distribution  functions  and  "summation 
over  states"  terms  and  this  was  simplified  by  Urey  (ll8) 
and  Bigeleisen  and  Mayer  (10).  Their  simplification  made 
it  possible  to  calculate  the  equilibrium  constants  of 
isotopic  reactions  solely  with  a  knowledge  of  the  vibrational 
frequencies  of  the  isotopic  molecules.  It  is  to  be  noted. 
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however,,  that  such  a  treatment  does  not  consider  the  time 
taken  to  attain  equilibrium  or  whether  there  in  fact  exists 
a  mechanism  by  which  equilibrium  may  be  attained. 

An  isotopic  exchange  reaction  may  be  written 

aA^  +  bB2  ^  aA2  +  bB^  (l) 

where  the  subscripts  1  and  2  refer  to  the  light  and  heavy 
isotopes  respectively  and  A  and  B  are  molecules  having  the 
element  under  consideration  as  a  common  constituent.  From 
thermodynamics  (72) ,  the  standard  free  energy  change  for 
the  above  reaction  is  given  in  terms  of  the  equilibrium 
constant  by 

-R'T  InK  =  AF°  (2) 

or 

-RT  In  K  =  aFA  +  bF^  -  aFA  -  bF^  (3) 

Ar)  -D-j-  A-^  -lV} 

In  terms  of  the  partition  function  "Q",  the  free 
energy  can  be  written  as  (F  =  Gibbs  free  energy) 

F  =  E  +  RT  InN  -  RT  InQ  (4) 

o 

where  N  is  AvagadroTs  number,  Eq  is  the  zero  point  energy, 
and  Q  is  the  total  partition  function  given  by 

-e  /kT 

Q  =  2  gne 

n 


(5) 
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where  gn  is  the  statistical  weight  factor  and  the  summation 
extends  over  all  the  quantum  states  of  energy  €  .  A  d-fold 
degenerate  level  would  count  as  d  states. 

If  only  one  molecule  is  considered  in  unit  volume 
instead  of  one  mole,  N  is  replaced  by  unity  and  the  second 
term  of  (4)  disappears.  Substitution  of  (4)  into  (3)  and 
simplification  gives 
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Since  they  are  determined  by  the  electrical  forces,  the 
potential  energy  curves  have  been  found  to  be  essentially 
identical  for  isotopic  molecules.  Referring  all  energies 
to  the  minimum  of  the  potential  energy  curves  rather  than 
the  zero  point  energy  simplifies  (6)  to 
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and  the  equilibrium  constant  becomes  the  ratios  of  the 
partition  functions  for  the  two  isotopically  substituted 


molecules . 
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B.  Evaluation  of  Partition  Function  Ratios 

The  partition  function  of  a  molecule  depends  upon 
its  mass,  moments  of  inertia,  fundamental  vibrational 
frequencies,  temperature,  and  nuclear  spins  of  the  consti- 
uent  atoms.  It  is,  however,  possible  to  evaluate  parti¬ 
tion  function  ratios  for  isotopic  molecules  using  only- 
vibrational  frequencies. 

Following  Herzberg  (49)  considering  a  volume  V, 


3/2 

^  (2irMkT) 

^  —  o 


3n-5  -u  /2 


V  .  kT 

i — 1 1 

J 

•H 

t 

CD 

si _ 

(21+1) 

hcB 

i 

-U.  0 

1-e 

(21+1) 

2 


2 


for  diatomic  (or  linear  polyatomic)  molecules,  and 
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for  polyatomic  molecules,  where 

n  =  number  of  atoms  in  the  molecule 
M  =  total  molecular  mass 
k  =  Boltsman  constant 

A,B,  C  =  rotational  comstants  of  the  molecules,  related 
to  its  principal  moments  of  inertia,  i.e.. 
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hv . 

4“  In 

Ui  =  ~kT  where  v  is  the  i  fundamental 


frequency  of  the  molecule. 
The  product 
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extends  over  all  i  fundamental 
frequencies  of  the  polyatomic  molecule  and  the  bracketed 
terms  are  symmetry  and  nuclear  spin  factors.  The  nuclear 
spin  factors  cancel  out  in  later  expressions  and  have  no 
effect  on  the  estimation  of  isotope  effects.  The  symmetry 
factors  are  discussed  later. 


Thus,  for  two  isotopic  polyatomic  molecules. 
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Urey  (ll8)  simplified  the  two  expressions  (  8  ) 

and  (  9  )  and  defined  new  partition  functions  which  were 

the  equilibrium  constants  for  exchange  reactions  between  the 
compound  considered  and  the  separated  atoms.  Both  sides  of 
equations  (  8  )  and  (  9  )  are  multiplied  by  (M^/M^)  ^//^n 
where  M-,  and  M0  are  the  atomic  weights  of  the  isotopic 
atoms  being  considered  and  n  is  the  number  of  the  isotopic 
atoms  exchanged.  Further,  the  right  sides  of  (  Q  )  and 
(  9  )  are  multiplied  and  divided  by  (u-j/u2)  and 


(u1./u21) 

respectively.  According  to  the  Teller-Redli ch  theorem  (70,87)* 


1 


47 


I2 

N 

3/2 

Til 

3/2n 

U1 

^2^2^2 

1/2 

N 

3/2 

Til 

3/2n  I 

ui.n 

h 

LhJ 

N 

/2_ 

LAiBTiJ 

KJ 

LM2j 

i 

i 

Lu2<  J 

i 


Carrying  out  the  steps  above  and  using  (  11),  equation 

(  10  )  becomes 
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for  diatomic  molecules  and  polyatomic  molecules  respectively. 


The  terms  (M^  /M-^  will  not  affect  the 

i  i 

expression  for  an  equilibrium  constant  since  they  will  appear 
in  both  numerator  and  denominator.  Thus  the  equilibrium 
constant  can  be  given  by 
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is  the  relation  between  the 
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new  partition  function  ratio  and  the  old  ratio.  K  can  now 
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be  calculated  with  only  a  knowledge  of  the  vibrational 
frequencies  of  the  isotopic  molecules.  These  frequencies 
are  obtained  either  by  direct  observation  or  by  calculation 
from  observations  on  a  more  abundant  isotope  by  means  of 
force  equations. 

This  expression  excludes  the  anharmonic  terms  in 
the  vibrational  energy.,  and  assumes  that  the  rotational 
partition  functions  have  reached  classical  values.  However, 
these  have  been  investigated  by  Urey  (118)  and  Herzberg  (49) 
and  the  approximation  is  not  greatly  in  error  when  consider¬ 
ing  the  ratios  of  partition  functions  for  isotopic  molecules, 
provided  the  observed  fundamental  frequencies  are  used  rather 
than  the  zero  order  frequencies  and  that  the  temperature  is 
low  enough  so  that  only  low  lying  vibrational  states  have 
appreciable  populations.  The  ratio  of  symmetry  numbers 
will  be  unity  if  the  molecule  under  consideration  contains 
only  one  atom  of  the  element  for  which  an  exchange  is  con¬ 
sidered  or  if  the  molecule  contains  more  than  one  such  atom 
but  these  atoms  occupy  indistinguishable  positions  in  the 
molecule  and  are  all  exchanged  in  the  reaction. 

Thus,  we  have  for  diatomic  molecules 
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and  for  polyatomic  molecules, 
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These  expressions  have  been  simplified  by  Urey 
(ll8)  and  Bigeleisen  and  Mayer  (10)  for  manual  calculations. 
In  the  work  of  this  thesis  IBM  7040  programs  were  set  up 
to  calculate  the  exact  expressions  directly. 

C.  Kinetic  Isotope  Effects 

The  theory  of  isotope  fractionation  in  kinetic 
chemical  reactions  is  considered  in  terms  of  the  theory  of 
absolute  reaction  rates.  This  theory  was  developed  by 
Eyring  (36)  and  Evans  and  Polanyi  (35)  from  a  statistical 
treatment  of  reaction  rates.  Bigeleisen  (9)  then  used 
their  theory  and  collision  theory  to  develop  formulae 
which  predict  the  kinetic  isotope  effect. 

This  theory  is  developed  around  two  postulates. 
Firstly,  there  exists  an  "activated  complex"  in  chemical 
equilibrium  with  the  initial  state  of  the  reactants  and 
intermediate  between  this  initial  state  and  a  final  state 
of  the  products.  Secondly,  this  "activated  complex"  de¬ 
composes  at  a  definite  rate  to  form  the  products  of  the 


reaction . 
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From  these  postulates,  the  rate  constant  k  can 
be  written  (40a)  as 


k '  T 

k  =  -fp  .  kTtK 


where 

k’  =  Boltzman  constant 
T  =  absolute  temperature 
h  =  Planck’s  constant 
k,p  =  transmission  coefficient 
t  =  quantum  mechanical  tunnelling  correction 
K  =  equilibrium  constant  for  the  equilibrium  between 
the  reactant  and  the  activated  complex 


Thus,  for  isotopic  molecules,  the  ratio  of  rate  constants 
becomes , 


k 


T1 


k 


T2 


Hirschfelder  and  Wigner  (52)  have  conducted  a 
theoretical  study  to  show  that  the  transmission  coeffi¬ 
cients  ratio  is  very  nearly  unity  for  systems  above  room 
temperature  having  a  distribution  in  the  velocities  of  the 
reacting  isotopic  molecules.  Also,  t^  is  set  equal  to  t^, 
the  assumption  of  which  is  discussed  by  Bigeleisen  ( 1 3 ) • 


Thus  we  can  write 


51 


Prom  the  previous  discussion  of  equilibrium 

constants , 


where  Q,  and  Qp  refer  to  partition  functions  for  the 

$  t 

reactants,,  and  Q  and  Qp  refer  to  partition  functions 
for  the  activated  complexes.  The  partition  function  ratios 
can  be  set  up  in  the  same  manner  as  was  done  in  part  A 
of  the  theory.  For  the  activated  complex,  however,  we 
have  to  consider  its  nature.  According  to  the  theory,  one 
of  the  vibrational  degrees  of  freedom  of  the  activated 
complex  is  non-genuine,  corresponding  to  the  cleavage 
(or  formation)  of  a  bond.  Thus,  supposing  that  the  L 
vibrational  mode  is  non-genuine,  so  that  has  a  pure 
imaginary  value,  the  vibrational  contribution  to  the  par¬ 
tition  function  of  the  activated  complex  will  then  take 
the  form 
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so  that 
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for  polyatomic  molecules. 


Again,  the  mass  ratios  will  cancel  so  we  can  write 
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Furthermore,  this  "imaginary"  frequency  is  the 
stretching  frequency  of  the  bond  being  broken  (or  formed). 
Thus,  this  frequency  can  be  expressed  as  the  square  root  of 
the  inverse  of  the  effective  masses  along  the  co-ordinate 
of  decomposition  of  the  activated  complex.  Our  expression 


then  becomes 


53 


Q, 


t' 


Q 


t 


1 


tJ  1 


M 


M 


t 


1/2 


3n^-7 


i/L 


U 


X 


-u 


21 


t 

21/2 


U 


* 


11 


-U 


X 

11/2 


"Uli 

e  11 ) 


(1  - 


-°L 

e  ) 


writing  this  as 
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it  can  be  seen  that  equation  (17) 


can  be  written  as 


(18) 


Prom  equation  (l8)  it  can  be  seen  that  the  light 
molecule  usually  has  a  greater  rate  constant  since  the 

4-4-  . 

factor  (Mg/lYLp  '  is  always  greater  than  unity  and  the 
reacting  molecule  is  usually  more  tightly  bound  than  the 
activated  complex.  That  is 
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Regarding  the  ratio  (M^/mJ),  two  different 
approaches  can  be  made.  Wolfsberg  (136)  gives  the 
following  reasoning.  Consider  a  molecule,  X-A-B-Y  in 
which  the  A-B  bond  is  broken.  If  the  A-B  bond  is  strong 
compared  to  the  X-A,  B-Y  bonds,  then  the  latter  will 
shorten  as  the  former  lengthens  and  the  reduced  mass  will 
be  closer  to  (^rg)  •  If,  however,  the  A-B  bond  is  the 
weaker,  then  the  reduced  mass  will  be  closer  to 

( ^XA * Mxb) / ( Mxa  +  MXB^  since  m&ss  centres  and 
are  being  separated. 

Since  very  little  is  known  about  the  nature  of 
the  activated  complex,  the  rate  constant  cannot  be  evaluated 
specifically.  However,  three  different  approaches  are 
feasible  which  follow  the  formalism  of  the  theory  and 
incorporate  the  idea  that  in  the  transition  state,  the 
vibrational  frequencies  of  the  activated  complex  are  the 
same  as  those  of  the  reactant  molecule  except  for  that 
corresponding  to  bond  cleavage  or  formation.  Approaches  I 
and  III  give  extreme  lower  and  upper  limits  respectively 
that  might  be  expected  for  the  ratio  of  rate  constants  and 
approach  II  gives  an  intermediate  and  best  description  of 
the  physical  situation. 
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I.  When  the  partition  function  ratio  term  is 
very  close  to  unity  which  occurs  when  the  temperature  is 
high.  In  the  high  temperature  limit. 
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M2 
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The  mass  ratio  can  either  be  the  atomic  mass  ratio  or 
mass  fragment  ratio. 


II.  Where  the  activated  complex  is  identical  to 
the  starting  material,  with  the  exception  that  the  fre¬ 
quency  associated  with  the  breaking  bond  is  missing.  In 
this  case. 
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where  this  ratio  is  the  same  as  the  partition  function 
ratio  for  diatomic  isotopic  molecules  with  frequencies 

an<^  ^11/  that  frequency  in  the  starting  material 

which  is  missing  in  the  activated  complex.  The  mass  ratio 
is  evaluated  as  before. 


III.  Where  the  activated  complex  has  some  struc¬ 
ture  which  is  strictly  speaking  unknown,  but  for  which 
the  value  of  the  partition  function  ratio  lies  between 
that  for  the  starting  material  (with  the  omission  of  a 
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frequency  term)  and  unity.  The  extreme  case  of  unity 
may  be  approached  if  in  the  transition  state  there  is  a 
complexing  of  the  reacting  molecule  with  other  material , 
so  that  a  complicated  structure  is  formed  in  which  not 
many  of  the  vibrational  frequencies  are  sensitive  to 
isotopic  substitution.  In  this  case. 


«1 


1/2 


These  three  approaches  give  for  a  particular 
model  of  the  activated  complex  an  intermediate  (II )  as 
well  as  extreme  (I,  and  III)  estimates  for  the  ratio 
of  rate  constants. 

D.  Isotope  fractionation  during  first  order  competing 
reactions 

In  general,  we  have  the  competitive  isotopic 
reactions  of  the  type 

kl 

A1  +  B  - - - >  X1  +  Y 

k2 

A2  +  B  - - - *  X2  +  Y 

where  subscripts  1  and  2  refer  to  the  light  and  heavy 
isotopic  species  respectively.  The  rate  constant  k 
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indicates  how  fast  A  is  converted  to  X. 

The  rate  of  the  reaction  is  given  in  terms  of 
the  rate  constant  as 


Rate  =  kC^Cg 


where  k  =  the  rate  constant 

=  concentration  of  reactant  A 
Cg  =  concentration  of  reactant  B 

A  first  order  reaction  is  one  where  the  time  dependence 
of  product  formation  can  be  written  as 


dX 

dt 


=  k(A0  -  X)B 


where  A^  is  the  initial  concentration  of  the  reactant. 

Thus,  for  the  competing  reactions  mentioned. 


dX 


dt 


1  =  k1(A01  -  X1)B 


dX 


dt 


2  =  k2(AQ2  -  X2)B 


At 


t  =  0,  X1  =  X2  =  0,  so  upon  integrating. 


k 


log 


1 


A 


01 


A01  "  X1 


k. 


A 


02 


A02  "  X2 


log 
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If  we  designate  the  fraction  of  molecules  which  have 
reacted  as 


X2  +  x2  _  X1(l  +  x2/x1) 

A01  +  A02 


and  the  integrated  isotopic  composition  of  the  product  as 

x2Ax 

A02/A01 


we  can  write 


log 
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1  +  A02//A01 

1  +  x2/x~ 


log 


r  f 


1  +  A02//A01 

i  +  x2/x1 


(19) 


Thus,  the  ratio  of  the  rate  constants  for 
competitive  isotopic  reactions  is  obtained  from  the 
determination  of  isotopic  ratios  in  the  initial  reactant 

A02//A01  and  in  the  product  X^/X^  after  fraction  f  of  the 

molecules  have  reacted. 


THEORETICAL  CALCULATIONS 


A.  Calculation  of  equilibrium  exchange  constants 

As  was  noted,  the  extent  to  which  one  would 

expect  Isotopic  fractionation  to  take  place  in  exchange 

reactions  is  found  from  the  partition  function  ratios. 

Thus,  the  partition  function  ratios  have  been  calculated 

ISO  122 

for  various  molecules  containing  Te  and  Te 

Appendix  I  contains  calculations  of  vibrational 

ISO  122 

frequencies  for  some  Te  and  Te  -containing  compounds 
for  which  molecular  data  such  as  symmetry  type  and  fun¬ 
damental  vibrational  frequencies  were  available.  For 

the  diatomic  species,  anharmonicity  corrections  were  made. 

ISO 

Since  Te  is  the  most  abundant  isotope,  and  also  since 

use  of  this  mass  number  gave  results  closest  to  the  observed, 

it  was  assumed  that  fundamental  frequencies  applied  to 

compounds  containing  this  isotope.  The  frequencies  of 
122 

the  Te  -containing  molecules  were  then  calculated  by 
means  of  force  field  equations  from  these  observed  values. 

In  most  cases,  sets  of  equations  describing 
the  force  fields  were  accompanied  by  force  constants 
already  calculated.  However,  in  the  case  of  H^Te  and  TeO^ 
the  force  constants  had  not  been  previously  calculated  so 
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these  were  calculated  by  substituting  the  known  funda¬ 
mental  frequencies  into  a  set  of  force  equations  describing 
the  molecule . 

From  the  isotopic  vibrational  frequencies,  the 
partition  function  ratios  at  various  temperatures  were 
calculated  using  equations  (15)  and  (l6)  for  diatomic  and 
polyatomic  molecules  respectively.  The  equilibrium 
exchange  constants  for  possible  isotopic  exchange  reactions 
were  then  calculated  using  equation  (l4).  Tables  III  and 
IV  give  the  partition  function  ratios  and  equilibrium 
exchange  constants  respectively.  Looking  at  Table  IV, 
the  values  in  the  body  of  the  table  are  the  equilibrium 
constants  for  exchange  reactions  between  the  compounds 
listed  along  the  top  and  the  compounds  listed  down  the 
left-hand  column  calculated  for  the  temperatures  listed 
in  the  far  right-hand  column.  A  constant  greater  than 
unity  means  that  the  heavier  isotope  is  concentrated  in 
the  compound  in  the  left-hand  column. 

For  example,  consider  TeFg  in  the  left-hand 
column,  and  TeO  at  the  top  of  the  table.  The  value  1.023^ 
is  found  at  20°C  in  the  table  where  the  TeO  column  and 
TeF^  row  intersect.  This  is  the  equilibrium  constant  for 
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Te13°0  +  Te122F£  zJL.  Te1220  +  Te13°F^ 

This  means  that  provided  the  reaction  is  possible,  a  2.3 % 
tellurium  isotope  effect  is  predicted  and  the  TeF^ 
species  will  be  enriched  by  2.3 %  in  the  heavier  isotope. 


B.  Cal culation  of  ratios  of  isotopic  rate  constants 

As  mentioned  in  the  theory,  intermediate  values 
as  well  as  upper  and  lower  extremes  could  be  determined 
for  the  ratio  of  kinetic  rate  constants.  The  lower  limit 
reached  at  high  temperatures  is  where 


t 

M ^  n 

T 

M1 


1/2 


4-4.  . 

and  (M^/M^)  '  represents  the  effective  mass  of  the 

complex  along  the  reaction  co-ordinate.  If  this  is  simply 
the  reduced  mass  of  the  atoms  whose  bond  is  directly 

t  t 

involved  in  the  reaction.  Mg  and  a.re  the  isotopic 

reduced  masses  of  the  molecules.  Table  V  gives  the  mini¬ 
mum  value  of  the  rate  constant  ratio  to  be  expected  in  a 
kinetic  isotope  reaction  in  which  the  bond  is  broken  between 
the  diatomic  molecules  listed.  These  are  temperature 
independent . 
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Minimum  Value  of  k-^Q 

/k122 

GeTe 

PbTe 

SnTe 

TeO 

SiTe 

Te  Te 

HTe 

1.012 

1.020 

1.016 

1.004 

1.006 

1.016 

1.000 

Table  V 


Considering  now  the  reduction  of  potassium 
tellurite,,  and  suppose  the  rate  determining  step  is  the 
cleavage  of  a  Te  -  0  bond: 

I.  The  high  temperature  limit  for  k^/k^  is  the  mass 

term  1.004  if  atom  fragments  are  considered  or  1.002 
if  mass  fragments  in  the  "di-fragmental  molecule" 
TeO^  -  0  are  used. 

II.  The  intermediate  value^and  closest  to  the  physical 


+ 1  /  -f  *  = 

/  Q-^  is  simply  that  for  TeO^  with  the 

frequency  contribution  of  the  Te  -  0  bond  omitted. 

This  mode  is  most  probably  one  component  of  the 

doubly  degenerate  E  type  asymmetric  stretching 
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1.0041  at  0°C 
1.0037  at  20°C 
1.0031  at  50°C 


The  mass  ratio  is  as  determined  in  part  I  for  either  atomic 
or  mass  fragments. 


III.  The  upper  limit  for  k-^/k^  is  given  by 


kx  Q(Te13°0~ 
k2  Q(Te1220p 


The  results  of  the  three  methods  of  calculation 
are  summarized  in  Table  VI 0  All  these  simple  models  give 
a  regular  decrease  in  the  value  of  k^/k^  with  increase  in 
temperature . 
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EXPERIMENTAL 


In  investigating  the  isotopes  of  tellurium,  tellurium 
samples  were  examined  to  see  if  the  isotope  ratio  could 
be  altered  by  natural  and  laboratory  processes.  Tellurium 
was  extracted  from  natural  telluride  samples  taken  from 
various  geographical  locations,  A  controlled  chemical 
reduction  of  K^TeO^  to  elemental  tellurium  was  carried 
out  in  the  laboratory  and  attempts  made  to  measure  the 
kinetic  isotope  effect,  K^TeO^  was  also  partially  re¬ 
duced  to  elemental  tellurium  by  micro-organisms  and  the 
isotope  composition  of  the  products  and  reactant  compared. 

The  collected  elemental  tellurium  samples  were  pre¬ 
pared  for  mass  spectrometric  analysis  by  converting  them 
to  tellurium  hexafluoride  (TeP^) .  This  method  was  cho¬ 
sen  because  of  the  simplicity  of  the  TeF^  spectrum  (there 
is  only  one  stable  fluorine  isotope),  the  stability  of 
the  gas,  and  because  of  its  ease  of  preparation.  HgTe 
could  also  be  used,  but  the  gas  is  unstable  and  its  mass 
spectrum  is  quite  complicated  by  overlapping  of  the  ion 
species,  eg.,  Te1^  and  Te^^^H^  both  at  mass  130.  Com¬ 
putations  on  such  a  spectrum  are  tedious  and  subject  to 


considerable  error. 
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A .  Extraction  of  Tellurium  from  Natural  Samples 

In  1832  it  was  discovered  by  P.  Berthier  that 
tellurium  could  be  precipitated  from  raw  materials  with 
sulphur  dioxide  and  hydrochloric  acid  (76) .  However, 
the  extraction  of  pure  elemental  tellurium  was  hampered 
by  the  presence  of  selenium  which  is  closely  related  to 
tellurium  being  of  the  same  chemical  group.  Keller  (58) 
in  1897  showed  that  selenium  and  tellurium  could  be 
separated  with  sulphur  dioxide  in  differing  acidities 
of  hydrochloric  acid,  and  Gooch  and  Peirce  in  1898  (42) 
were  able  to  separate  tellurium  from  selenium  on  the 
basis  of  the  different  volatilities  of  their  tetra- 
bromides,  but  this  was  somewhat  inefficient. 

To  further  add  to  the  difficulty,  certain 
heavy  metals  are  partially  precipitated  by  sulphur 
dioxide  (58),  and  combinations  of  bismuth,  antimony, 
and  selenium  with  tellurium  make  separations  even 
more  complicated  because  of  mutual  interferences  (88) . 
Hillebrand  and  co-workers  in  1953  (51)  gave  the  dif¬ 
ferent  hydrochloric  acid  concentrations  for  pre¬ 
cipitation  of  each  of  selenium  and  tellurium  with  sulphurous 
acid.  However,  they  dried  both  elements  in  air  at 
110°C;  whereas,  in  the  same  year  Duval  showed  that  tel- 
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lurium  oxidizes  at  temperatures  above  4o°C  (30).  Then 
in  I960,  a  method  was  described  by  Reed  (88)  by  which 
pure  elemental  tellurium  could  be  separated  from  all 
other  elements.  This  method  was  employed  in  this  in¬ 
vestigation  . 

About  0.5  grams  of  sample  are  powdered  and  dis¬ 
solved  in  10  ml,  of  nitric  acid  and  evaporated  to  dry¬ 
ness,  The  residue  is  dissolved  in  100  ml,  of  concen¬ 
trated  hydrochloric  acid  with  25  ml,  sulphurous  acid. 
This  step  precipitates  any  selenium  present,  which  is 
then  filtered  out. 

To  the  filtrate,  2  grams  of  tartaric  acid  are 
added  to  prevent  any  antimony  from  hydrolyzing.  Water 
is  then  added  until  the  hydrochloric  acid  concentration 
is  3M.  To  speed  up  the  reaction,  10  ml,  of  hydrazine 
hydrochloride  solution  is  added  and  the  tellurium  is 
then  precipitated  with  25  ml,  of  sulphurous  acid. 

After  boiling  for  1  or  2  minutes,  the  tellurium  is  sep¬ 
arated  from  the  supernatant  solution  by  centrifuging. 
This  is  then  washed  and  re-centrifuged  several  times 
to  remove  soluble  ions. 

The  method  works  as  follows:  the  nitric  acid  acts 
as  an  oxidizing  agent  producing  H^TeO^  . 


The  sulphurous 
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acid  then  reacts  according  to  the  equation 

H2Te03  +  2H2S03  3M  HC1-  Te  +  2R2S0^  +  H20 

(The  selenium  was  precipitated  in  like  manner  in  the 
concentrated  hydrochloric  acid  solution.) 

According  to  the  literature,  only  gold  and  silver 
would  interfere  in  the  selenium-tellurium  determination 
if  present,  but  the  gold  not  being  dissolved  in  nitric 
acid  can  be  filtered  out,  and  any  silver  is  easily 
separated  as  chloride  (58). 

The  efficiency  of  this  method  was  tested.  The 
precipitated  tellurium  was  filtered  through  a  fritted 
crucible,  washed  with  0.2N  hydrochloric  acid  and  finally 
water.  It  was  dried  by  placing  the  crucible  and  tellur¬ 
ium  in  a  desiccator  at  20  mm.  pressure  or  less  at  room 
temperature  over  anhydrous  calcium  sulphate  for  two  hours. 
Of  13^6.5  mg.  of  sample,  1323.3  mg.  or  98.28 %  was  re¬ 
covered  . 

B.  Chemical  reduction  of  K^TeO^ 

In  order  to  study  kinetic  isotope  effects,  the 
desireable  conditions  for  the  experiment  are: 

1.  a  gaseous  reaction  or  reaction  in  solution. 
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2.  a  method  of  halting  the  reaction  at  any 
point,  and 

3.  the  ability  to  recover  the  reduced  product. 
A  solution  or  gaseous  reaction  is  preferable  since  larger 
isotope  effects  are  found  under  these  conditions.  Very 
little  isotope  fractionation  is  expected  in  a  solid  phase 
reaction  where  the  process  proceeds  layer  by  layer. 

The  second  condition  is  necessary  if  the  kinetic  effect 
is  to  be  measured.  The  reaction  must  be  stopped  and  the 
isotope  ratios  determined  in  the  unreduced  reactant  and 
the  reduced  product.  Therefore,  a  chemical  reduction  of 
K^TeO^  in  solution  was  chosen  for  study. 

Halting  the  reaction  is  possible  in  either  of  two 
ways.  A  method  of  complete  reduction  can  be  employed 
and  the  reaction  quickly  halted  by  neutralization  at  the 
required  time;  or,  only  a  percentage  of  the  reducing 
agent  required  for  complete  reduction  could  be  used. 
Halting  by  neutralization  has  disadvantages  in  that 
much  time  and  many  trial  runs  are  required  to  determine 
the  exact  time  to  be  allowed  for  the  desired  reduction 
to  take  place. 

Furthermore,  in  the  case  of  tellurium,  using  the 
hydrochloric  acid-sulphur  dioxide  method  for  instance,  the 
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reaction  takes  several  days  for  complete  reduction  un¬ 
less  hydrazine-hydrochloride  is  used,  in  which  case  the 
reduction  takes  place  too  rapidly  to  control.  Therefore, 
5$  -  10$  reductions  were  affected  by  adding  only  5$  - 
10$  of  the  stoichiometrically  required  reducing  agent. 

Using  the  reduction  equation 

K2Te03  +  2H2S03  - »  K2SC>4  +  H2S04  +  Te  +  H20 

253.81  164.03  174  98.02  127.6  18.02 

0.65  grams  of  H2S03  are  required  for  the  complete  re¬ 
duction  of  one  gram  of  K2Te03#  For  a  10$  reduction, 

10$  of  O.65  grams  or  65  mg.  are  required  per  gram  of 
K2Te03 .  Six  grams  of  K2Te03  were  reduced  10$  at  room 
temperature . 

The  precipitated  tellurium  was  separated  from  the 
supernatant  solution  by  centrifuging,  was  washed,  re¬ 
centrifuged  to  remove  soluble  ions,  and  then  dried  at 


C.  Reduction  of  K2TeC>3  by  Scopulariopsis  brevlcaulis 
That  micro-organisms  are  capable  of  reducing 
a  complex  tellurium  compound  to  elemental  tellurium  has 
been  known  for  several  years.  As  early  as  1914,  W.  E. 

King  and  L.  Davis  observed  that  nearly  all  of  the  more 
common  micro-organisms  reacted  with  potassium  tellurite 
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to  form  "black  compounds"  depending  on  reduction  of  the 
tellurite  (59).  T.  L.  Ch’in  in  1938  discovered  that  living 
colonies  of  pathogenic  fungi  could  be  identified  by  the 
"black  color"  produced  upon  injection  of  potassium  tellurite 
(l8).  A  year  later.  Bird  and  Challenger  (13)  discovered 
that  upon  ingestion  of  potassium  tellurite,  Scopulariopsis 
brevicaulis  (Penicillium  brevicaule)  reduced  the  tellurite 
to  black  elemental  tellurium  and  evolved  di -methyl  telluride. 
Other  micro-organisms  were  later  discovered  to  reduce  potas¬ 
sium  tellurite  to  elemental  tellurium  such  as  Coryne bacterium 
diphtheria  by  Morton  and  Anderson  (79)  in  1941,  Myco¬ 
bacterium  tuberculosis  by  L.  Sula  in  1947  (102),  and  living 
tissues  (bacteria,  yeast,  human  leucocytes)  by  M.  Wachstein 
in  1949  (125).  Others  can  also  be  found  in  later  liter¬ 
ature  . 

To  examine  the  isotopic  variations  by  such  a  reduction, 
it  was  decided  that  Scopulariopsis  brevicaulis  should  be 
used  because  it  is  non-pathogenic .  Also,  a  culture  of  the 
organisms  was  readily  obtainable  from  Dr.  Carmichael  of 
The  Provincial  Laboratory  of  Public  Health,  University  of 
Alberta.  Furthermore,  growth  of  the  culture  is  quite 
simple  and  may  be  carried  out  at  room  temperature. 


75 


1 .  Reduction  of  K^TeO^  by  culture  on  liquid  medium 
The  culture  was  grown  on  a  liquid  medium 
(Trypticase  Soy  Broth)*  in  an  aerating  flask.  500  ml. 
of  TSB  was  innoculated  with  S.  brevicaulis .  The  air  was 
sterilized  by  passing  it  through  H^SO^,  sterile  cotton 
wool,  then  a  1:1000  solution  of  mercuric  chloride,  and 
again  through  sterile  cotton  wool  (17).  After  passing 
through  the  culture  flask,  it  was  aspired  through  an 
alcoholic  iodine  solution  to  form  di -methyl  telluride  di¬ 
iodide,  (l6)  from  any  di -methyl  telluride 

evolved  from  the  culture. 

The  culture  was  allowed  to  grow  until  a  good  growth 
covered  the  medium  surface  (12  days).  It  was  then  in¬ 
noculated  with  1.5  grams  of  potassium  tellurite  in  solution, 
and  samples  were  taken  every  12  hours  until  enough  of  the 
tellurite  had  been  reduced  to  yield  sufficient  tellurium 
for  analysis,  which  took  72  hours.  At  this  point,  10  cc. 
samples  were  taken  every  6  hours  and  the  tellurium  ex¬ 
tracted  from  the  broth  as  outlined  later.  Finally,  after 
100  hours,  the  yield  of  reduced  tellurium  was  too  small 
for  analysis  and  the  experiment  was  halted.  The  tellurium 
was  then  extracted  from  the  yet  unreduced  tellurite,  the 
broth,  and  from  the  di -methyl  telluride  di -iodide  compound. 

*This  broth  was  prepared  by  mixing  30  grams  Trypticase  Soy 
Broth  powder  from  Baltimore  Biological  Company  with  1  litre 
of  water.  Hitherto  this  will  be  referred  to  as  TSB. 


2.  Extraction  of  elemental  tellurium  reduced  from 


K2Te03  in_  the  biological  investigation 

In  order  to  obtain  only  the  tellurium  reduced 
by  the  bacteria  for  analysis,  it  had  to  be  separated  from 
the  rest  of  the  broth  which  still  contained  unreduced 
potassium  tellurite  in  solution  and  which  would  also  be 
reduced  by  the  extraction  procedure.  Thus,  the  10  cc. 
sample  was  washed  with  water  and  filtered  to  separate  the 
tellurium  reduced  by.  the  bacteria  from  the  unreduced  tel¬ 
lurite  solution.  (Elemental  tellurium  is  insoluble  in 
water.)  The  filter  paper  was  then  dissolved  in  nitric  acid 
and  the  method  of  extraction  carried  out  as  outlined  on  page 
72.  The  filtrate  was  treated  with  concentrated  hydrochloric 
acid  and  sulphur  dioxide  and  let  stand  for  three  hours. 

After  this,  the  hydrochloric  acid  concentration  was  brought 
to  3M,  10  ml.  of  hydrazine  hydrochloride  was  added,  and 
sulphur  dioxide  was  bubbled  through  which  precipitated  any 
tellurium  present.  Thus,  at  any  sampling  time  the  bac- 
teriologically-reduced  tellurium  and  tellurium  represent¬ 
ative  of  the  unreduced  tellurite  were  available  for  analysis. 

The  tellurium  was  separated  from  the  di -methyl  telluride 
di-iodide  by  a  method  outlined  by  Vernon  (122).  Tellurium 
and  methyl  iodide  react  according  to  the  equation 
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Thus,  the  di -methyl  telluride  di-iodide  was  distilled  in  a 
sulphuric  acid  bath  at  100°C  evolving  methyl  iodide  leaving 
a  residue  of  elemental  tellurium  which  was  recovered. 


D.  Preparation  of  TeF^ 

As  stated,  TeF^  gas  was  chosen  as  the  form  into 
which  the  tellurium  samples  were  converted  for  mass  spectro- 
metrlc  analysis.  Fluorination  of  elemental  tellurium  can 
result  in  the  production  of  both  TeF^  and  TeF^.  According 
to  J.  H.  Junkins  et  al  (56),  TeF^  is  formed  at  a  temperature 
of  200°C  in  the  presence  of  excess  F2;  whereas,  TeF^  is 
formed  at  a  temperature  of  l80°C  with  excess  Te.  In  the 
method  used,  fluorine  was  present  in  a  large  excess  and  this 


resulted  in  a  quantitative  conversion  of  tellurium  to  TeF 
at  200°C . 
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The  fluorine  line  built  for  this  purpose  is  shown  in 
figure  5.  The  fluorine  source  was  a  commercial  cylinder  of 
fluorine  obtained  from  Matheson  Company  of  New  Jersey.  Many 
precautions  are  required  since  fluorine  is  most  reactive 
and  will  react  with  everything  except  inert  gases  and  already- 
fluorinated  compounds.  The  reaction  of  fluorine  with  some 
metals  such  as  copper  and  monel  is  retarded  after  an  initial 
fluoride  coating  is  formed  so  these  metals  can  be  safely  used 
for  construction  of  a  fluorine  line.  Therefore,  the  tellurium 
was  fluorinated  and  the  TeF^  trapped  in  a  monel  line. 


U  i  \  r  o  cj  ©  n 
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Dry  TeF^  is  fairly  non-reactive  so  that  after  being 
purified  several  times  to  remove  any  unreacted  fluorine,  the 
TeF^  could  be  collected  into  glass  break-seal  tubes.  All 
glass  stopcocks  were  greased  with  Kel-f#90  grease  since 
violent  reactions  could  occur  with  ordinary  greases  upon  con¬ 
tact  with  fluorine. 

Safety  precautions  included  locating  both  the  fluorine 
cylinder  and  the  entire  fluorine  line  under  a  fume  hood,  and 
the  wearing  of  a  rubber  coat,  polyethylene  gloves,  and  face 
mask  during  f luorinations . 

The  fluorination  procedure  is  as  follows.  The  tellurium 
sample  is  placed  in  a  monel  boat  and  placed  in  the  furnace 
which  is  heated  to  200°C.  With  closed,  V^,  Vg,  and 
open,  the  line  is  evacuated  back  as  far  as  V^.  After  pump¬ 
ing  on  the  sample  to  dry  it  thoroughly  and  testing  the  line 
for  vacuum  tightness,  T^  is  cooled  with  liquid  oxygen,  valves 
and  are  closed  and  V^,  V^,  and  are  opened  to  admit 
a  small  flow  of  nitrogen  through  the  line  controlled  by  use 
of  a  flowmeter.  Trap  TQ  is  now  cooled  with  liquid  oxygen. 

With  the  needle  valve  of  the  fluorine  cylinder  closed, 
the  main  valve  is  opened  and  closed,  trapping  a  small  volume 
of  fluorine  between  these  valves.  is  then  opened  and  the 

needle  valve  opened  allowing  the  fluorine  to  pass  through  the 
line.  The  fluorine -nitrogen  mixture  is  passed  through  a 
NaF  trap  to  remove  HF.  T1  freezes  out  any  remaining  HF  and 
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water  vapour  present ,  and  the  TeF^  formed  is  trapped  in  . 

The  flow  of  nitrogen  is  maintained  for  a  further  30  minutes 
to  ensure  scavanging  of  elemental  fluorine  from  the  system. 

After  scavenging,  and  are  closed  and  the  liquid 
oxygen  dewar  removed  from  T^.  This  vapourizes  the  frozen 
TeF^  and  releases  any  F^  which  might  have  been  trapped  by  the 
TeF^  molecules.  The  TeF^  is  again  frozen  down  in  T^  and 
is  opened  to  evacuate  the  line.  This  is  repeated  several 
times  to  ensure  that  no  F^  is  present  before  admitting  the 
TeF^  to  the  glass  line. 

With  and  closed,  is  opened  and  the  liquid 
oxygen  dewar  moved  from  T£  to  the  break-seal  tube.  After 
the  sample  is  frozen  down,  is  closed  and  the  liquid  oxy¬ 
gen  is  replaced  by  a  liquid-solid  alcohol  bath.  is  then 

opened  allowing  a  final  evacuation  of  the  break-seal  tube  to 
remove  any  traces  of  silicon  fluorides  formed.  is  then 

closed  and  the  break-seal  tube  removed. 

E .  Mass  Spectrometric  Analysis 

ioq  122 

In  order  to  determine  the  Te  J  /Te  ratio,  the 
TeF^  gas  was  analyzed  in  a  mass  spectrometer  using  a  12-inch 
radius  90°  magnetic  analyzer  to  separate  ion  currents  of  dif¬ 
ferent  masses.  The  TeF^  molecules  were  ionized  by  an  electron 
beam  travelling  normally  to  the  inlet  flow  direction.  The 
ions  formed  were  then  linearly  accelerated  by  a  3500-volt 


. 


81 


potential  difference  and  subsequently  separated  in  the 
magnetic  field.  The  ion  currents  are  introduced  in  turn 
into  the  collector  slit  by  magnetic  scanning.  Each  ion 
current  is  amplified  and  displayed  as  a  peak  on  a  continuous 
recording  chart  (see  Fig.  6).  Since  fluorine  has  only  one 
stable  isotope,  F1^,  the  spectra  yield  directly  the  abun¬ 
dances  of  the  tellurium  isotopes. 

In  single  collection  with  the  tracing  of  peaks,  a 
realizable  precision  is  0.1  percent.  This  limitation  results 
because  the  peaks  are  traced  at  different  times  and  so  fluc¬ 
tuations  in  the  mass  spectrometer  source  limit  the  precision. 
The  precision  can  be  improved  by  the  simultaneous  collection 
of  two  ion  beams  containing  the  isotopes  of  interest  and 
making  a  direct  ratio  measurement  of  the  two  ion  currents. 
Fluctuations  in  the  ion  production  effect  both  ion  currents  in 
the  same  way  and  so  the  ratio  remains  unaltered.  By  intro¬ 
ducing  a  standard  and  an  unknown  sample  alternately  into  the 
mass  spectrometer  by  a  magnetic  valve  system  (McKinney  et  al, 
Wanless  and  Thode  (75))  a  precision  of  0.01  percent  can  be 
obtained  in  comparing  isotope  ratios  of  different  samples. 

In  the  ionization  process,  TeF^  was  split  into  many 

singly  and  doubly  charged  ion  species  as  listed  in  Table  VII. 

The  TeFl  ion  species  are  the  most  abundant  and  so  ion  cur- 
5 

rents  of  masses  217  (Te^^F^)  and  225  (Te^  ^F^)  were  simul- 
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taneously  collected  with  two  slits  in  the  collector.  A 
narrow  slit  (0.024")  was  used  to  collect  the  mass  217 
ion  current  while  a  wide  slit  (0.040")  was  used  for  the 
mass  225  ion  current.  With  the  former,  the  resolution 
of  the  instrument  was  1/350  while  with  the  wider  slit, 
the  resolution  was  1/200.  The  lower  resolution  is  per- 
missable  in  the  latter  case  since  the  closest  ion  current 
is  of  mass  223  (Te  F^)  .  The  spectrum  in  Figure  6  was 

traced  using  the  narrow  collector  slit. 

170  122 

The  Te  /Te  ratio  is  obtained  directly  by  means  of 
a  five-figure  integrating  digital  voltmeter-ratiometer  and 
a  digital  printer  in  the  measuring  circuit  of  the  mass 
spectrometer.  The  Te^^F^  and  Te1^°F^  ion  beams  are  col¬ 
lected  simultaneously  and  are  amplified  by  vibrating  reed 
electrometers.  The  currents  are  converted  to  voltages  by 
means  of  lO1^  and  10^  ohm  grid  leak  resistors  for  the  122 
and  130  mass  currents  respectively.  The  integrating  digital 
voltmeter  integrates  an  input  voltage  over  a  selected  time 
and  converts  it  to  a  proportional  frequency.  In  passing 
through  a  counter  section,  this  frequency  is  referred  to  a 
time  base  frequency  so  that  the  number  printed  is  really  a 
ratio  of  the  converted  frequency  to  the  time  base  frequency. 
If  an  additional  voltage-to-f requency  converter  is  used,  a 
second  voltage  can  be  converted  Into  a  frequency  and  used 
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in  place  of  the  time  base  frequency  in  the  IDVM.  The  IDVM 
then  is  a  ratiometer  and  displays  the  ratio  of  the  two  volt¬ 
ages  . 

180  122 

The  Te  J  /Te  ratios  of  a  standard  sample  and  an  un¬ 
known  sample  can  be  calculated  directly  from  the  printed  output 
providing  the  major  ion  currents  of  the  two  samples  are  matched. 
Because  of  memory  effects,  it  was  necessary  to  allow  at  least 
five  minutes  after  switching  of  the  magnetic  valves  before 
recording  the  ratio  in  order  to  get  reproduceable  results. 
Approximately  six  sets  each  of  ten  printed  ratios  for  each 
sample  were  obtained. 


TABLE  VTI 


Ion 

Abundance 

TeFt 

5 

58.3 

TeF^ 

2.8 

TeF* 

5.5 

TeF* 

2.8 

TeF+ 

0.4 

Te+ 

6.5 

H — b 

TeFj^ 

4.2 

TeF++ 

7-9 

TeF^ 

4.6 

TeF++ 

3.6 

m 

Te 

3.3 

RESULTS  AND  DISCUSSION 


The  results  of  the  theoretical  calculations  as 
summarized  in  Tables  IV  and  VI  indicate  that  the  Te13°/Te122 
ratio  may  be  altered  by  over  3^  in  hypothetical  reactions. 

The  experimental  observations  are  consistent  with 
some  of  the  theoretical  predictions.  The  measured  Te1^°/Te122 
ratios  for  the  samples  analyzed  were  referred  to  the  isotopic 
composition  of  a  standard  sample  of  telluride  from  Copper 
Cliff,  Ontario.  The  isotopic  deviation  from  the  standard 


is  expressed  in  6^0  units  where 


>130 


(Te^O/Tel22) 
v  e  /ie  '  sample  -  1 


(Te^/Te122)  H  , 

.  '  'standard 


x  1000 


Generally  speaking,  a  standard  deviation  of  0.02  percent 
could  be  realized  in  the  measurements.  The  results  are 
summarized  in  Table VIII  and  depicted  graphically  In  Fig.  7. 

A.  Chemical  reduction  of  potassium  tellurite 

Analysis  showed  that  the  tellurium  precipitated  in 
a  10  percent  reduction  of  tellurite  was  0.71  percent  depleted 
in  the  heavier  isotope  in  comparison  with  the  initial  tel¬ 
lurite.  Using  equation  (19)  developed  in  the  theory,  this 
value  can  be  extrapolated  with  the  assumption  of  first  order 
kinetics  to  give  the  ratio  of  the  rate  constants  ^122/^130 
for  the  competitive  reactions 


Table  VIII  Te13°/Te122  FRACTIONATION  STUDY 
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122  = 

Te  0 

kl  22 

— »  Te122 

3 

Te^^0= 

kl30 

^  leo 

.  Te130 

le  u3 

^  ieo 

This  extrapolation  gives  k122//k130  =  1*00^1  "  0*0°°3  meaning 
that  the  lighter  isotopic  ion  ruptures  0.7  percent  faster. 

This  is  to  be  compared  with  the  theoretically  predicted  re¬ 
sults  of  Table  VI.  Looking  at  Table  VI.,  "most  reasonable" 
values  of  1.00568  and  1.00718  for  mass  and  atom  fragments 
respectively  are  predicted.  TeO^  has  a  pyramidal  XY^  structure 
as  shown  in  Figure  12,  Appendix  II.  The  bond  being  broken 
is  equally  as  strong  as  the  other  bonds  so  that  by  Wolfsberg's 
reasoning  (136),  the  reduced  mass  term  will  be  closer  to  the 
reduced  mass  of  the  Te  -  0  diatomic  bond.  Thus  the  value  of 
1.0071  ±  0.0003  compares  very  closely  to  the  predicted  re¬ 
sult  of  1.00718  for  the  atom- fragment  model.  This  is  also 
consistent  with  results  found  for  selenite  reduction  (89) . 

B.  Microbiological  reduction  study 

For  measurement,  the  elemental  tellurium  fraction 
precipitated  between  4l  and  53  percent  reduction  was  used. 

This  value  can  also  be  extrapolated  to  give  the  ratio  k:]_22//k130 

The  extrapolation  is  however  more  difficult  and  is  shown 
graphically  in  Figure  8.  The  quantity  "r"  in  equation  (19  ) 
is  for  the  integrated  isotopic  composition  of  the  product. 


89 


r~i  t  <  r— ■ 

i  iwUi  » 


v/ 1  r>?  r  "tAM  n~  'CnTHDir  p  r>  *  * OAC!~tr> V 

V  s-**!  tilt  i  i  i  VI  '*  U i  i  w'  j  V  J  1  V  u  vj  i V J  i  'w.  »  I  8  v1 1\ 


CF  PRO 


rvwOUC  s 


WITH 


DcprcM' 

k  w.  i  \  v  w  ♦  « 


REACTION 


1/ 

822 


ISOTOPIC  COMPOSITION  OF  FRACTIONS 


90 


Three  plots  of  "r"  as  a  function  of  percentage  reaction  are 
shown  in  Figure  8  as  dashed  lines  for  values  of  k122/kl30  = 
1.005*  1.006,  and  1.007.  The  point  X  indicates  that  if  a 
reaction  whose  k]_22/k]_3Q  =  1*007  has  gone  40  percent  to  com¬ 
pletion,  the  Te122/Te1'^0  ratio  of  that  total  40  percent  is 
1.0053  times  that  of  the  initial  reactant. 

For  the  microbiological  study,  however,  a  plot  of 
isotopic  composition  of  fractions  vs.  percentage  reaction 
is  necessary.  Five  plots  for  ki22Aq3Q  =  1*0°50j  1.0055, 
1.0060,  I.OO65*  and  1.0070  are  shown  by  solid  lines  in  Fig.  8. 
The  information  given  by  such  a  plot  is  represented  by 
construction  Y.  If  k]_22/ki3Q  =  l-00?*  and  It  is  desired  to 
know  the  average  isotopic  composition  for  the  fraction 
between  10  percent  and  20  percent  reduction,  the  horizontal 
line  is  drawn  at  Y  so  that  the  shaded  areas  above  and  below 
the  line  are  equal.  In  this  example,  the  Te  /Te  D  ratio 
of  the  10  to  20  percent  fraction  of  the  product  is  1 . OO58 
times  that  of  the  initial  reactant. 

The  above  construction  is  performed  in  reverse  for 
the  experimental  result.  At  Z,  the  horizontal  line  is  the 
mean  value  of 


('T>r:  ^reduced  fraction 

(Tel22/Te1^°) 


=  1.0021 


initial  reactant 
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for  the  fraction  between  4l  and  53  percent  reaction.  It  is 
seen  that  the  line  for  ^122/^130  =  1*0060  very  closely  equal¬ 
izes  the  shaded  areas.  The  measurement  error  of  10.02  percent 
is  represented  by  the  vertical  line  at  Z  and  its  extremities 
touch  the  curves  corresponding  to  ^22/^130  =  1.0055  and 
I.OO65.  Therefore,  for  the  bacterial  reduction,  ^i22///^<:130  = 
1.0060  t  0.0005.  This  compares  very  favourably  with  the 
result  of  the  chemical  reduction  study  and  suggests  that  in 
the  bacterial  reduction,  the  effect  noted  arises  essentially 
in  tellurium-oxygen  bond  breakage. 

^122/^130  can  ^etermined  more  accurately  by 
measurements  made  of  fractions  produced  near  the  beginning 
of  the  reaction.  Such  samples  are  available,  but  measure¬ 
ments  are  being  withheld  until  the  mass  spectrometer  is 
adapted  to  handle  smaller  samples. 

C.  Natural  .isotopic  abundance  studies 

Two  of  the  natural  samples  and  the  commercial  tel¬ 
lurite  have  the  same  Te122/Te1^0  ratio  within  the  measurement 
errors.  Similarly  in  the  survey  of  natural  selenium  isotopes 
(60)  a  large  number  of  samples  had  a  similar  Se^/Se^2  ratio. 
In  both  cases,  the  constancy  of  isotopic  composition  can  be 
related  to  high  temperature  processes  where  the  isotope  frac¬ 
tionation  is  very  small.  (The  isotope  fractionation  generally 
decreases  with  temperature  as  exemplified  in  the  calculations 
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of  Tables  IV  and  VI.)  In  the  selenium  study  (60),  large 
alterations  to  the  Se^/Se^  ratio  occurred  only  in  soil  and 
plant  samples  where  the  temperature  of  the  processes  would 
be  very  low.  In  the  tellurium  study,  however,  the  sample  of 
tetradymite  from  the  west  shore  of  the  Okanagan  Lake  was 
0.4  percent  depleted  in  Te  J  .  This  result  is  surprising. 

A  measurement  of  yield  rules  out  any  possibility  that  the 
sample  was  isotopically  altered  in  the  chemical  extraction. 
This  sample  behaved  the  same  as  the  others  in  hexafluoride 
preparation  and  in  mass  spectrometric  analysis.  It  would 
seem  that  the  result  is  real  but  it  should  be  interpreted 
as  one  measurement  of  one  sample.  No  attempt  can  be  made 
to  explain  this  result  until  more  natural  samples  are  ex¬ 
amined. 

In  conclusion,  therefore,  this  preliminary  survey 
shows  that  the  Te^^/Te1^0  ratio  can  be  altered  chemically 
and  microbiologically  in  agreement  with  theoretical  pre¬ 
dictions.  Further,  it  would  seem  that  a  more  extensive  study 
of  natural  samples  might  be  beneficial  in  elucidating  natural 


processes . 
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APPENDIX  I 

CALCULATIONS  OF  ISOTOPIC  VIBRATIONAL  FREQUENCY  SHIFTS 
FOR  TELLURIUM  -  CONTAINING  MOLECULES 

Values  for  the  normal  vibrational  frequencies  of 
many  tellurium-containing  molecules  are  available  from 
infra-red  spectra  and  can  be  found  in  various  publications. 
However,  since  no  consideration  is  given  for  differing 
isotopic  masses  therein,  it  is  assumed  that  these  experi¬ 
mental  values  apply  to  the  most  abundant  tellurium  isotope, 

ISO  122 

Te  .  The  frequencies  of  the  Te  containing  molecules 

were  then  calculated  from  these  known  frequencies  using 

normal  force  equations  which  describe  molecules  of  different 

symmetries . 

A .  Diatomic  Molecules 

In  considering  diatomic  molecules,  the  anharmonic 
oscillator  is  chosen  as  a  model.  The  potential  energy  function 
is  then  given  in  the  form 

U  =  f(r  -  re)2  -g(r-re)3+  .  .  .  (i.l) 

where  r^  is  the  equilibrium  internuclear  distance  and 
(r  -  r  )  is  the  displacement  from  this  equilibrium  position. 


A2 


Schroedinger ’ s  equation  for  such  a  model  then  becomes 

2 

+  2EJL  (E  .  u)^  =  0  (1.2) 

h 


If  the  anharmonicity  if  small,  (g«f),  the  energy  values 
of  this  equation  will  be  the  eigenvalues  given  by 


E 


vib 


=  he 


cdq  (v+i) 


-  0}  x  (v+i) 
e  e  d ' 


+  (v+4)3  •  •  • 


where  v  is  the  vibrational  quantum  number,  is  the 
classical  vibrational  frequency  in  cm-1,  and  q  is  the 


reduced  mass  of  the  two  atoms,  ( 


m1m2 

m^+m^ 


). 


(1.3) 


Classically,  for  a  harmonic  oscillator 


v 

osc 


-1 

sec 


where  v  is  the  vibrational  frequency  and  the  restoring 

osc 

force  is  -K(r  -  r  )  during  a  displacement  of  (r  -  r  ). 

The  potential  energy  curves  for  isotopic  molecules 
are  identical  to  a  high  degree  since  they  depend  only  on 
the  motions  of  the  electrons  and  the  coulomb  repulsion  of 
the  nuclei  (48).  This  means  that  the  force  constant  K  can 


be  assumed  independent  of  isotopic  substitution.  There¬ 
fore  the  ratio  of  the  classical  frequencies  between  two 
isotopic  molecules  would  be  given  as 
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Or  in  terms  of  the  most  abundant  isotope. 


v . 


isotope 
v 


M- 


M- 


isotope 


From  the  relationship  that  v  =  cu>, 


=  P 


P 


(1.4) 


Thus  from  equations  1.3  and  1.4  the  energy  levels  of  the 
isotopic  molecule  are  given  by 


(Ev)x  =  he  [  ptt>e(v+i)  -  p2“exe(v+i)2  +  .  .  ]  (1.5) 


The  only  vibrational  energy  transitions  possible 
are  those  dictated  by  the  selection  rule,  i.e.,  Av  =  ±1 . 
Thus  the  vibrational  energy  change  in  a  transition  from 
v  =  0  to  v  =  1  is  given  by  equation  1.3  as 


A4 

AEv  =  he  [05g  -  20)exe  +  .  .  .  ]  (1.6) 

=  hc<^> 

Prom  1.4,  the  energy  change  for  av=ltov=0  transi¬ 
tion  for  the  isotopic  molecule  is  given  in  terms  of  the 
transition  for  the  abundant  molecule  as 

(AEv)1  =  ho  [  p<»e  -  2p2coexe  +  .  .  .  ]  (1.7) 

=  hc<u. 
i 


Prom  equation  (l.7)j  the  expression  for  the 
vibrational  frequency  of  the  isotopic  molecule  in  terms  of 
the  frequency  of  the  more  abundant  molecule  as 


O), 
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2p 


co  x  + 
e  e 


(1.8) 
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Exoerimental  values  for  and  cu  x  for  the  Te 

e  e  e 

containing  molecules  are  obtained  from  observed  data 
(48),  (67)*  Then  from  1.4, 
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and 


Table  IX  gives  the  calculations  for  the  vibra- 

122 

tional  frequencies  of  Te  -containing  molecules  from  the 
available  data.  These  were  computed  with  the  IBM  7040 
computer  program  found  on  page  A37,Appendix  II. 


B.  Polyatomic  Molecules 

In  these  molecules,  a  simple  model  cannot  be  chosen 

for  the  vibrational  frequencies.  Detailed  information  must 

be  obtained  about  the  structure  of  the  molecule  in  order 

to  describe  the  force  fields  and  set  up  equations  for  the 

vibrational  frequencies.  The  force  equations  which  seemed 

to  best  describe  the  observed  spectra  were  chosen.  The 

force  constants  were  solved  in  terms  of  the  most  abundant 
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Te  isotope  and  experimentally  observed  frequencies. 

Then  the  other  isotopic  masses  were  substituted  to  find  the 
corresponding  isotopic  frequencies. 

In  several  cases  even  the  best  description  of  the 
molecule  did  not  give  results  in  perfect  agreement  with  the 
observed  vibrational  frequencies.  Therefore,  linear  correc¬ 
tion  factors  are  used  when  needed  to  align  the  calculated 
frequencies  to  those  observed. 


A  6 


A7 


TePg 

Tellurium  hexafluoride  Is  of  octahedral  structure 
of  0^  symmetry.  The  octahedral  molecules,  XYg,  of  symmetry 
0^  each  show  six  fundamental  vibration  frequences.  Of  these 
six  modes,  03^  Is  non-degenerate,  03^  Is  double  degenerate, 
and  03^,  cd^,  ao^,  and  03^  are  triply  degenerate.  The  form  of 
these  are  shown  in  Fig,  9  in  which  the  vibrational  frequen¬ 
cies  are  labelled  03^.  .  .  03^.  cd^,  oo^,  and  are  asso¬ 
ciated  with  essentially  valency  vibrations  and  03^,  and 
03g  with  essentially  deformation  vibrations,  od^,  03^,  and  03^ 
appear  in  the  Raman  spectrum,  03^  an(i  ^  appear  in  the  infra¬ 
red  spectrum,  03^  does  not  appear  in  either  because  the  trans- 
tion  is  forbidden  for  both  Raman  and  infra-red  spectrum. 

This  means  that  only  the  first  five  fundamental  vibrational 
frequencies  can  be  directly  determined  spectroscopically. 

03g  has  to  be  obtained  either  from  combination  tones  or  from 
specific-heat  data  and  it  is,  therefore,  less  accurate.  With 
TeFg,  03^  was  not  directly  observed  but  was  also  derived  from 
overtones.  Heath  and  Linnett  (46)  feel  that  any  value  be¬ 
tween  340  and  370  cm-1  would  account  satisfactorily  for  the 
overtones  observed  and  that  the  value  cannot  be  fixed  more 
closely.  However,  more  recent  publications  (40)  and  (66) 
assign  a  value  of  327  cm 


for  this  vibration. 
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Heath  and  Linnett  (46)  have  examined  the  applica¬ 
tion  of  two  different  force  fields  to  octahedral  XY^  molecules, 
and  in  particular  TeF^.  They  found  the  "Orbital  Valency 
Force  Field”  equations,  which  are  based  on  the  view  that 
valency  bonds  are  related  to  the  form  of  the  atomic  orbitals, 
to  give  superior  results  to  the  "Simple  Valency  Force  Field" 
equations.  The  following  table  shows  a  comparative  agree¬ 
ment  between  observed  and  calculated  frequencies  using  the 
OVFF  and  SVFF  approach  respectively. 


Vibrational  Frequencies  of  TeF^  in  cm 


TeF^ 

(JO 

1 

(0 

2 

or 

d 

“4 

“5 

“6 

Observed  (46  ) 

701 

67U 

751 

327 

313 

165 

Calculated  OVFF 

698 

66  9 

775 

351 

313 

388 

Calculated  SVFF 

704 

674 

780 

277 

320 

224 

It  is  noted  that  the  normal  vibration,  is  in 
fairly  good  agreement  with  the  calculated  value;  whereas 
the  value  is  not  satisfactorily  considered  in  the 
OVFF  approach. 


. 

■ 
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Heath  and  Linnett  (46)  give  the  following  equations 
for  TeF^  from  the  OVFF  approach : 


\  =  (kx  +  8a) /nip  (1.9) 

=  (kp  -  f  |  +  3)  rap  (I-IO) 

A3  +  A4=  fll  (•*'/raF  +  2/,mTe^  +  f22^/mTe  +  1'/mF^ 

-  4  2(f12/mTe)  (1. 11) 

A3Ait  =  (fllf22  ‘  f122)  (1/mF2  +  6/raFraTe)  C1*12) 

where  ,, 

fn  =  <ki  -E§+  4A) 

f12  =  (2  2  A  -  2  B/R) 

f22  =  ^kB  +  ^R  +  2A^ 

A5  =  (kB  +  +  4A)/mp  (1.13) 

A6  =  (kB  +  2A^/raF 


(1.14) 


All 


where, 


2  2  2 

=  4tt  c  ^  =  frequency  in  cm 

mp  =  mass  of  fluorine  atom  (amu) 
mTe  =  mass  of  tellurium  atom  (amu) 

=  the  bond  stretching  constant  of  the 
Te-F  bond 

kg  =  the  angle  bending  constant  of  the 
F-Te-F  bond 


-B  and  2A  are  the  first  and  second  differentials  of  the 
potential  energy  function  representing  repulsion  between 
non-bonded  atoms. 


The  constants  have  been  calculated  by  Heath  and 
Linnett  to  be : 

k^  =  4. 85  x  10^  dyne/cm 
k^.  =  0,78  x  10^  dyne/cm 

XD 

A  =  0.075  x  10^  dyne/cm 
~  =  0.0115  x  10^  dyne/cm 

The  values  used  for  the  fluorine  and  tellurium  masses  are 
those  given  by  Semat  (94). 
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It  is  immediately  seen  from  the  equations  that 
A^,  A,-,  and  A^  are  independent  of  the  mass  of  tellurium 

and*  therefore,  remain  the  same  for  isotopic  substitution 
of  tellurium.  For  this  reason,  the  poor  agreement  of  the 
calculated  value  of  ^  with  the  observed  value  can  be 
overlooked  in  calculating  isotopic  shifts.  Using  equations 
I. 11  and  1.12,  A^  +  A^  and  A^A^  are  obtained.  Then, 

A^  and  A^  and  subsequently  co  and  ^  can  be  obtained  by 
using  the  equation: 

A2  -  (A^  +  A^)  A  +  A^A^  =  0  (1.15) 

whose  roots  are  given  by 

2A  =  ?*3  +  ±/(A 3  +  A^)2  -  )  (I.l6) 

An  IBM  7040  computer  program  was  set  up  to  cal¬ 
culate  and  ^  for  a11  isotopic  masses  of  tellurium  using 
the  foregoing  considerations.  (See  pagps  Ab3,  65;  Appendix  II ) . 

The  values  of  an(^  ^4  were  found  to  be  the 


following : 
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Isotope 

“3 

-1 

cm 

“4 

-1 

cm 

m  120^ 

Te  F^ 

775.^2 

356.47 

Tel22F£ 

773.78 

355.79 

Tel23F£ 

772.98 

355.46 

Tel2^F^ 

772.19 

355.13 

Te125F^ 

771.42 

354.80 

Te12^F^ 

770.66 

354.48 

Te128F^ 

769.17 

353-85 

Te130F£ 

1  767.72 

353.22 

130 

The  computed  values  for  Te  J  are  now  compared 
with  the  observed  values  as  follows : 


Tel3°F£ 

“3 

-1 

cm 

“4 

-1 

cm 

Observed 

752 

327 

Computed 

767.72 

353.22 

Correction  factors 


752 


and 


327 


are  used  to 


767-72  353.22 

align  the  calculated  vibrational  frequencies  for  the  isotopic 

130 

mass  with  the  observed  vibrational  frequencies  for  Te  F^. 


« 

have  : 


Al4 


Then  in  terms  of  the  corrected  frequencies,  we 


0) 

1 

cm"  ^ 

CD 

3 

cm”'*' 

°3 

„  -1 
cm 

CD 

4 

-1 

cm 

“5 

-1 

cm 

“6 

-1 

cm 

m  120- 

Te  F^ 

701 

674 

758.52 

330.00 

313 

164 

Te122F^ 

701 

Gib 

756.92 

329.38 

313 

l64 

Te123F^ 

701 

Gib 

756.14 

329.07 

313 

l64 

m  124,, 

Te  F^ 

701 

Gib 

755.37 

328.77 

313 

164 

Te125F^ 

701 

Gib 

754.61 

328.47 

313 

164 

Te12^F^ 

701 

Gib 

753-87 

328.17 

313 

164 

Te128F^ 

701 

Gib 

752.41 

327.58 

313 

164 

Te13°F£ 

701 

Gib 

751.00 

327.00 

313 

i64 

TeBr^ 

TeBr^  belongs  to  the  group  XY^  having  a  tetra¬ 
hedral  structure.  There  are  four  normal  modes  of  vibration; 
one  totally  symmetric  (non-degenerate)  mode  (^)  of  symmetry 
type  one  doubly  degenerate  mode  (uy>)  of  type  E,  and 

two  triply  degenerate  modes  and  ^ )  of  type  F^. 

These  normal  vibrations  are  shown  in  Fig. 10. 
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Various  workers  have  set  down  force  equations  for 
this  molecule.  Venkateswarlu,  Somasundaram,  and  Krishna 
(l2l)  give  a  method  of  evaluation,  but  to  quote  C.W.F.T. 
Pistorius,  (85) *  "  .  .  .they  seem  to  have  made  a  numeri¬ 

cal  error,  because  they  list  real  values  while  the  correct 
results  according  to  their  method  undoubtedly  are  imaginary. 
I  found  the  same  to  be  true. 


The  best  results  were  found  to  be  those  of  the 
OVFF  approach  of  Heath  and  Linnett  (45).  They  give  the 
force  equations  for  some  XY^  tetrahalides  as: 

\  =  (k  +  8a)/mf 

Ag  =  (3k/(r0)2  +  2A  +  B/R° ) Mp 
A  +  =  (k  +  8A/3  -  4B/3R°)(1/MF  +  4/3MTe) 

+  (2k/(r0)2  +  4A/3  +  10B/3R0) (1/Mp  +  8/3MTe) 
-  (8A/3  -  4B/3R°)(8/3MTe) 

=  (R  +  8A/3  -  4B/3R0) (2k/(r0)2  + 

4A/3  +  10B/3R0)  -  i/2(8a/3  -  4B/3R0)2 
l/Mp2  +  4/MTe.MF 

Miss  R.  Thanalakshmi  (103)  calculated  the  force 
constants  of  the  OVFF  equations  of  Heath  and  Linnett  (45) 
and  gives  the  following  values : 


■ 

'A 
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k  =  1.4500  x  10^  dyne/cm. 
k/(r0)2  =  0.0192  x  105  dyne/cm. 

B/R°  =  0.0040  x  10^  dyne/cm. 

A  =  0.0300  x  108  dyne/cm. 

It  is  noted  that  >  +  A^  are  independent  of  the 

mass  of  tellurium  and  therefore  remain  the  same  for  iso¬ 
topic  substitution  of  tellurium.  The  method  of  determining 
the  roots  and  thus  the  frequencies  ^  and  ^  is  exactly 
the  same  as  was  done  for  TeF^.  See  equations  1.15  and  I.l6. 

Again  an  IBM  7040  computer  program  was  set  up 
(Pages  A63  and  6 5)  and  results  are  given  as  follows: 


Isotope 

“3 

-1 

cm 

“4 

-1 

cm 

Te120Brij 

244.92 

61.49 

Tel22Brz+ 

243.96 

6l .  36 

Te123Br/| 

243.50 

61.30 

Te12l|BrZ| 

243.04 

61.23 

Te125BrZ| 

242.58 

61.17 

Te12^Br^ 

242.14 

6l.ll 

Te128Br4 

241.26 

60.99 

Te13°BrZ| 

240.41 

60 . 86 
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Again  a  discrepancy  is  found  between  the  computed 
and  observed  frequencies  for  Te18°Br^ 


Te13°Br2| 

-1 

cm 

“4 

-1 

cm 

observed 

209 

64 

calculated 

240.91 

6o .  86 

Upon  correction,  the  final  normal  frequencies 

are  then: 


CO 

1 

,  -i 
cm 

“2 

-1 

cm  . 

“3 

cm  1 

“4 

-1 

cm 

Tel20BrZ| 

190 

51 

212.92 

64.66 

1  P  ° 

Te  Br^ 

190 

51 

212.09 

64.52 

Te123BrZ| 

190 

51 

211.68 

64.46 

Te12Z+BrZj 

190 

51 

211.29 

64.39 

Te125BrZ| 

190 

51 

210.89 

64.32 

Te12^Br;4 

190 

51 

210.50 

64.26 

Te128BrZ| 

190 

51 

209.74 

64.13 

Te18°BrZ| 

190 

51 

209.00 

64.00 

SCTe 

This  molecule  is  of  the  triatomic  linear  type 
having  a  symmetric  stretching  mode  w  ,  a  bending  mode 
and  an  asymmetric  stretching  mode  ^ •  Of  these,  the 


asymmetric  stretching  mode  (u>2)  is  doubly  degenerate,  and 
^  and  non -degenerate  (49)  .  The  forms  of  the  normal 
vibrations  would  look  like  the  following: 


X  Y  Z 


Wentink  (120)  gives  the  following  force  equations 
for  this  molecule: 


A1  +  >3  =  K^l/Mg  +  1/MC)  +  K2(1/Mc  +  l/MTe)  -  2K3/Mc 

+  Mr 

12 


AiA3 


° _ 5!®.  (k  k  -  k 


S  C  Te 


A2  -  (l/r1  r2  )  •  (r1  /MTe  +  r2  /Mg  +  r^2  /Mc^Ka 


where 

K1  =  stretching  constant  for  the  S-C  bond 

K0  =  stretching  constant  for  the  C-Te  bond 

Kg  =  stretching  interaction  constant 

K  =  bending  constant 
a 

r1  =  S-C  internuclear  distance 
rp  =  C-Te  internuclear  distance 
r S-Te  internuclear  distance 
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He  gives  the  values  for  these  constants  as  follows: 

=  7*04  x  10^  dyne/cm. 

K^>  =  4. 58  x  10^  dyne/cm. 

=  0.30  x  105  dyne/cm. 

Ka  =  0.52  x  lO^  dyne/cm-rad 
rx  =  1.557  A 
r2  =  1.904  A 
r12  =  3.461  A 

As  noted,  each  A  is  dependent  on  the  mass  of  tellurium 
and  therefore,  each  is  affected  by  isotopic  substitution. 
Knowing  A^  +  A^  and  A^A^  one  uses  equations  1.15  and 
I.l6  as  before  for  obtaining  A^  and  A^  and  subsequently 
and  03  A^  is  calculated  directly  and  from  this 
is  obtained. 

Using  the  computer  program  (pagss  a63.>  64  )  the 
values  of  0^,  and  03  were  found  to  be  as  follows: 
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CL) 

1 

cm  2 

“2 
cm  1 

“3 
cm  2 

SCTe120 

1347.64 

340.05 

426.35 

SCTe122 

1347.54 

340.00 

425.45 

SCTe123 

1347.48 

339.98 

425.00 

SCTe12*4 

1347.43 

339-95 

424.57 

SCTe125 

1347.38 

339.93 

424.14 

SCTe228 

1347.33 

339.91 

423.72 

SCTe128 

1347.23 

339-86 

422.89 

SCTe180 

1347.14 

337.00 

422.08 

Again  a  discrepancy  is  found  between  the  calculated  and 
observed  frequencies  for  SCTe232. 


SCTe130 

0) 

1 

-1 

cm 

0) 

2 

cm  2 

0) 

3 

-1 

cm  ..  .. 

Observed 

1347.00 

337-00 

423.00 

Calculated 

1347.14 

337.00 

422.08 

When  this  has  been  corrected  for.,  the  final  normal 
frequencies  are  then: 
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i — i 

i 

i  £ 

3  o 

-1 

cm 

0) 

3-1 

cm 

SCTe120 

1347.51 

340.50 

427.28 

SCTe12Z| 

1347.40 

340.00 

426.37 

SCTe123 

1347.35 

339.98 

425.93 

SCTe12^ 

1347.30 

339.95 

425.49 

SCTe125 

1347.24 

339.93 

425.06 

SCTe128 

1347.19 

339.91 

424.64 

SCTe128 

1347.10 

339.86 

423.81 

SCTe180 

1347.00 

337.00 

423.00 

H2Te 

This  molecule  is  of  non-linear  XYQ  structure 
having  the  form  of  normal  vibrations  shown  below.  Each 
vibration  is  non-degenerate. 


X 

i  1 


/ 


\  J 
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Due  to  the  chemical  instability  of  this  gas  (93) 
high  resolution  spectroscopy  has  been  limited  and  the 
vibrational  frequencies  given  are  peak  values  for  bands 
of  unresolved  spectral  lines.  Some  overlapping  of  bands 
is  also  reported. 


Triatomic  non-linear  molecules  have  been  studied 
by  Glockler  and  Tung  (4l)  with  the  assumption  of  general 
force  fields  which  produce  three  equations  with  four 
force  constants  for  isosceles  non-linear  triatomic  molecules 
as  follows : 


A1  +  A2  -  L11  +  L12  +  L33 


=  L33  (kjx  +  ~  2L^L 


13  31 


(1.17) 


^3  ~  L11  L12 


where 


-L 


11 


c  C~  c 

—  +  -ft  Sln  “  -  -M 


4 


cos  a 


1  c. 


C  c3 

-H  008  a  M 


sin  a  + 


1 


-l13  =  -M  sin  “  -  2C3 


-L3i  =  -TJ  sin  a  -  C3 


1 

P 


cosa 

M 


+ 


M 


1  ,  cosa 
(I  +  ~M~ 


sin  a 


L33  -  2C; 


1_  cosa 
p  "  M 


!C 


M 


3 


sin  a 
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where 

M  =  mass  of  central  atom 
m  =  mass  of  terminal  atoms 


1 

P 


1 

m 


a  is  the  apex  angle  as  shown 


The  observed  angle  a  for  the  H-Te  bond  is 
close  to  90°  (lOl)  (Rossmann  and  Straley  in  1956  gave  it 
as  89O30’  (93))*  The  two  observed  symmetric  frequencies 
are  2047  cm"1  and  2133  cm"1  and  the  anti -symmetric 
frequency  observed  is  860.79  cm  1  (50,  102),  Since  a  is 
so  close  to  90°,  the  equations  are  simplified  to 


+ 


P 


2C3  C4 
H  M 


P 
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giving 


A 


4C. 


1  +  *2  -  jl  t2C2  +  Cij  +  Cl  ]  -  ^ 


~  t2Cp(C1  +  C„)  -  40-^ ]  [-^p  -  ^p]  . 

0  1 1  M 


2V"1  'r  "4 
1 


(I 


A3  =  (Cl  -  C4>  £> 


It  is  seen  from  I.l8  that 
cies  for  <b  are  related  by 


the  isotopic  frequen- 


and  are  thus  independent  of  the  force  constants. 
Also* 


»  i 

A1A2 


A.  A 


1''2 


which  is  also  independent  of  the  force  constants.  Thus 
only  the  calculation  of  A^  +  AQ  depends  on  the  values  of 
the  force  constants  and  being  of  much  smaller  magnitude 
than  the  product  A^A,-,,  any  error  in  the  values  chosen  for 
the  force  constants  is  very  nearly  negligible  as  affecting 
isotopic  vibrational  frequency  computations. 


.18) 
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It  Is  seen  by  the  set  of  equations  (I.l8) 
that  only  three  equations  are  available  for  determining 
four  unknown  force  constants.  Thus  it  is  impossible  to 
determine  the  force  constants  uniquely.  However,  if 
is  assumed  known  and  the  other  constants  calculated  in 
terms  of  C^,  there  exists  only  a  limited  range  of  values  of 
which  give  a  ’ non- imaginary 1  solution.  The  relationship 
between  the  four  constants  is  plotted  in  Pig.  11,  Thus  the 
constants  have  a  limited  range  of  possible  values,  namely, 

=  0.6l  +  0.50,  C3  =  0  +  0.50,  C2  =  0.72  +  0.50,  and 
C1  =  2.05  +  0,50  (x  10^  dyne/cm.).  Therefore,  if  a  unique 
set  of  force  constants  is  to  be  found,  one  of  the  constants 
must  be  determined  by  some  other  means,  o>3  is  the  stretch¬ 
ing  frequency  of  the  hydrogen-tellurium  bond.  In  the 
polyatomic  molecule  its  value  is  affected  by  interaction 
with  the  other  forces,  but  one  can  get  an  approximate  value 
for  the  H-Te  stretching  constant  C±  by  considering  a 
diatomic  H-Te  molecule  where  the  interaction  forces  are 
absent.  Lippincott  and  Dayhoff  (67)  give  the  value  of 
for  TeH  as  2250  cm"1.  Using  the  relation  that 


where 


c  =  velocity  of  light 
\i  =  reduced  mass 
k  =  bond  constant 
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k  is  deduced  to  be  2.66  x  10^  dyne/cm.  for  the  diatomic 
H-Te  bond.  This  would  thus  delimit  the  values  of  the 
constants  to  those  values  appearing  to  the  left  in  Pig.  10 
where  has  the  maximum  value  of  2.55  x  10^  dyne/cm. 

Smith  and  Linnett  (99)  have  calculated  the 
constants  for  the  other  members  of  the  group  VIA  elements 
using  similar  considerations  along  with  the  relation: 

Clog  C-^  =  -6  log  r  +  a] 

from  excited  state  spectroscopic  data  to  deduce  and 
therefore  deduce  C^3  and  .  They  give  the  following 
set  of  force  constants:  (xlO^  dynes/cm) 


ci 

C2 

C3 

cit 

H20 

8 . 36 

0.75 

0.20 

-0.10 

H2S 

4. 4 ±0.2 

0.55±0.05 

0.2 ±0.2 

-0.15±0.5 

H2Se 

3 . 4  ±0 . 1 

0.38±0.06 

0.15*0.03 

-0. 15±0. 15 
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It  is  noted  that  the  value  of  the  constants 
decreases  with  increasing  order  in  the  periodic  table. 
Thus,  this  further  delimits  the  values  of  the  constants  to 
<  0.18  which  in  turn  limits  C2  <  0.25,  C1  >  2.52, 

>  -0.14.  It  is  noted  also  that  is  positive.  In 
terms  of  this  value  we  can  then  quite  accurately  assign 
the  following  values  to  the  H2Te  force  constants. 

C-^  =  2.54  -  0.01  x  10^  dyne/cm 

C2  =  0.23  t  0.01  x  10^  dyne/cm 

=  0.13  1  0.05  x  10^  dyne/cm 

=-0.12  +  0.01  x  10^  dyne/cm 


Having  thus  determined  the  force  constants,  the 

isotopic  frequencies  are  calculated  and  corrected  to  line 

130 

up  with  the  observed  frequencies  for  H0Te  .  These 
are  as  follows  below: 


cm 

^-1 

cm 

1 — 1 

1 

on  a 

3  0 

H?Te120 

2047.49 

861.13 

2133.68 

H2Te122 

2047.39 

861.39 

2133.54 

H2Te123 

2047.33 

861.02 

2133.47 

H2Te12Zj 

2047.28 

860.99 

2133.40 

H2Te125 

2047. 24 

860.95 

2133.33 

H2Te12^ 

2047. 19 

860.92 

2133.26 

H2Te128 

2047.09 

860. 85 

2133.13 

H2Te130 

2047.00 

860.79 

2133.00 
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Tellurite  ion  (TeO^  ) 

The  fundamental  vibrational  frequencies  have  been 
reported  by  Siebert  (95)  to  be  as  follows  (cm-1) 


CO 

1 

“2 

“3 

“4 

758 

364 

703 

326 

The  tellurite  ion  has  pyramidal  XY^  structure  (see  fig.  12) 
and  a  force  constant  analysis  must  take  into  account  the 
apex  angle  of  the  pyramid.  No  experimental  report  of  the 
apex  angle  for  this  ion  has  been  reported  in  the  litera¬ 
ture,  Therefore j  the  expression  from  Page  l63  of  Herzberg  (49) 


cos 


1 


4 


2  2 

V3  V4 

2 — 2 

v  v  ~ 

1  2 


3My  -  Mx 
3MY  -  Mx 


where  sin  5  -  sin  p- 

73 

and  a  =  apex  angle 
t  h 

v  =  i  vibrational  frequency 
My  =  mass  of  oxygen 


mass  of  tellurium 
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was  used  to  calculate  the  angle.  It  is  noted  from  the 
expression  that  the  angle  is  dependent  on  the  tellurium 
mass.  Therefore,  a  was  calculated  for  each  isotope  but 
the  computations  showed  a  range  in  alpha  from  101. l6°  in 
Te^^  0^  to  101.15°  in  Te^22  so  one  could  use  a  value  of 
101.15°  and  be  well  within  the  limits  of  experimental 
error.  As  to  the  accuracy  of  such  a  calculation,  the 
observed  angle  for  SeO^-  has  been  determined  to  be  100.6° 
(128)  whereas  the  calculated  value  for  ^SeO^-  was  99*7°* 

g0  _ 

Assuming  the  mean  observed  angle  to  apply  to  SeO^  , 
this  represents  an  error  of  less  than  0.8$. 


The  Wilson  F-G  Matrix  method  has  been  used  to 
determine  the  force  constants.  Following  Nakomoto  (82) 
the  F  and  G  matrix  elements  are: 

F  matrix  elements : 


For  A^  vibrations 

F11  =  fr  +  2frr 
^12  =  r(^ra  + 


F 


2 


22 


r  ( f a  +  2fnn) 


ra 
>  ' 
acr 


For  E  vibrations 

F11  =  fr  -  frr 
F12  =  r^"fra  +  ^raT ^ 
F22  ~  r 
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G  Matrix  elements: 
For  vibrations 


II 

1 — 1 

1 — 1 
o 

LL  + 

y 

(1  + 

2  cos  a)  [i 

X 

II 

CVI 

rH 

O 

2 

(1+2 

cosa) ( l -cosa) 

r 

sin  a  ^x 

G22  " 

2 

r 

1+2  cosa~|  r  ,  o  \ 

l+cosa  +  2^x(l-cosa) 

For  vibrations 

G11  =  +  ^x^1  ”  cos 

2 

r  =  -  —  (l-cosa) 

u12  r  sina  ^x 

1  2 
G00  =  — q - [(2+cosa)|j,  +  (l-cosa)  n  ] 

r  (l+cosa)  y 


where 

Te-0  bond  distance 
1 

Oxygen 

The  vibrations  consist  of  the  two  non-degenerate  vibra¬ 
tions  ^  and  ^  anc^  the  E  vibrations  consist  of  the  doubly 
degenerate  and  ^  vibrations. 
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The  secular  equation  for  both  and  E 
vibrations  is 


11F22  F12  ^ (G 


11G22  G12 


)  =  0 


There  are  six  force  constants  but  only  four  distinguish¬ 
able  frequencies  so  in  performing  the  analysis  it  has 
been  assumed  that  f  and  f  may  be  neglected  leaving 

X  LX  X  LX 

f  ,  f ,  fa,  and  f  to  be  determined.  It  is  to  be 
noted  also,  that  because  F^.  is  always  multiplied  by  G^., 
the  "r's"  cancel,  so  for  this  particular  molecular  structure, 
the  Te-0  bond  distance  need  not  be  known  if  the  apex 
angle  is  known.  Using  the  frequency  data  and  the  IBM  7040 
program  (pageA58)  the  following  set  of  force  constants 


f  =  4.377  x  10^  dyne/cm. 
f  =  0.493  x  10^  dyne/cm. 

CL 

f  =  0.321  x  10^  dyne/cm. 
rr 

f  =  O.O78  x  10^  dyne/cm. 
aa 


was  selected: 
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These  were  then  used  to  calculate  the 
isotopic  frequencies  given  below: 


OD 

1 

cm-'*' 

“2 

-1 

cm 

“3 

cm  ■*■ 

“4 

-1 

cm 

m  120 

Te 

760. 56 

366.97 

706.89 

326.79 

Te122 

760.01 

366.35 

706.07 

326.62 

Te123 

759.74 

366.04 

705.66 

326.54 

Te12^- 

759-48 

365.74 

705.26 

326.46 

Te125 

759.22 

365.44 

704.87 

326.38 

Te126 

759.97 

365.14 

704.49 

326.30 

Te128 

758.40 

364 . 56 

703.73 

326.15 

Te180 

758.00 

364.00 

703.10 

326.00 

. 
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APPENDIX  II 

IBM  7040  computer  programs  for 
calculation  of  isotopic  vibrational 
frequencies  and  partition  function 
ratios . 
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A.  Calculation  of  isotopic  vibrational  frequencies 

and  partition  function  ratios  for  diatomic  species 


In  this  program,  %eavy  and  Q)light  are  com¬ 
puted  for  diatomic  species  for  which  values  of  and 
^  x  are  available.  The  partition  function  ratio  for 
temperatures  0°  -  100°  in  steps  of  10°  and  from  100°  - 
1000°C.  in  steps  of  100°  is  then  computed  according 
to  equation  15  on  page  48. 

The  following  symbols  were  employed  for 


data  which  must  be  supplied: 


needed . 


X130, X122 


WE130(I) 


WEX130 ( I ) 


X1(I) 


=  masses  of  the  two  Y-atom  iso¬ 
topes  for  which  isotopic  data 
is  desired  (atomic  weight  units) 

=  known  value  of  co^  for  the  1^ 
diatomic  molecule  applying  to 

the  most  abundant  Y  isotope, 

m  130 

m  this  case  Te  . 

t  h 

=  known  value  of  x  for  the  I 
molecule  and  most  abundant  iso¬ 
tope  . 

t  h 

=  mass  of  the  I  X  atom  in  XY 

74  / 

diatomic  molecules,  eg.,  Ge  (awu 


In  this  work. 


the  ratio  of  Te 


130 


/Te 


122 


was 


is  the  most 


abundant 


isotope 


so  compu¬ 


tation  of  co  and  ci)  x  for  the  heavy  isotope  was  un- 

necessary.  Should  the  heavy  isotope  not  be  the  most 

abundant,  then  statements  must  be  included  computing 

a)  and  co  x  for  the  heavy  isotope, 
e  e  e 
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$  IRSYS 

%  JOR 

733003  PARTITION  FUNCTIONS  niATOMj 

$ T  IMF 

003 

$  T  R  JOB 

TE 

NODECK 

$  I  BETC 

D I  ATOM 

NOLIST 

DIMENSION  WF130(7),  WEX130(7)*  RM130(7)*  RM122(7),  Xl(7),  P(7)* 
1  S  0  R  P ( 7 )  *  WEI  2  2 ( 7 ) ♦  WEX122(7)>  W130(7)*  W122(7)*  T(20)*  UU(20)* 

2 U 130(7, 20)  *  U 1  2  2  (  7  *  2  0  )  ,  A(7*20)*  B(7*20)*  PR(7,20) 

XI 30=1 29. R4853 
XI 22=1 21 .94193 
WEI  30 (  1 )  =  3 2 3  •  4 
WF1 3 0 ( 2 )  =  2 1 1  .  « 

WE130(3)=259.3 

WEI  30 ( 4 ) =796.0 
WE  1 30 ( 5 ) =  48  0 • 4 
WE  1 30 ( 6 )  =  25 1 • 5 
WE130(7)=2250.0 
WEX130 ( 1 ) =1 .00 
WEX130 ( 2 ) =0. 12 
WE  X 1 30 ( 3 ) =0.50 
WE  XI 30 ( 4) =3.50 
WE  XI 30 ( 5 ) =1  .  16 
WE  X 1 30 ( 6 ) =0. 55 
WE  XI 30(3) =41 . 7 
XI  (  1  )  =7*a.  QA26 
XI  (  2  )  =20R  .0*t7^4 
XI  ( 3  )  =  1  1 o . 940FQ 

XI ( 4 ) =1 6. 

XI ( 5 ) =27.985825 
XI ( 6 )=129. 94853 
XI ( 7 ) = 1 .008 
DO  20  1=1*7 

RM 1 30 (  I  )  =  X 1 (  I  ) * X 1 3 0 / (XI  (  I ) +X 1 3 0 ) 

RM 1 2 2 (  I)  =  X1(  I ) *X1 22/ ( XI (  I )+Xl 2  2 ) 

SQRP  (  I  ) =P M 1 3 0 (  I  ) / R M 1 2 2 (  I  ) 

P ( i ) =SORP ( I ) **0.3 
W  E 1 2  2 (  T  )  =  p(  T)*WE130(  I  ) 

WE  X 1 2  2 (  T ) =  SQRP (  I  ) *WFXl ° 0 (  T ) 

W130( I )=WE130( I )-2.*WEXl 30 ( I ) 

20  W122( I ) =WE 1 22 ( I )-2.*WEXl 22(1) 

DO  28  1=2*11 

T ( 1 ) =0. 

28  T (  I  )  =T(  1-1 )+10. 

DO  29  T  =  1 2*20 

29  T (  I  ) =  T (  1-1 ) +1 00 • 

E  =  2 • 7 1 828 

V  =  1  .  4  3  P  o  2  8 


•3**1  T  ? 

qnt«T* 

MOT  M  n  T  P  T  ' 

( r ) «  « {  .  •  i  ♦  '  w  « ( r ) 

•  on  1  ♦  (  " )  .  r " )  <  r ,i  )  r  r  ,,  #  (  t  j  ■*:  s:  r  x?  w  »inc^n»  •  m  qso  r 

.  .  ...  .  .  r  -  '  .  ■ 

1  o .  ,  o  <;  r  -  c  r  f  ' 

. rr  r  r  i  r  r  ‘  • 

O  .  ;  s  f  c  )  r  r  r 

.  >or  -  (  •  r.  -  r 

A  •  "  A  -  f  -  )  '  ' 

ras*( ajoo  nw 
%c?ss«<T)oniw 
VM  =  (  firf  rx^ 

s  r  .o«($  )Oorx3w 

r  3  .  -  (  r  )  p  r  r  x  ”  v 

'** a  .  0  ! »  (  A  )  pry  -  v 
r  #  r  &  s  ( r  )  0 «  *■  *  3W 
r  )  rx 

.V  » r  A  r  :  #  O  c  _  (  r  )  r ' 

.  •>  r  -  (  A  )  r  < 

r  #  f  r  r  f  )  r 

•  '  q^o#oc r  =  ( n r x 

«co.f*(T>rx 

V, r=T  o^  nn 

(  op  r  x-m  I )  rx)V'prx*(  i )  r  x  =  c  i  »rp  tma 
<ssrx*><  n  rx  j \ss  x*c  i )  rx*(  i  ks  tma 

(  j  )  cc  r  mo  \  (  i  )  rp'C.Ji  ) 

3  .  *  *  f  t  )  <  a  >  -(  »  )<: 

{MOP  nW#(  T  )  r  rr  (  )  Cr  r  -J  . 

( t  )  r r  r x  >  *  ( j )  qon  *  *  {  t  )  <;<;  r  X qw 
f  MOP  fX  *,-(')  rw 

( : )  cs  rx*1**.*-  f  )  ■  s  nv.  (  *  >**  rw 

rr«sn  «* 

•o-c i >t 

•  o  r  +  (  £-!  )Ts(I  )  T  B* 

?*•*  f»!  °*  on 

.  0  M(  r-T  >T  =  (  T  )  T  rr 

PC-0  rr^.l 

ocopp^#  r  =v 
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no  21  T  =1  *  7 

DO  21  J  =  1 ♦ 20 

UU( J)=V/( T ( J ) +2  73 •  ) 

U 1 30 (  I *J)=UU( J)*W130(  I  ) 

U 1 2 2 (  I » J)=UU( J)*W122 (  I  ) 

A (  I  *  J  ) =-U 130 (  I  *  J  )  /  2  • 

B (  I  *  J ) =-U 1 22 ( I >  J  )  /  2  • 

21  PR (  I  * J )  =  ( U1 30 (  I  *  J ) * ( F**  A ( I  * J )  ) * ( 1  •  -  ( F** ( 2 • *8 ( I  *  J ) )  )  )  ) 


1  / ( ui 22 (  I 

*  J  )  * 

( E**B (I*J))*(l.-( F** ( 2.* A ( I ♦ J )  )  )  )  ) 

WR  I  TF 

(  6 

*  4  ) 

WR  I  TF 

(  6 

*  3  ) 

( SQPP ( I ) ♦  1=1*7) 

WR  I  TF 

(  6 

*5  ) 

WR  I  TF 

(  6 

♦  3  ) 

(  P  (  I  )  *  I  =  1  *  7  ) 

WRITE 

(  6 

♦  6  ) 

WR  I  TF 

(  6 

*  3  ) 

(WF122 (  I  ) *  1  =  1 *7) 

WR  ITE 

(  6 

*  7  ) 

WRITE 

(  6 

*  3  ) 

(WFX122U)*  1  =  1*7) 

WR  ITE 

(  6 

*  8  ) 

WR  ITE 

(  6 

*  3  ) 

( W 1 3  0 (  I  )  *  1  =  1,7) 

WR  ITE 

(  6 

*  9  ) 

WR  ITE 

(  6 

*  3  ) 

( W1 22 ( I ) *  1=1*7) 

WR  I  TF 

(  6 

*  1 1  ) 

no  10 

J  = 

1  *20 

WR  I  TF 

(6 

♦  12) 

( PR (  I  * J  )  *  1  =  1*7) 

3  FORMAT  (2X»(r15.P) 

4  FORMAT  (  1  RH  THIS  IS  U130/U122) 

5  FORMAT  (10H  THIS  IS  P) 

6  FORMAT  (14H  THIS  IS  WF122) 

7  FORMAT  (16H  THIS  IS  WFXF122) 

8  FORMAT  (13H  THIS  IS  W130) 

9  FORMAT  ( 1 3H  THIS  IS  W122) 

11  FORMAT  ( 4 OH  THESE  ARE  THE  PARTITION  FUNCTION  RATIOS) 

1 2  FORMAT  ( 7 ( 2X  *  FI  5.8  )  ) 

CALL  FXIT 

PNO 


$  FNTp  Y 


ni ATOM 


I,frs+IUT)\V*(LI  uu 
UU*»  (  L  •  I  )  Of  r 

= '  •  i  k  <:  r 

.  \{  ..  !  )  f U-- ( l • !  )  a 

.  l#,i  K«;ru--rL.  *  ) 

.  ■  .  *  •  .  t* i  r 


.  r  *  T  )  r  *  ■ 

♦  t  )  r  c  r  )  \  r 

•  ^  ) 

f  r ,  r  •  •  ( T  ) noo? } 

.  A  ) 

*  ) 

.  '  -  T  *  (  T  )  ) 

• 

(  0 

♦ 

* 

• 

•  A  ) 

.  .  ■  , 

*  ■'  ) 

•  d  ) 

("•1*1  •  m  6  e  x  wi 

• 

•  *  > 

IT!  qw 

(  T  «  f  « 

•  A  ) 

{ r  r 

• 

( s  r 

♦  -■>  ) 

t ru\o* ru 


'  '  •  r  •  <  c  i 

IHT  Hdf) 
2  I HT  Hfl) 

n  !  'Mt  HO  A) 


•  ’.c  5  r  ) 


,  .  *,  •  -  • 

*  M.'OO 

T  AMqOl 

TAMqo^ 


r  r 
sr 


r  4 


YOTl^T? 
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B.  Computation  of  isotopic  vibrational  frequencies  for 
XY^  molecules 

This  program  was  set  up  to  compute  isotopic 
vibrational  frequencies  according  to  the  OVFF  equations 
of  Heath  and  Linnett  (46).  The  masses  of  the  (i)  iso¬ 
topes  (awu)  for  which  vibrational  data  is  desired  are 
used  in  the  program  as  data  cards.  The  following  symbols 
apply  to  data  which  must  be  supplied: 

AK  =  force  constant  (dynes/cm) 

BK  =  force  constant  k^ 

A  =  force  constant  A 

C  =  force  constant  B/R 

YM  =  mass  of  the  Y  atom  (awu) 


A  4  1 


S  TRSYS 

$ JOR  7R500R  VTRPATIONAL  WAVFLFNGTHS  X  Y6 

SJPJOR  TFF6  NODFTK 

$  I  RFTC  OCTA  NOLTST 

DIMENSION  XM (  8  )  *  H ( 8 )  *P(8)*Q1(8)*Q2(8)* 
1W1(8)»W2(8)»CW1(8) *CW2 ( 8 ) 

AK  =  4,  85F+05 

BK=0.78F+05 

A=0,075E+05 

C=0.01 1 5E  +  05 

YM  = I 9 . 004444 

FI  1  =  ( AK-2 • *C  +  4 • *A ) 

FI  2= ( 2 ,8284*4-1  , 4 1 4 2 *r ) 

F22= ( BK  +  R  • *C  +2 , *A ) 

R  F  AD  (5*1)  N 
DO  22  I =1  *  N 
READ  (5*2)  XM( I) 

H (  I  )  =  ( F  1 1  ) * ( 1 • / YM  +  2 • / XM (  I  )  )+F22*(4,/XM(  I  )  + 1 ,/YM) 
1-5.6568*F12/XM( I ) 

P  (  I  )  =  (  Fn*E22-Fl  2**2  )  *(  1  •  /YM**2+6.  /  (  YM*XM  (  I  )  )  ) 

01 (  I )  =  . 5* ( H(  I  )  +SQP T  ( H (  I )**2-4.*P(  I )  )  ) 
02(I)=.C*(H(I ) -SORT  ( H ( I )**2-4.*P( I ) ) ) 

W 1 (  I  ) =  SOP T  (01  (  I  )/,0RPP94) 

22  W2 (  I ) =  SOR  T  ( 02 (  I  )/• 0^8894 ) 

no  r o  i  =  i ,n 

CW 1 (  I  ) =  75 1 • * W 1  (  I  ) / W 1  (°) 

R 0  fW2(  I  )=727,*W2(  I  ) /W2 ( 8 ) 


WR  TTF 

(  6*  R  ) 

WR  I  TE 

( ft*  11  ) 

(  H  (  I  )  ,  I 

=  1  *  N  ) 

WR  ITE 

(6*4) 

WRITE 

(6*11) 

(  P  (  I  )  *  I 

=  1  *N  ) 

WRITE 

(6*5) 

WR  ITE 

(6*11) 

(OKI)* 

1  =  1  *  N ) 

WR  ITE 

(6*6) 

WR  I  TF 

(6*11) 

(02(1)* 

1=1 ,N) 

WR  I  TF 

(6*7) 

WR  TTF 

(6*11) 

(  W 1  (  T  )  * 

1  =  1  *  N ) 

WR  I  TF 

(6*8) 

WR  TTF 

(6*11) 

( W2 ( I ) * 

I  =  1  *  N  ) 

WR  I  TF 

(6*9) 

WR  T  TF 

(6*11) 

(  CW1  (  I  ) 

*  T  *  1  *  N  ) 

WR  ITE 

(6*10) 

WR  I  TE 

(6*11) 

(CW2 ( I ) 

*  I  *  1  *  N  ) 

FORMAT 

(14) 

2  FORMAT  (2X*F10,5) 

R  FORMAT  (ROH  THFSF  ARE  THE  WAVELENGTH  SUMS) 

4  FORMAT  (R4H  THFSF  A R F  THF  WAVFLFNGTH  PRODUCTS) 

5  FORMAT  (24H  THFSF  ARF  THF  SUM  ROOTS) 

ft  FORMAT  (  R 1  H  T  HFS r  ARF  THE  D  T  FFrprNjTF  ROOTS) 
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7 


P 


9 

10 

11 


format  (pth  this  is  onf  sft  of  wavflfngths) 
format  ( p 7 h  this  is  thf  other  sft  of  wavflfngths) 

FORMAT  (2RH  THFSF  ARF  THE  CORRFFTFD  W1S) 

FORMAT  ( 2  8H  THFSF  ARF  THE  CORRECTED  W2S ) 

FORMAT  (2X*F14.8) 

CALL  EXIT 
END 


SENTRY  OCTA 

8 

119.94288 

121.94193 

122.94368 

123.9427° 


124.94460 

125.94420 


127.94649 

129.94853 
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C.  Computation  of  isotopic  vibrati onal  frequencies  for 
XY/j  molecules 

This  program  Is  a  solution  of  the  OVFF  equa¬ 
tions  of  Heath  and  Linnett  (45)  .  In  the  program,,  the 
mass  of  the  (i)  X  isotopes  (awu)  for  which  frequency  data 
is  desired  are  put  in  as  data  cards  to  be  read  and^ 

B  =  K 

A  =  A 

Z  =  2k/(rof 

YM  =  mass  of  the  Y  atom  (awu) 
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$ IBSYS 

$  JOB  73P00B  VIRPATTONAL  WAVFLFNGTHS  TFTR AHEORAL 

$  I R JOB  TFRP4  NOOECK 

SIBFTC  T  FTP A  NOLIST 

0  I MENS  ION  Xm(8)*H(8)*P(8)»Q1 ( 8 )  *02 ( B )  * 

1 W 1  ( 8 )  ,W2 ( 8 )  *CW1 ( P )  *  CW2 ( 8 ) 

R= 145. 0F3 
A  =  3  •  0  F  3 
Z=0.0384E+05 
W  =  400 • 0 
YM=79.916 
R  E  AO  (5*1)  N 
00  22  I  =1  *  N 
P  E  AO  (5,2)  XM (  I  ) 

H  (  I )=(B +2 .6666667 *A-1 .7333337*W)*(1 ./YM+1.3333737/XM(  I)) 

S+( 1 . /YM+2 . 6666667 /XM ( T ) )*(Z+1 . 7733737*A+3 . 7373B3*W ) 

S-( 2. 666666 7*A~1 #33 733° R*W)*( 2.666666 7/XM( I ) ) 
P(I)=U./(YM**?)+4./(XM(I)*YM))*((R+?.6666667*A~1.7B37P73*W) 
S* ( Z+l . 3333333*A+3.3333333*W ) - . 5* ( 2 • 6666667*A-1 . 333373 3*W ) **2 ) 
01  (  I  )  =  .5*( H(  I ) +SQRT  ( H (  I  )**2~^.*P(  I )  )  ) 

02 (  I )  =  .5*(H(  I ) -SORT  (H(  I  )**2-4.*P(  I )  )  ) 

W 1 (  I ) =  SQR  T  ( Q1 (  I  )/ .058894 ) 

22  W2 (  I  ) =  SQR  T  ( 02 (  I  ) / . 0 5 8 8 94 ) 

00  30  1=1  , N 

CW 1(1) =209. *W1 ( I ) /Wl ( 8 ) 

30  CW2 ( I ) =64.*W2 ( I ) /W2 ( 8 ) 


WP  I  TF 

(6*7) 

WP  I  TF 

(6*11) 

(  H(  I 

) *  1  =  1  *N ) 

WP  I  TF 

(6*4) 

WP  I  TE 

(6*11) 

(  P  (  I 

) *  I  =  1  *  N ) 

WR  I  TE 

(6*5) 

WR  I  TE 

(6*11) 

(01  ( 

I )  *  1  =  1 *N ) 

WR  I  TE 

(6*6) 

WRITE 

(6*11) 

(  02  ( 

I  )  *  I  =  1  *  N  ) 

WR  I  TE 

(6*7) 

WR  I  TE 

(6*11) 

(  Wl  ( 

I  )  *  1  =  1 *N ) 

WR  I  TE 

(6*8) 

WR  I  TE 

(6*11) 

(  W2  ( 

I  )  *  I  =  1  *  N  ) 

WP  I  TF 

(6*9) 

WP  I  TE 

(6*11) 

(  cwi 

(  i  )  *  i  =  i  ,  N  ) 

WR  I  TF 

(6*10) 

WP  I  TE 

(6*11) 

( CW? (  I )  *  t  =  1  *  N ) 

1  FORMAT  (14) 

2  FORMAT  (2X*F10.5) 

3  FORMAT  ( 3 OH  THESE  ARE  THE  WAVFLENGTH  SUMS) 

4  FORMAT  ( 34H  THESE  ARE  THE  WAVELENGTH  PROOUCTS) 

5  FORMAT  (24H  THESE  ARE  THE  SUM  ROOTS) 

6  FORMAT  ( 3 1 H  THFSF  ARE  THE  DIFFERENCE  ROOTS) 
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1  3 


( 3 i h  this  is  onf  sft  of  wavflfngths) 

( 37H  this  is  thf  other  set  of  wavelengths) 


7  FORMAT 

8  FORMAT 

9  FORMAT  (  2  8  H  ThFSF  ARE  THE  CORRFCTFD  W1S) 
10  FORMAT  ( 2  8 H  THESE  ARE  THE  CORRECTED  W2S ) 

(2X*E14.8 ) 


FORMAT 
CALL  EXIT 
END 


SENTRY 

8 


TETRA 


1 1 9*9428° 
121 .941 9? 
1 22.9436° 
323.9427° 
124.94460 
125.94420 
127.94649 
129.94853 
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D.  Computation  of  vibrational  frequencies  for  linear 
XYZ  molecules 

This  program  computes  the  isotopic  vibra¬ 
tional  frequencies  using  the  force  field  equations  and 
force  constants  of  Wentink  (138).  The  (i)  X-isotopic 
masses  for  which  vibrational  frequencies  are  desired 
are  put  into  the  program  and  read  as  data  cards.  The 
symbols  used  are  the  following: 

A  =  k-^  (dynes/cm) 

B  =  k2 

c  =  k3 

D  =  ka 
E  =  r^  (cm) 

P  =  rp 

G  =  rl2 

YM  =  mass  of  the  Y  atom  (av.ru) 

The  program  computes  k^  for  the  most  abundant 

isotope  (in  this  case  the  eighth).  If  a  known  value  of 

k  is  desired  to  be  used,  instead  of  using  D  =  A 2/DM,  put 
a 

D  =  (value  of  k  )  and  remove  statements  for  A2,  R,  and  DM. 

CL 


A47 


S I RSYS 

$  JOB 

735003  VIBRATIONAL  WAVFLFNGTHS  TRIATOMIC  LINEAR 

$  I BJOB 

TRI 

NODFCK 

$ IRFTC 

SCTE 

NOLIST 

DIMENSION  XM(8)  *  H  (  8  )  *  P  (  8  )  »Q1(8)  *02(8)  *03(8)  * 

1 W1  ( 8  )  *W2 ( 8 ) *W3 ( 8 ) *CW1 ( 8 )  *CW2 ( 8 ) *CW3 ( 8 ) 

A=704.0F3 
p=458.0F3 
0=30. OF^ 

R  =  3  3  7  , 

A2=0.0pp8B4*(p**2 ) 

F=1557.0F-11 
F=  1  904  • OF- 1 1 
G=  3461  • OF-1 1 
ZM  =  31 .98 
YM  =  1  2 • 

READ  (5*1)  N 
DO  32  I =1  *  N 
READ  (5*2)  XM (  I  ) 

DM= ( ( F**2 / XM ( p ) )+(F**2/ZM)+( G**2 / YM ) ) / ( ( F**2 ) * ( F**2 ) ) 

D=  A2/DM 

H ( I ) =A* ( 1 ./ZM+l./YM) +  P* ( 1 • /YM+1 • /XM ( I ) )-2.*C/YM 
P( I )=(ZM+YM+XM( I ) )  *  ( ( A*R ) -C**2 ) / ( ZM*YM*XM ( I ) ) 

02 (  I )  =  (  ( F  **2 / XM (  I)  )  +  ( F**2/ZM )  +  ( G**2 / YM )  )*D/(  (F**7)*(F**7)  ) 
01  (  I )  =  .5*(H(  I )+SQRT(H( I )**2-4.*P(  I  )  )  ) 

Q3(I)  =  .C*(H(I) -SORT  ( H (  I  ) **2~4,*P (  I  )  )  ) 

W 1  (  I  ) = SQR  T  ( 0 1 (  I ) / 0 • 058894 ) 

W  2 (  I ) =  SOP  T  (02 (  I  ) /O. 058894) 

32  W  3 (  I ) =SOR  T  ( 03 (  I  ) /O. 058894 ) 


DO  3  3  1=1  *  N 

CW1 (  I  ) =1 347.*W1 (  I  )  /W1 ( 8 ) 
CW 2 (  I  )  =337. *W2 (  I  ) / W2 ( 8 ) 


33  rW  3(1) 
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WR  ITF 

(  6 

*  3  ) 

WR  TTF 
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N  ) 

WR  I  TF 

(  6 
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WR  I  TF 

(  6 
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N  ) 

WR  I  TF 

(  6 
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WR  ITF 
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*14) 

(03(1 
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WR  ITF 
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*  3  )  )  \ (  ( * '  A  Sv  *0 )  +  (  MS \ r  n)  +  ((o)vx\c*»i)),Mn 

Mn\^A=o 

#r- ,  (  r  )mx\ .  r+MY\  *a+  {wy\#  r*MS\.  f  )  *A»  ff  Ji¬ 
ff  T'*  *m\  ■  *\  )  \  (  <;t  f  o*  a  )  ;  f  i  J  )  MX+MY+»»S  )  *  (  !  )  Q 
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(A  n  '  *C\ (  I)  fO  )  T  ap l = ( I  ) I M 
<  * .  )\  (  !  K0)  10^  =  (  !  )  S* 

( ‘  a  ' . 0  \  (  )  r  P  )  t  qc  =  (  I  )  ? 

At  r - [  F  l 

(  '  )  '  1  )  *•  ' l'  r  -  r  i  )  r 

< 0 )  \  (  r  *  .T‘Cr  -  (  I  )  r  *) 

f  )  r  ‘  \  in 


(  P  t  A  ) 

BTT 

( \A  • 

f -  T  t  (  T  >*’ ) 

(  A  '  •  A  ) 

(  A  «  A  ) 

(H* 

r  *  I  *  (  T  )  q  ) 

(  a  r «  ^  ) 

( ’/  •  r 

-It  r  )ro ) 

(  a  f  t* ) 

t  a  #  * ) 

'T  i  rw 

(  •  [ 

*  I t  (  !  ) If  ) 

(  A  f  *  A  ) 

f  r  t  A  ) 

(  k  •  r 

=  l  .  (  T  )FO) 

(  A  r  t  a  ) 

(  o  t  A  ) 

i  m  •  r 

“  I  t  (  I  )  *  v  ‘ ) 

(  A  f  t  A  ) 

A4<> 


WR  I  TE 

(  6*9  ) 

WR  I  TE 

(6*14) 

( W2 ( I ) *  1  =  1 ,N ) 

WRITE 

(6*10) 

WRITE 

(6*14) 

(  W3 (  I ) *  I =1 ,N  ) 

WR  I  TE 

(6*11) 

WR  ITE 

(6*14) 

( CW1  (  I )  ,  1  =  1  , 

N  ) 

WRITE 

(6*12) 

WR  I  TR 

(6*14) 

( rw2 ( I ) *  1  =  1  , 

N  ) 

WRITE 

(6*13) 

WR  ITF 

(6*14) 

( r W7 ( I )  ,  T  =  1  , 

N  ) 

WR  I  TR 

(6*15) 

WR  ITE 

(6*14) 

0 

1 

FORMAT 

(14) 

2 

FORMAT 

(2X*F10.5 ) 

3 

FORMAT 

(  30H 

thfsf  arf  the 

wavelength 

SUMS  ) 

4 

FORMAT 

(  34H 

these  are  the 

wavelength 

PRODUCTS ) 

5 

FORMAT 

(  1  1  H 

this  IS  02 ) 

6 

FORMAT 

(  24H 

THFSR  arf  the 

SUM  ROOTS) 

7 

FORMAT 

(  3  1  H 

these  are  the 

difference 

ROOTS ) 

P 

FORMAT 

(  1  1  H 

this  IS  W1 ) 

9 

FORMAT 

(  1  1  H 

THIS  IS  W2 ) 

1  0 

format 

(  1  1  H 

this  is  w*m 

1  1 

format 

(  2  R  H 

T  H  R  S  r  ADR  THR 

CORRRGTRD 

W1S) 

1  2 

rormat 

(  28H 

ThRSE  ARR  the 

CORRRCTED 

W2S  ) 

1  3 

FORMAT 

(  2  8H 

these  arf  the 

CORRRCTED 

W3S  ) 

14 

format 

(2X*E15.8 ) 

15 

FORMAT 

(  16H 

this  is  k-alpha) 

CALL  EXIT 

END 

SENTRY  SCTE 
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1 19.94288 
12 1 • 94198 
] 22 . Q4?6P 
123.94278 
124.94460 
125.94420 
127.94649 
129.94853 
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E.  Computation  of  force  constants  for  non-linear  XY^> 
molecules 

The  force  equations  used  are  those  given  by 
Glockler  and  Tung  (4l).  Only  three  equations  are  avail¬ 
able  to  determine  four  force  constants.  Therefore ,  pro¬ 
grams  were  set  up  to  determine  (l)  the  possible  range  of 
values  of  the  constant  (other  values  give  imaginary  re¬ 
sults  ),  and  (2)  the  values  of  the  other  three  constants 
for  values  of  the  permissible  range  of  . 

Program  (l)  gives  the  lower  limit  to  the  per¬ 
missible  range  of  and  program  (2)  gives  the  upper  limit. 

In  these  programs ^ 

TM  =  mass  of  the  most  abundant  isotope  to 

which  the  observed  fundamental  frequencies 
apply  (grams). 

U  =  the  reduced  mass  of  the  X-Y  atom 
ie.,  1/u  =  l/'M  x  +  1/My 

A1  =  A-j-  +  A  2 

A  2  =  AxA2 

a3  =  y 

C4  =  force  constant  (dynes/cm) 

The  observed  fundamental  vibrational  frequencies 
are  put  in  as  data  cards  to  be  read  in  order  ox,,  ay,  a>  (cm"1). 
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SIBSYS 

$ JOB  735003  FORCE  CONSTANTS  TRIATOMIC  (1) 

$  I B JOB  H?TE  NODECK 
SIBFTC  EORCON  NOLIST 
DIMENSION  WW ( 3 ) 

TM=13  0.*(  1 • / ( 6 • 02  5* ( 10.**?)  )  ) 

U=  1  •  /  (  (6.02F*(10.**3)  )*  (  ( 1 • / 1 30 • ) +1 •  )  ) 

READ  ( 5  *1  ) N 
DO  22  I =1  *  N 
22  READ  (5*2)  WW( I ) 

A1 =355. 306* ( WW (1 ) **2+WW ( 2 ) **2 ) 

A  2  =  355. 30 6**2 *( WW( 1  ) **2*WW ( 2 ) **2 ) 

A3=355.306*( WW( 3 )**2 ) 

B  =  A 1  -  A  3 

H= A2  /  ( 1  ./ ( U**3 )“!•/( T M**2 )  ) 

E  =  D/2 • 

P= ( 2.* ( U**2 ) *A3 ) /TM 

R= ( 4 • *U ) / TM 

S= ( ( U**2 ) *A3*R ) / 2. 

T  = (  ( 2.*P*R ) /8. ) - ( U*A3 )  +  ( U*B ) 

V  = (  ( R**2 ) / 8  * )-2. 

X= ( (P**2) / 8 • ) -S 

C4= ( (-T-( ( T**2 ) ~ ( 4.*X*V ) ) **0 • 5 ) / ( 2 .*V) ) 

WRITE  (6*9) 

WRITE  (6*7)  C4 

1  FORMAT  (14) 

2  FORMAT  (2X*c10.5) 

7  format  (2X*E15.8  ) 

q  FORMAT  ( 1 1 H  THIS  IS  C4 ) 

CALL  EXIT 
END 

SENTRY  FORCON 

3 

2047. 

860.79 

2133. 
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r  3  r  r 


°ot? 
■t<;h  ° I ? 


t>\  iop  *'T39t? 

(p)W'.v  T  ^  ^  T  ^ 

(((  f**,^r)*a^O.M\if)#.ni 
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.  ' 

k#  r*i  ss  on 
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( c**( p  )  )*^ce.:,;,r  =  c^ 

r  a- r  a  =  c 
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$ IBSYS 

$  JOB  735  03  FORCE  CONSTANTS  TRIATOMIC  (2) 

$  I B JOB  H 2 TE  NODECK 
$  I BFTC  FORCON  NOLIST 
DIMENSION  WW ( 3 ) 

TM= 1 ?0.*( l./(6.02^*(10.**3) ) ) 

U=l./(  ( 6.025* ( 1 0.**3 ) )*( Cl. /I  30. )+l. ) ) 

READ  (5*1  )  N 
DO  22  1=1, N 
22  READ  (5*2)  WW( I ) 

A1=355.306*(WW(1 )  **2  +  WW ( 2 ) **2 ) 

A2  =  355 • 30  6**2* ( WW ( 1 ) **2*WW ( 2 ) **2 ) 

A 3= 3 5 5. 306* ( WW ( 3 ) **2 ) 

B= A1-A3 

D=  A  2 / ( 1 • / ( U**2 ) -1 . / ( TM**2 )  ) 

E=D/2 • 

P=  (  2.*  ( IJ**2  )  *A^  )  /TM 

R  = (  4.*U ) /TM 

S= (  ( U**2  ) *  A  3  *  R ) / 2 • 

T  = (  (2.*P*R)/P.  >-(U*A3)+(U*B) 

V= ( (R**2) / 8 • ) -2 • 

X= ( ( P**2 ) / 8 . )-S 

C 4  = (  ( -T+(  ( T **2 ) - ( 4.*X*V )  ) **0. 5 ) / ( 2. *V)  )  +  ( 0. 01*1 0.**5 ) 
WRITE  (6*9) 

WRITE  (6*7)  C4 

1  FORMAT  (14) 

2  FORMAT  (2X,F10.5) 

7  FORMAT  (2X,E15.P) 

9  FORMAT  ( 1 1 H  THIS  IS  C4) 

CALL  EXIT 
FND 

SENTRY  FORCON 

3 

2047. 

860.79 

2133. 
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•  S-<  *8\ (S**S ) ) =V 
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Programs  (3)  and  (4)  then  give  the  values  of 


the  constants  C-^,  and  for  the  permissible  range  of 

Cv  C4  going  in  steps  of  0.01  x  10  from  maximum  value  to 


minimum  value.  Two  programs  are  necessary,  since  for  each 


value  of  has  values 


4 

where  H  and  G  are  defined  as  in  the  program. 

Defining  the  symbols, 

TM  =  mass  of  the  most  abundant  isotope  to 

which  the  observed  fundamental  frequencies 


apply  (grams) 

U,  Al,  A2,  A3j  04  as  previously  defined 

Cl  =  force  constant  C1  (dynes/cm) 

032  =  force  constant 

C22  =  force  constant  C0 


The  observed  fundamental  frequencies  are  again 
'  read  in  as  data  cards  in  order. 
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S TBSYS 

$ JOB  7SFQ03  FOR^F  CONSTANTS  TPIATQMir  (2) 

$ T  I MF  OOF 

$  I R JOB  H2TF  NODECK 
SIRFTC  FORCON  NOLIST 
DIMENSION  WW ( 3 ) 

TM=130.*(l./(6.025*{ 10.**3)  )  ) 

U  =  1  *  /  (  ( 6. 02  5* ( 10.**3 )  )*(  ( 1  •  / 1 30 • ) +1 •  )  ) 

READ  ( 5  ,l  )  N 
DO  22  1=1  *  N 
22  READ  (5*2)  WW ( I ) 

A1=355.306*(WW(1 ) **2+WW ( 2 ) **2 ) 

A2=355 . 306**2* ( WW ( 1 ) **2*WW( 2 ) **2 ) 
A3=3FF.706*(WW(3)**2) 

B  =  A  1 - A  3 

D=  A  2 / ( 1 • / ( U**? ) - 1 • / ( Tm**2 )  ) 

E=D/2. 

P= ( 2 .* ( U**2 ) * A3 ) /TM 

R= ( 4.*U ) /TM 

S= (  ( U**2 ) *A3*B ) / 2  • 

T= (  ( 2 • *P*R ) / 8 • >  —  ( U  *  A  3 )  +  ( U*9 ) 

V= ( ( P**2 ) / 8 • ) -2 • 

X= ( ( P**2 ) / 8 • ) -S 
C4=-0. 0^*1 0.**5 
] 0  C1=U*A3+CA 
F=B/2.-C4/U 

G=  F- ( ( U**2 ) *A3*F ) -2 • *U*F*r  4 
H= ( 2 .* ( U**2 ) *A3/Tm+4 ,*U*C4/TM ) 

C3  2= ( H- (  (H**2)-8.*G)**0.5)/4. 

C22=U* ( F+ ( 2 .*C32/TM ) ) 

WRITE  (6*0) 

WRITE  (6*7)  C4 
WRITE  (6*4) 

WRITE  (6*7)  C32 
WRITE  (6*6) 

WRITE  (6*7)  C22 
WRITE  (6*8) 

WR I T  F  (6*7)  Cl 

C4  =  C4- ( 0. 01*1 0.**F ) 


IF ( C  4 . 

G  F  .  - 1 

. 1 5*1 0 

.**5  ) 

1 

FORMAT 

(14) 

2 

format 

(  2  X  * 

F 1  0  •  5  ) 

4 

FORMAT 

(  1  1H 

THIS 

IS 

C3  ) 

6 

FORMAT 

(  1  1H 

this 

IS 

C2  ) 

7 

FORMAT 

(  2  X  * 

E15.8  ) 

8 

FORMAT 

(  1  1  H 

THIS 

IS 

Cl  ) 

9 

FORMAT 

(11H 

this 

IS 

C  4  ) 

CALL  EXIT 
FND 

SFNTbY  FORCON 
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SIBSYS 

$JOP  735003  FORCF  CONSTANTS  TRIATOMIC  (4) 

STIMF  005 

$  I B JOB  H 2 TE  NODFCK 
SIBFTC  FORCON  NOLIST 
DIMENSION  WW ( 3 ) 

TM=130.*(1./(6.025*(10.**3) ) ) 

U=1  ./(  ( 6 • 02  5* ( 1 0 . **3 ) )*( ( 1 . /130. )+l. ) ) 

R F AD  ( 5  *  1 ) N 
DO  22  1=1* N 
22  R  F  AD  (5,2)  WW (  I  ) 

A1=3C'F.^0F*(  WW  (  I  )  **2+WW  (  2  )  **2  ) 

A2=355 , 706**2* ( WW( 1 ) **2*WW( 2 ) **2 ) 

A  3  =  355 . 3  0  6* ( WW ( 3 ) **2 ) 

B  =  A 1  -  A  3 

D=A2/ ( 1 ./ < U**2 ) -1 • / ( TM**2 ) ) 

E=  D/2 • 

P= ( 2 .* ( U**2 ) *A3 ) /TM 

R  = ( 4.*U  )  /TM 

S= ( ( U**2 ) * A3*R ) /2. 

T  = (  ( 2 • *P*R ) / 8 • ) - ( U*A3 )  +  ( U*B ) 

V= ( (P**2) / 8 • ) ~2 • 

X= ( ( P**2 ) / 8 • ) -S 
C4=-0. 05*1 0.**^ 

]0  C1=U*A3+C4 
F=R/2«-C4/U 

G=F- ( ( U**2 ) *A3*F ) -2, *U*F*F4 
H= ( 2 • * ( U**2 ) *A3/TM+4.*U*C4/TM ) 

C  3  2  = ( H+ (  ( H**2 ) -8 .*G ) **0 . 5 ) /4. 

C2  2  =  U* ( F+ ( 2 • *C3  2 / TM )  ) 


WR  I  TE 

(  6*9  ) 

WR  I  TE 

(6*7) 

C4 

WRITE 

(6*4) 

WR  I  TE 

(6*7) 

C37 

WRITE 

(6*6) 

WP  I  TF 

(6*7) 

C22 

WRITE 

(6*8) 

WR  I  TE 

(6*7) 

Cl 

C4=C4- 

( 0. 01*10. **5 ) 

IF (C4. 

GF.-l 

.15*10 

.  **  5  ) 

1 

format 

(  14  ) 

2 

FORMAT 

(  2X* 

F10.5 ) 

4 

FORMAT 

(  1  1H 

this 

IS 

C3  ) 

6 

FORMAT 

(  1  1H 

this 

IS 

C2  ) 

7 

FORMAT 

(  2  X  * 

F15.8  ) 

P 

format 

(  1  1H 

this 

IS 

Cl  ) 

9 

FORMAT 

(  1  1H 

this 

IS 

C  4 ) 

CALL  ^XTT 
FND 
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F.  Computation  of  isotopic  vibrational  frequencies  and 

partition  function  ratios  for  non-linear  XY^  molecules 

Knowing  the  value  of  the  four  force  constants 

C-^  .  C^j,  this  program  computes  the  isotopic  vibrational 

frequencies  and  partition  function  ratios  for  non-linear 
XYg  molecules  according  to  the  equations  given  by  Glockler 
and  Tung  (4l). 

The  partition  function  ratios  are  given  for 
temperatures  0°  -  100°C  in  steps  of  10°  and  from  100°  - 
1000°C  in  steps  of  100°.  The  isotopic  masses  (awu)  are 
read  in  as  data  cards s  and  symbols  for  which  data  must  be 
supplied  are : 

W(l)  =  u)(l)  observed  (cm 

Cl  =  force  constants  C-^  ...  C^  (dyne/cm) 

YM  =  mass  of  the  Y  atom  (awu) 

t  h 

DEGI  =  degree  of  degeneracy  of  the  I  vibration 

U(L  etc)  =  hc*CWW(l, j)/kT(L)  where  CWW(l,J)  is  the 

t  h  t  h 

corrected  I  '  frequency  of  the  J  isotope 

for  which  the  partition  function  ratio  is 

desired 

The  partition  function  ratio  is  computed  accord¬ 
ing  to  equation l6,  page 49  and  will  be  the  ratio  for  those 
masses  included  in  the  U(l)  statements,  eg.  CWW(l,l)  will 
give  the  ratio  of  the  7th  isotope  to  the  first  (see  page  A58) 
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$  I BSYS 

SJOB  735003  WAVELENGTHS-PART I T I  ON  FUNCTIONS  XY2  NON-LINEAR 

$  I B JOB  H2TE  NODECK 
SIBFTC  PARRAT  NOLIST 

DIMENSION  T  (  2  0  )  *  W ( 3 ) *  XM(10),  UX(10)*  A(3,10),  AL(3.10)* 

1 WW (3*10)*  CWW (3*10) ♦  U ( 140) *  PR (20),  ZM(10) 

T (  1  )  =0. 

DO  2P  1=2*11 

28  T (  I  )  =  T (  I - 1  )+10. 

DO  29  1=12*20 

29  T (  I  )  =  T (  I  —  1  ) +100. 

W (  1  )  =2047. 

W( 2  )  =860.79 
W ( 3  )  =2133. 

READ  (5*1)  N 
DO  30  I  =1  *  N 

30  READ  (5*2)  ZM(  I  ) 

YM=1 .008 

Cl =2. 5425891 *10.**^ 

C2=0. 23362281*10. **5 
C  3  =  0 . 1  3090571*1  0 . **  5 
C4=-0. 12*10. **5 
DO  31  I =1  *  N 

XM (  I  ) =ZM(  I  )  / ( 6. 025*10.** 3) 

UX (  I  )  =  (  ( 6 . 025*10. **3)/YM)  +  (l./XM (I)) 

A ( 1*1 )=(UX( I )*( ( 2.*C2 )+C4+Cl ) )-(4.*C3/XM( I) ) 

A ( 2  *  I )  =  (2.*C2* ( C1+C4 )-( 4.*(C3**2 )  ) ) * (  ( UX (  I ) **2 ) “ ( 1 . /XM (  I  )**2 )  ) 

A ( 3  *  I )  =  (Cl  —  C 4 ) *UX (  I  ) 

AL(l*I)  =  (A(l*I)  +  (  ( A ( 1  ,  I  ) **? )-4.*A (2*1)  )**0.5)/2. 

AL(2*I)  =  CA(1  *  I  )  -  (  ( A ( 1  *  I )**? )-4.*A (2*1)  )**0.5 ) /?. 

31  AL ( 7 , I ) =A ( 3*  I  ) 

DO  36  K  =  1 , 3 

DO  36  L  =  1  *  N 
DO  35  1=1*3 
DO  35  J  =  1  *  N 

35  WW (  I  ,J)  =  (AL(  I  , J ) / ( 4 • * ( ( 2  2 • /7 • )**2 )*9. )  )**0. 5 

36  CWW(K,L)=(W(K)/WW(K,1) )*WW(K*L) 


WR  I  TE 

(  6 

*19) 

WR  I  TE 

(  6 

*16) 

{ WW( 1 , 

J  )  * 

J  =  1 

*  N  ) 

WR  I  TE 

(  6 

*  20  ) 

WR  I  TE 

(  6 

*16) 

( WW ( 2  * 

J  )  * 

J  =  1 

*  N  ) 

WR  I  TE 

(  6 

♦  21  ) 

WR  I  TE 

(  6 

*  16  ) 

( WW ( 7 , 

J  )  * 

J  =  1 

*N) 

WR  ITF 

(  6 

*23  ) 

WR  I  TE 

(  6 

*16) 

(CWW( 1 

*  J  ) 

*  J  = 

1  *N) 

WRITE 

(  6 

*24) 

WR  I  TE 

(  6 

*16) 

( CWW( 2 

*  J  ) 

*  J  = 

1  *  N  ) 

WR  I  TE 

(  6 

*25  ) 

WRITE 

(  6 

*16) 

( CWW ( 3 

*  J  ) 

*  J  = 

1  *N  ) 

>  3  A 
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V  = 1  .438828 
DO  37  L  =  1  *  20 
UU  =  V / ( T ( L ) +273.  ) 

U  (  L.  )  =UU*CWW(  1,7) 

U( L+20 ) =UU*CWW( 1,1) 

U ( L  +  40 ) =UU*CWW (2*7) 

U  (  L.+60  )  =UU*rWW  (7*1  ) 

U ( L  +  «0 ) =UU*f WW ( 3  >7 ) 

37  U< L+100 )=UU*CWW( 3*1 ) 

E=2. 71828 

DEG 1=1 
DEG2  =  1 
DEG3= 1 

DO  38  T  = 1  *  20 
H=-U( I ) / 2 . 

R  =  -U (  1+20 ) /2. 

C  = - U (  1+60 ) /2. 

D=-U ( 1+40 ) /2. 

F  =  -U (  I  +  100 ) / 2 • 

G=-U( 1+80 ) /2. 

P= -U ( 1+140 ) / 2 • 

0= -U ( 1+120 ) /2 . 

38  PR (  I)  =  (  (U(  I +2  0 ) * ( E**B )* < 1 .- ( F** ( -U(  I  )  )  )  )  ) 

1/ ( U ( I ) * ( E**H )*(].-( E**( -U ( 1+20 )))))) **DEG1 
2*(  (U(  1+60 )*( E**C )  *  ( 1 •-( E** ( —  U (  1+40)  )  )  )  ) 

3/ ( U(  I +4  0 ) * ( E**D ) * ( l.-(E**(-U(  1+60)  )  )  )  )  )**DEG2 
4* ( ( U ( 1+1 00 ) * ( E**F ) * ( l.-(E**(-U( 1+80) ) ) ) ) 

5  /  (  U  (  1  +  8  0  )  *  (  E**G  )  *  (  1  (  F**(  -U  (  I  +  1  00  )  )  )  )  )  )  **r>FG3 

WRITF  (6*27) 


WR  I  TE 

(6*16) 

( PR (  I  )  *  I 

=1*20) 

1 

format 

(14) 

2 

FORMAT 

( 2X  ,F1 0. 5 ) 

16 

FORMAT 

(2X*F15.8 ) 

19 

FORMAT 

(  22H 

THFSF 

ARE 

THE 

OMEGA  IS ) 

20 

FORMAT 

(  2  2  H 

these 

ARE 

THE 

OMEGA2S ) 

21 

FORMAT 

(  2  2  H 

these 

ARE 

THE 

OMEGA3S ) 

23 

FORMAT 

(  32H 

these 

ARE 

THE 

CORRECTED 

OMEGA  1 S ) 

24 

FORMAT 

(  32H 

these 

ARE 

THE 

CORRECTED 

OMEGA2S ) 

25 

FORMAT 

(  32H 

thfse 

ARE 

THE 

CORRECTED 

OMFGA3S ) 

27 

format 

(  31  H 

THFSF 

ARE 

THF 

PART  IT  ION 

RATIOS) 

GALL  EX  T  T 
FND 

SENTRY  PAPRAT 

8 

129.94853 
127.94649 
125.944 20 
124.94460 
123.94278 
122.94368 
121.94193 
119.94288 
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G.  Computation  of  force  constants,  apex  angle,  isotopi c 
vibrational  frequencies ,  and  partition  function  ratios 
for  pyramidal  XY^  molecules 

This  program  uses  the  Wilson  F-G  matrix  method 
for  finding  the  force  constants.  The  apex  angle  is  deter¬ 
mined  using  Herzberg's  formula  (see  page  A30).  Once  these 
are  known,  isotopic  vibrational  frequencies  and  thus  the 
partition  function  ratios  can  be  computed.  These  are  com¬ 
puted  for  the  range  of  temperatures  defined  in  parts  A  and 
F. 

In  the  program,  the  following  must  be  supplied: 

YM  =  mass  of  the  Y  atom  (awu) 

W(l)  =  observed  fundamental  vibration  frequencies 
in  cm”'*'  in  order 

t  h 

DEGI  =  degree  of  degeneracy  of  the  I  vibration 

Isotopic  masses  are  read  in  as  data  cards,  so 
for  example,  if  the  partition  function  ratio  of  the  molecule 
containing  the  first  isotope  to  the  molecule  containing  the 
seventh  isotope  is  desired,  the  values,  CWW(l,7),  CWW(l,l), 
CWW (2,7),  CWW ( 2 , 1 ) ,  CWW(3,7)j  CWW(3*l)j  etc.  are  used  to 
define  successive  U's.  If  the  ratio  of  the  fourth  isotope 
to  the  fifth  is  desired*  one  would  use  CWW(l,5),  CWW(l,4), 
CWW(2,5),  CWW (2,4) ,  CWW(3,5),  CWW(3>4),  etc.  in  that  order 
designating  the  U’s  from  U(l)  .  U(l40). 
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$  I PSYS 
$  JOP 
FT  IMF 
$  T  P JOP 
$ IBFTC 


T  F04 
FOPOON 
DIMENSION 

1  A  L  PHA ( 10) 

2  G 1  1  E  <  1  0  )  * 


737007  FOP f F  CONSTANTS 
00^  ,700 
NODFCK 
NOL 1ST 

W  (  4  )  >  XM ( 1 0  )  ♦  UX ( 1 0 )  » 
*  ANGLF ( 10) *  G1 1 A ( 10 ) , 


PAPTITTON  FUNCTIONS,  PYRAMIDAL  —  XY3 


AR ( 4 ) ,  BFTA ( 10 ) , 

G 1  2  A ( 1 0 )  ,  G  2  2  A ( 10) , 


Gl 2  F ( 1 0 )  ,  G22E ( 1 0 ) ,  AA(10),  RA(10),  AE(10) 


2  P 


2P 


3BE (10),  A L ( 4 , 1 0 )  ♦ 
4C WW (4,1 0 ) ,  RP ( 20 ) 
DO  2  8  1=2,11 
T (  1  )  =0. 

T (  I)=T(I-1)+10. 

DO  20  1=12,20 
T  (  I  )  =  T  (  I  - 1  [  )  + 1  0  0  . 


WW ( 4, 1 0 ) 
SR ( 20 ) , 


T ( 20 )  ,  U ( 160 )  ,  PR ( 20) , 

CO A ( 2  0 ) ,  C  OM (20),  HRA(20) 


HPM ( 20 ) 


YM= l 6* 

UY= (6,025*10.**? ) /16. 
W( 1 ) =758. 

W ( 2  )  =364. 


W( 3 ) =703. 

W( 4)=326. 

READ  (5,1)  N 
DO  30  1=1, N 
30  READ  (5,2)  XM  (  I  ) 

DO  31  1=1  ,  N 

Pi  UX ( I ) = ( 6. 02^*1 0.**7 ) /XM ( I ) 

no  32  I  =1  ,  A 

77  AB  (  I  )  =  4  •  *  (  (  2  2  .  /  7  •  )  *"*  2)*9.*(W(  I  )**2) 

C  A  =  AB ( 1  )  +  AB ( 2 ) 

DA  =  AB ( 1 ) - AB ( 2 ) 

CE  =  AB( 3  )  +  A B ( 4 ) 

DE  =  AB ( 3 ) - A  R ( 4) 

DO  33  1=1, N 
R=  1  . 

PET  A l  I  ) =  A  R  C  0  S ( 1 ,/(4.*(W(3)**2)*(W(4)**2)/((W(1  ) **2 )*( W (?)**? )  )  ) 
1+( (?.*YM)-XM( I ) )/( (7.*YM)*XM( I ) ) ) **0 • 5 
AL  PHA (  I  )  =  2.*< ARSIN(  (  (7. **0,5)/?.  ) *S  T  N ( PFT  A (  I  )  )  )  ) 

ANGLF  (  I  )  =  ALPHA  (  I  )  *57  •  20*^8 

611 A( I ) =UY+ ( ( 1 .+?.*C0S( ALPHA ( T ) ) )*UX ( T ) ) 

G12A(I)=-(2.*(1 .+2.*COS( ALPHA ( I ) ) >*( 1 . -COS ( ALPHA ( I ) ) ) *UX ( I ) ) 

1/ ( R*S I N (ALPHA ( I ) ) ) 

G2  2  A (  I  )  =  ( 2 • * ( 1 • +2 • *COS ( ALPHA (  I  )  )  )*(UY+2.*UX(  I ) 

!*( i .-COS ( ALPHA ( I )))))/ ( (R**2 >*( 1 .+COS( ALPHA ( I ) ) ) ) 

Gl IE ( I )=UY+(UX( I )*( 1 .-COS (ALPHA ( I ) ) ) ) 

Gl ?  E (  I  )=-(((  1  .-COS (ALPHA (  I )  )  ) **2 ) *UX (  I )  ) / ( R*S I N (ALPHA (  I  )  )  ) 

G2  2  E (  I ) = (  ( 2.+COS ( ALPHA (  I  )  )  )*UY+ (  ( 1 .-COS ( ALPHA (  I  ) )  ) **2 ) *UX (  I )  ) / 

1 ( ( R**2 ) * ( 1 • +COS ( ALPHA ( I ) ) ) ) 


o  , 

T?  T  ICW 

.  r  f)'  1*  •(  r>>  *(  n^K  .'A)*,  ;."M  /*v 

.  . 

•  f  r  r  >  ■  .  (  r  )  ■ a  •  <  r )  •'  *<  r  >  -  ♦  ( o  r ) *  c  ro  •  (  r )  ~  r  'os 

*  (c%)^  q  •  ( oa  r )  ♦  ( oc; )  t  ♦  ( o  i  *  a )  ww  •  ( r.  r  *  a  )  j  a  •  ( c  n  ir p 

.  )  5 ) A  CO  # ( 0  S )  •  (  )  «  ( 0  I • 

r  r  *  c  -  i  c  t; 

•  -  ( r  >  t 
.Of  +  ir-T)T«(!)T 

c  r  T  On 

.  )Cr4-(r~f)T=(T)T  oc 

.  *  r  S  M  v 

•  M\(p**.or**sc.A>*Yi 

,«3tS|  -  )  w 

.^rs(  5  )  w 

. 

•  dSCM  A  )  W 
n  (  r .  a )  ono 
V*  ♦  f  S  ! 

mi  x  ( •> )  ob ■  o 
r=  t  r  c  on 
<  T)MX\ (  r**#o  r*^<:  cA)e(T)X< 

A  •  r  =  T  r  r  n 

(  r  *  *  (  i  ■  )  .  r>#  (  ^  «  *  <  #  r  \  .  <;  *  )  )  #  .  * .  (  |  )  ^  A  c  ^ 

is )  i '.-( r )  r  a  =  * 

(  A  )  <*A  +  t  F 

Htl*i  er  oa 
•  r  =  9 

fr)  )  )  \  c  c»  *  (  4\ )  •  )  {  q**  (  r  )  w  )  >  .A )  \  .  r  )  (  T  >  *  T 

*>  .  »  {  I  (  I  )  ' >  *  I  my#  .  p  )  )  \  f  {  T  )  mx-  (MY*  .D  )  *  r 

(  (  i  (  l  !  r  -  -  )  m  r  .  c  \  (  -  .  0  * :  *  .  p  )  )  > 14  1  ,*  )  *  #  c  =  (  t  )  'no  (  a 

p3o<:#r?*n  i  ahq j a  =  ( i ) i  ioi-* a 

t  i  f  )  H  r :  t  •  )  «,'./)]■  •  f  •  )  *  r 

.  •  )  .  <  r  3 

(  )  <  '  ,C.*  )  *  <  (  (  it  <  )  ~  *  9C  »  .  r  )-.')=(  !  )  A  c  r 

r  J  '  .  (  r  *  >  fl  ))'(((((  T  )  r  c  JA  )  r  )  r 

(  (  m  t  )  *  j  a  -  .')»■(  '  x  -  v  -  (  )  r  r ' 

i ^ii \(  ( i  )xu*(  c##(  ( ( r ) ahqja ) 20'-. n  >  )-*d  nsr 

- .  n  )  f  :m  m  !  )  ■'  |A)  '  .  I 

( ( ( ( i )  *  wq  ja  )  r  r )  *  ( p**°  > ) » 


A60 


E  A  = ( 4.*CA* ( G12A ( 1 ) **2 ) /Gl 1 A ( 1 )  )- (4**CA*G22A ( 1 )  ) 

HA= ( 4.* (G22A ( 1 ) **2 ) ) - ( 4 . *G22 A ( 1 )*(G1?A<1)**2))/G11A<1) 
P A  = ( C  A  *  *  2 ) - ( DA**2 ) 

EE  = ( 4 •  *CE* ( G1 2E ( 1 ) **2 ) /Gl 1 E ( 1 )  ) - ( 4. *CE*G22E ( 1 )  ) 

HE= ( 4.* (G2  2E (  1  ) **2 )  ) - ( 4.*G22E ( ]  ) *( G12E ( 1 ) **2 )  ) /Gl IE ( 1 ) 
PE= ( CE**2 )- ( DE**2 ) 

E22A= ( -EA-( ( EA**2 ) ~4.*HA*PA ) **0.5 ) / ( 2.*HA ) 

E22E=(-EE-( ( EE**2 )~4.*HE*PE ) **0. 5 ) / ( 2.*HE) 

FI  1A=( CA-F22A*G22A C 1  )  ) /Gl 1  A ( 1  ) 

FI  1 E= ( C F-F2  2E*G22P( 1 ) ) /Gl IE ( 1 > 

FPP= (FI  1  A- FI  IE) /5. 

FP  = ( F] 1  A  +  2 • * F 1 1 E ) / 3  • 

FA  LPHA= (F22A+2.*F22E ) / ( 3.*(R**2 )  ) 
FALPAL=(F22A-F22E)/(3.*( R**2 )  ) 

A A (  I )  =  ( FI  1 A*Gi 1 A ( I )  )  +  ( F22A*G22A(  I )  ) 

BA ( I )=(F11A*F22A)*(G11A( I )*G22A( I )-G12A( I )**2) 

AE (  I )  =  (E11E*G11E(  I )  )  +  ( E22E*G22E (  I  )  ) 

33  BE ( I >=(F1 1E*F22E)*(G11E( I >*G22E( I )-G12E( I )**2 ) 

DO  34  J  =  1 * N 

AL ( 1 » J )  =  ( A A ( J )  +  (  ( A  A ( J )**2 ) ~4 • *B A ( J)  )**0.5)/2. 

AL(2*J)  =  (AA(J)-((AA(J)**2)-4.*PA(J))**0.c')/2. 

AL(3*J)  =  (AE(J)  +  (  (AE(J)**2 ) - 4 • *  B  E (J)  )**0.5)/2. 

54  AL(4*J)  =  (AE(J)-(  ( A E ( J )**2 ) ~4.*BE ( J) )**0.c )/2. 

DO  36  K=l*4 
DO  36  L  =  1 *  N 
DO  35  1=1*4 
DO  35  J  =  1  *  N 

35  WW(I,J)=(AL(I,J)/(4.*((22./7.)**2)*9.))**0.5 

36  CWW(K*L)=(W(K)/WW(K*1 ) )*WW(K,L) 

V= 1.438828 

DO  37  L  =  1 *  20 
UU  =  V/( T (L ) +2  73 •  ) 

U ( L ) =UU*CWW( 1*7) 

U( L+20 ) =UU*CWW (1*1) 

U( L  +  40 ) =UU*GWW( 2  *7 ) 

U( L  +  6  0 ) =UU*C WW (2*1  ) 

U ( L+80 ) =UU*CWW (3*7) 

U( L+100 )=UU*CWW( 3*1 ) 

U( L+120 )=UU*CWW( 4*7 ) 

37  U( L+140 )=UU*CWW( 4, 1 ) 

E=2. 71828 

DO  38  1=1*20 
A  =  -U(  I  )  / 2 • 

B=-U( 1+20  )  /2 • 

C  =  -U( 1+60  )  /2 • 

D  = - U (  1+40 ) /2 • 

F  = - U (  I +  100) /2. 

G  = - U ( I  +  80  )  /2 • 

P=-U( 1+140 ) /2. 

Q  =  — U (  1  +  ] 20  )  /2 . 
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DE  G 1  =  1 
DEG2= 1 
DEG3=2 
D  E  G  4  =  2 

PR(  I  )  =  (  ( U (  1+20 ) * ( F  ) * ( 1 .-{ E**  f-U(  I)  )  )  )  ) 

1 / ( U ( I )*(E**A)*(1.-(E**(-U(  1  +  20)  ))) ) ) **DEG1 
2* (  (U(  1  +  60 )*( E**C)*( 1 .-( E**(-U(  1+40)  )  )  )  ) 

3/ ( U  C  I+40)*(E**D)*(1.-(E**(-U(  1+60)  )))  )  )**DEG2 
4*(  (U(  I+100)*(E**F)*( l.-(E**(-U( I +80)  )  )  )  ) 

5 / ( U ( I+80)*( E**G)*( l.-(E**(-U( I +100) ) ) ) ) )**DEG3 
6* ( ( U ( 1+140 )*( E**P ) *( l.-(E**(-U( 1+120 )))) ) 

7 / ( U ( I+120)*(E**Q)*(1.-(E**(-U( 1+140) ))) ) )**DEG4 
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WR  I  T8 

(  6 

♦  16) 
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WR  I  TF 

(6*23) 

WR  I  TF 

(6*16) 

(CWW( 1 , J ) , J=1 ,N) 

WR  ITE 

(6*24) 

WR  I  TF 

(6*16) 

( CWW ( 2  *  J ) , J  =  1  ,N ) 

WR  I  TF 

( 6*25 ) 

WR  I  TF 

(6*16) 

( CWW( 3  * J ) * J=1 ,N ) 

WR  I  TF 

(6*26) 

WR  I  TF 

(6*16) 

( CWW ( 4  *  J ) * J=1  , N ) 

WRITE 

(6,27) 

WR  I  TF 

(6*16) 

(PR(I)*I=1*20) 

1 

FORMAT 

(14) 

2 

FORMAT 

(2X,F10.5  ) 

3 

FORMAT 

(13H 

this  IS  G 1 1 A ) 

4 

FORMAT 

(  1  3H 

THIS  IS  G1 2A ) 

5 

FORMAT 

(  13H 

THIS  IS  G22 A ) 

6 

FORMAT 

(  1  1H 

THIS  IS  CA) 

7 

FORMAT 

(  1  1  H 

THIS  IS  DA) 

8 

FORMAT 

(  1  3H 

THIS  IS  F22A) 

9 

FORMAT 

(  1  3H 

THIS  IS  F 1 1  A ) 

10 

FORMAT 

(  1  3H 

this  is  f22F) 

11 

format 

(  1  3H 

This  IS  FI 1 F ) 

1  2 

FORMAT 

(  1  2H 

this  is  frrj 

13 

FORMAT 

(  1  1 H 

THIS  IS  FR) 

14 

FORMAT 

(  1  5H 

this  is  falpha) 

15 

FORMAT 

(  16H 

THIS  IS  FALPALP) 

16 

FORMAT 

(2X*E14.8 ) 

17 

FORMAT 

(  14H 

this  is  alpha) 

18 

FORMAT 

(  2  5  H 

THIS  IS  ALPHA  in  DEGREES) 

19 

FORMAT 

(  22H 

THESE  ARE  THE  OMEGA  1 S ) 

20 

FORMAT 

(22H 

THESE  ARE  THE  OMFGA2S ) 

21 

FORMAT 

(22H 

THESE  ARE  THE  OMEGA 3 S ) 

22 

FORMAT 

(22H 

THESE  ARE  THE  0MEGA4S ) 

23 

FORMAT 

(  32H 

3HFSF  ARE  THE  CORRFCTFO 

OMEGA  IS ) 

24 

FORMAT 

(  32H 

THFSF  ARE  THE  CORRFCTFO 

OMFGA2S ) 

25 

FORMAT 

(  32H 

THESE  ARE  THE  CORRECTED 

OMEGA  3  S ) 

26 

FORMAT 

(  32H 

THESE  ARE  THE  CORRFCTFO 

0MEGA4S ) 

27 

FORMAT 

(  3  1 H 

THESE  ARE  THE  PARTITION 

RATIOS ) 

CALL  EXIT 

END 

SENTRY  FORCON 
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123.94278 

122.94368 

121.94193 

119.94288 
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H.  Computation  of  partition  function  ratios  for  general 
molecules 

This  program  computes  the  partition  function 
ratio  for  any  molecule  for  which  isotopic  vibrational  data 
is  known,  using  equation  16,  page  49.  The  isotopic 
vibrational  frequencies  (cm  )  are  read  in  as  data  cards. 

If  Q2/Q2  is  wanted  where  subscript  1  refers  to  the  light 
molecule,  and  2  refers  to  the  heavy,  the  data  cards  are 
ordered  as  follows:  cd  (light),  (heavy),  (light), 
(heavy),  etc. 

The  degeneracy  of  the  frequency  is  indicated 
by  DEGI  corresponding  to  to(l).  For  example,  the  degeneracy 
of  for  SCTe  is  2,  thus  DEG2  =  2  etc.  Should  there  be 
more  or  less  than  three  frequencies  per  isotopic  molecule 
in  the  partition  function  ratio,  the  statements  of  PR(l) 
must  be  altered  accordingly.  The  program  for  TeF^  (page  A65) 
is  an  example  where  two  of  the  normal  frequencies  are 
altered  by  isotopic  substitution;  whereas,  the  program  for 
TeO^  is  an  example  where  four  of  the  normal  frequencies 
are  altered  by  isotopic  substitution. 
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$  I BSYS 

5 JOB  735003  PARTITION  FUNCTION  RATIOS 

S I B JOB  SCTE  NODECK 
SIBFTC  PARRAT  NOLIST 

DIMENSION  T ( 20  )  *  W(6),  PR(20),  U(120) 

A= 1 .433829 

T ( 1 ) =0. 

DO  24  1=2*11 

24  T ( I)=T( I-l)+10. 

DO  25  1=12*20 

25  T (  I ) =  T (  I - 1 )+100. 

READ  (5,1)  N 

DO  26  1=3  *  N 

26  READ  (5,2)  W  (  I  ) 

K  =  0 

DO  28  J  =  1 *  N 
DO  27  1=1,20 
K=  K+ 1 

27  U ( K  )  =  ( A/ ( T (  I  )  +273. )  ) *W( J ) 

28  CONTINUF 
F  =  2. 71828 
DE  G 1  =  1 
DEG2=2 
DEG3= 1 

DO  29  1=1*20 
H=-U( I ) / 2 • 

B=-U( 1+20 ) /  2  • 

C  =  -U(  1  +  60  )  /2 • 

D  =  -U (  1+40  )  /2. 

F  =  -U ( 1  +  100 ) /2. 

G  =  -U(  1  +  80  )  /  2  • 

P  =  -U (  1  +  140  )  / 2 . 

0=  -U  (  1  +  120/2. 

2Q  PR (  I  )  =  (  ( U ( I  +  20>*(E**B)*(1.-(E**(-U( I  ) )  )  )  ) 

1 / ( U (  I  )*(E**H)*( l.-( E**(-U(  1+20 )))))) **DEG1 
2* (  ( U(  1+60 )*( E**C ) * ( 1 .-( E** ( ~U {  1+40 )  )  )  )  ) 

3/ ( U (  I  +  40 ) * ( E**D ) * ( l.-(E**(-U(  1+60)  )))  )  )**DEG2 
4* (  (U( I+]00)*(E**F )*( l.-(E**(-U(  1+80)  )>>  ) 

5 / ( U (  I  +  8 0 ) * ( F**G ) * ( l.~(F**(~U(  I +  100)  )  )  )  )  )**DEG 
DO  36  1=1,20 
36  WRITE  ( 6*8 )  PD ( i ) 

1  FORMAT { 14 ) 

2  FORMAT ( 2X  *  F7 . 2 ) 

8  FORMAT  (2X*F1^.8) 

CALL  EXIT 
END 

SENTRY  PARRAT 

6 

1347.4004 

1347.0 

340.0 

337.0 

426.3703 

423.0 
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$ IBSYS 

5 JOB  735003  PARTITION  FUNCTION  RATIOS 

$  I B JOB  TEF6  NODECK 
SIBFTC  PARRAT  NOLIST 

DIMENSION  T ( 2  0 )  *  W ( 6 ) *  PR(20)*  U(120) 

A= 1 .438829 
T  (  1  )  =0  • 

DO  24  1=2*11 

24  T (  I  ) =  T (  1-1  )+10. 

DO  25  1=12*20 

25  T(  I  )=T(  I  —  1 J+100. 

READ  (5*1)  N 

DO  26  I  =1  ♦  N 

26  READ  (5*2)  W ( I ) 

K=  0 

DO  2  8  J  =  1 *  N 
DO  27  1=1*20 
K  =  K  +  l 

2  7  U(K)  =  (A/(T(  I ) +273 . )  )*W(J) 

29  CO  N  T I NUF 
E=2. 71828 
DE  G 1  =  3 
DEG2  =  3 

DO  29  1=1*20 
H=-U( I ) /2. 

B=-U( 1+20 ) /  2  • 

C=-U ( 1+60 ) /2 • 

D=-U ( 1+40 ) /2 • 

F  =  - U (  1  +  100 ) /2. 

G=-U ( I +90 ) /2. 

20  PR (  i )  =  (  (U(  I  +  2 0 ) * ( F**B ) * ( 1 • - ( E** ( -U ( I  )  )  )  )  ) 

1 / ( U ( I)*(E**H)*(1 .-(E**(-U( 1+20) ))) ) )**DFG1 
2*(  <U( 1+60 )*( E**C)*( 1 • - ( E** ( ~U ( 1+40) ))) ) 

3/ ( U ( 1+40 ) * ( E**D ) * ( 1 ( E** ( — U ( 1+60 ) ) ) ) ) )**DEG2 
DO  36  1=1*20 
36  WRITE  ( 6*8 )  PR (  I  ) 

1  FORMAT (14) 

2  FORMAT ( 2X,F7. 2 ) 

9  FORMAT  (2X*E15.8) 

CALL  EXIT 
END 

SENTPY  PARRAT 
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